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(ABSTRACT) 


An experimental study was conducted to measure the effects of oxygen entrainment 
on the transport of CO in building fires, and to develop a procedure for estimating CO 
levels during a building fire. Experiments were performed with an insulated 1/4-scale 
room connected to the side of a 1/4-scale hallway forming a L-shape. Measurements of 
CO, unburned hydrocarbons (UHC), CO,, and O, concentrations and temperature were 
performed within the compartment, the hallway and post-hallway in the exhaust duct. 

The level of CO transported to remote locations from the burning room was 
hypothesized to be most significantly affected by the oxygen entrainment into the 
compartment fire gases entering the hallway. With a fixed size opening connecting the 
compartment to the hallway, the oxygen entrainment was varied by changing the depth of 
the oxygen deficient hallway upper-layer. In experiments where compartment fire gases 
entered the hallway completely surrounded by oxygen deficient combustion gases, post- 
hallway CO yields were measured to be as much as 23% greater than CO yields measured 
inside the compartment, despite the presence of external burning. With deep upper-layers 
in the hallway, geometric effects were not observed to significantly affect the transported 
level of CO. Instead, the CO level was a function of the compartment stoichiometry and 


the occurrence of external burning. 


With the compartment on the side of the hallway, the movement of combustion 
products within the hallway was measured to be non-uniform. Gases containing high CO 
concentrations (>1.9%) were measured first traveling across the hallway and then down 
the side of the hallway opposite the compartment. Only 1.0% CO was measured in the 
gases on the compartment side of the hallway. The CO concentrations were 2.0-2.3% 
when gas concentrations in the hallway reached a steady-state. This level of CO would 
result in death with less than 2 minutes of exposure. 

A procedure for predicting CO levels in building fires was also developed. The 
procedure accounts for the effects of external burning, the non-uniform transport of toxic 
gases in the hallway, the hallway upper-layer depth and the stoichiometry of the system 


on the CO levels at remote locations. 
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Chapter 1. Introduction 


CHAPTER 1. 


INTRODUCTION 


1.1 MOTIVATION 


The total number of civilian deaths in home fires throughout the United States in 
1994 dropped to an all time low of 3,425 people (Karter, 1995). Though this news is 
encouraging, the statistics for the cause of deaths in building fires reveals some less 


encouraging aspects regarding smoke inhalation deaths. 


The cause of death in building fires continues to be dominated by smoke inhalation. 
Smoke inhalation was responsible for 76% of the deaths in building fires in 1990 (Hall 
and Hardwood, 1995). From 1980 to 1990, the percentage of smoke inhalation deaths in 
structure fires has risen 1% each year. If these numbers are extrapolated, smoke 
inhalation was responsible for 80% of the deaths in structure fires in 1994. Since the 
level of carboxyhemoglobin in a victim’s blood stream (>50%) was the measure of 
whether or not death was caused by smoke inhalation, the smoke inhalation data conveys 
that an increasing percentage of deaths in building fires are the result of carbon monoxide 


poisoning. 


Two-thirds of the smoke inhalation victims were found at locations distant from the 
room of fire origin (Gann et al., 1994). This situation is demonstrated vividly in two 
unfortunate fires during the last decade, both in nursing homes. On October 5, 1989 at 
the Hillhaven Nursing Home in Norfolk, Virginia, a fire in a patient’s room resulted in 
the death of 13 people (Nelson and Tu, 1991). Each victim died of carbon monoxide 
poisoning with 12 of the victims found in a room or position down the hallway from the 
room containing the fire, see Fig. 1.1. Twenty-three (23) patients resided along the wing 


containing the burning room. Nine (9) of the victims were found in rooms on the 
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opposite side of the hallway from the burning room while only one (1) victim was found 
on the burning room side of the hallway. A similar type fire in a Southern Michigan 
hospice on December 15, 1985 claimed the lives of 8 people. Six (6) victims died of 
carbon monoxide poisoning and were found in rooms down the hallway from the burning 
room (Nelson, 1988). These incidents are examples of types of fires which claim the 


lives of hundreds of people annually. 


The statistics on CO poisoning and the above cited incidents show a clear need in the 
fire science community for models which accurately predict the transport of CO to 
locations distant from the burning room. Further insight must first be gained on the 
phenomena controlling the transport of CO to remote locations before such a model can 
be produced. With this knowledge, measurements and analysis can be performed to 
produce models enabling fire protection engineers to predict the levels of CO at remote 
locations. Engineers could use such tools to design more fire-safe buildings, from a 
smoke inhalation hazard perspective, and to prevent tragic incidents such as those 


mentioned above from continuously reoccurring. 


—_ > 2 
Yes 

q 5 
= 

=] 


X = victim of fire 


Room of Fire Origin 


Figure 1.1 The location of the victims from the fire in the Hillhaven nursing home 
(Nelson and Tu, 1991) 
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1.2 PREVIOUS WORK 


The fire science community has begun to address the problem of CO poisoning in 
building fires. Initial research in the area focused on understanding the formation of CO 
inside the burning room (compartment). More recently, research has begun to address the 


transport of CO away from the compartment. 
1.2.1 Species Production in Compartment Fires 


The “hood” experiments were the first investigations focusing on understanding the 
production of combustion gases in an idealized room fire containing a two-layer 
environment (cool air in the bottom portion of the room and hot combustion gases in the 
top portion of the room). In the experiments performed by Cetegen (1982), Beyler 
(1983), Zukoski (1985, 1991), Toner (1987) and Morehart (1990), combustion gases from 
a fire were caught inside a hood located some distance above the fire. The gases inside 
the hood simulated a compartment fire upper-layer. The air entrainment into the fire 
plume rising to upper-layer was varied by adjusting the distance between the fire source 
and the upper-layer. The global equivalence ratio (defined as the mass of the gas in the 
upper-layer from the fuel divided by the mass of the gas in the upper-layer from the air all 
normalized by the fuel stoichiometric ratio) could be varied by altering the mass loss rate 
of the fuel or the air entrainment rate into the fire plume. The global equivalence ratio 


(GER) was used to predict species yields in the upper-layer. The plume equivalence 


sal (1.1) 
Mair 


st 


ratio, 


is defined as the ratio of the mass loss rate of fuel and the air entrainment rate into the fire 


plume normalized by the stochiometric fuel-to-air ratio. This is equivalent to the global 
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equivalence ratio during the steady-state burning of the fire (assuming the fuel and air 
enter the upper-layer via the fire plume). The plume equivalence ratio is typically used to 
determine the global equivalence ratio since the global equivalence ratio is 


experimentally difficult to measure. 


Gottuk (1992a,b,c), Tewarson (1984), Pitts (1994a) and Bryner et al. (1995) used the 
GER to predict species formation in compartment fires. From experiments by Gottuk 
(1992a), it was found that the hood experiments over-predicted the CO yields at global 
equivalence ratios from 0.4 up to 1.5, see Fig. 1.2. Pitts (1994 a,b) attributed the 
difference in the CO yields to the higher upper-layer temperatures in the compartment 


experiments. 


When wood is placed on the ceiling of the compartment, CO concentrations of as 
high as 14%-dry have been measured inside the compartment. Such high concentrations 
are much greater than what would be predicted through the current GER concept (Pitts, 
1994a). Some of the situations where the GER concept has difficulties predicting the CO 
levels in compartment fires, i.e. when wood is in the upper-layer of the compartment, are 


addressed and compensated for in an algorithm developed by Pitts (1994b). 
1.2.2 Species Transport to Locations Remote from Burning Compartment 


The transport of toxic gases away from the burning room has received minimal 
attention in the past. The majority of the smoke movement studies were performed to 
measure the temperature distribution, smoke layer depth and gas temperature distribution 
in the adjacent space (e.g. Heskestad and Hill (1987)). Studies on toxic gas movement 
away from the burning room have primarily been performed to verify that high levels of 
CO are transported to a location remote (a target room) from the burning room. More 
recently, experiments have been performed to investigate the evolution of compartment 


fire gases in the adjacent space. 
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Figure 1.2 The CO yields produced in n-hexane compartment fires performed by 
Gottuk (1992a) when sampling inside the compartment (C0). Beyler’s (1983) 
(A) n-hexane hood experiment data is also shown. 
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Full-scale room-hallway experiments were conducted by Fardell et al. (1986) to 
investigate the concentrations of different species produced in an underventilated fire 
environment. A gas chromatograph was used to analyze grab bag samples of gases taken 
from both the burning room and at the end of the hallway. In the course of measuring the 
levels of nearly forty different compounds, Fardell et al. (1986) assumed CO 
concentrations to be approximately 3.5-4.0%-dry in the gases entering the hallway. 


Experiments were performed by Morikawa et al. (1990, 1993) inside a full-scale two 
level facility. A fire, fueled by furnishings typically found in a home, was present in a 
large room in the downstairs portion of the facility. Temperatures and gas concentrations 
were measured in the burn room and in an upstairs target room at different times during 
the fire. The majority of the fire experiments produced gases containing approximately 
5.5-6.0%-dry CO, but levels as high as 8.0%-dry were measured. With the door open to 
the target room upstairs, CO levels reached as high as 4.5%-dry with levels typically in 
the 3-4%-dry range. The degree of toxicity of CO compared with HCN and acrolein was 
also studied, see Fig. 1.3. The toxicity equivalency curve shown in Fig. 1.3 is based on 
the lethal concentrations of CO, HCN and acrolein for 5-10 minutes of exposure (5000 
ppm, 350 ppm and 30 ppm respectively). With the majority of the data points on the 
lower side of the toxicity curve, CO is seen in Fig. 1.3 to be equivalently or more toxic 


than both HCN and acrolein. 


Levine and Nelson (1990) performed a full-scale simulation of a townhouse fire 
which occurred in Sharon, Pennsylvania. The fire resulted in the deaths of three people. 
One of the victim’s bloodstream contained a carboxyhemaglobin level of 91%. This is 
much higher than the levels measured in most CO poisoning victims (65%). Large 
quantities of wood present in the kitchen (the room of fire origin) were found to be 
responsible for the production of extraordinarily high levels of CO. Carbon monoxide in 


excess of 8%-dry was measured in the burn room while 5%-dry CO was measured in the 
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Figure 1.3 The toxicity of CO relative to (a) HCN and (b) acrolein in gases produced in 
a building fire (Morikawa et al., 1993). 
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upstairs bedroom. The results of the study prompted researchers to perform the 


previously mentioned compartment fires with the wood ceiling. 


The events of a fire which occurred in the Southern Michigan hospice fire mentioned 
earlier were simulated experimentally by Nelson (1988). Nelson hypothesized that the 
lack of oxygen in the hallway adjacent to the burning room was the mechanism which 
allowed the transport of high levels of CO. To test the hypothesis, a full-scale experiment 
was performed with the burning room connected to a hallway containing a target room. 
The burning room contained two openings, one connecting the room to the hallway and 
one connecting the room to the ambient surroundings (situated under a fume hood). The 
hallway contained a target room, with no openings, located approximately 12 m away 
from the burning room. The air flow into the hallway was restricted to the air which 
could enter from the burn room opening. The methane gas fire in the room replaced the 
air initially in the hallway and target room with combustion gases containing low O, 
levels. Carbon monoxide concentrations were measured to be as high as 4.5%-dry only 
0.5 m above the floor in the target room approximately 200 seconds after ignition. The 
rapid cooling of the hot fire exhaust gases was also cited as being responsible for 


allowing high levels of CO to be transported to the target room. 


The studies briefly described above verified that high concentrations of CO are 
generated and transported to locations remote from the fire origin, and that CO is the 
primary toxicant formed in these situations. None of these studies, however, investigated 
the evolution of the fire exhaust gas composition as the gases are transported to remote 
locations from the burning compartment. Only one study speculated on the phenomena 
which may be responsible for the transport of high levels of CO to remote locations from 


the fire origin (Nelson, 1988). 


A study by Ewens (1994a) investigated the applicability of the compartment GER for 
the prediction of CO yields in the gases exiting the hallway, termed post-hallway CO 
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yields. The experiments performed by Ewens (1994a) contained a hallway (with 
unlimited available oxygen) connected to the compartment used in the study of Gottuk 
(1992a). In all experiments the compartment fires were underventilated and the burning 
inside the compartment extended out into the hallway, a phenomena termed external 
burning (i.e. the oxidation of fuel rich gases as they exhaust from a burning compartment 
and mix with ambient air). The presence of external burning was noted in a previous 
study by Gottuk (1992a,c) to significantly reduce CO levels downstream of the 
compartment. Ewens (1994a) found that the CO post-hallway yield was a function of 
both the fluid mechanics in the hallway and the stoichiometry of the gases entering the 
hallway. The conditions at the opening between the compartment and the hallway 
governed the hallway fluid mechanics. The fluid mechanics determined the degree of 
oxygen entrainment into the compartment fire gases entering the hallway. The air 
entrainment where the gases enter the hallway, which is mainly governed by the soffit 
size and the opening size at this location, was attributed to having the most significant 


effect on the CO levels post-hallway. 


Studies by Ellison and Turner (1959), Alpert (1972), Hinkley et al. (1984) and 
Babrauskas (1980) demonstrated that the entrainment rate into a buoyant plume is much 
higher than the entrainment into a ceiling jet. In Ewens (1994) work, no soffit at the 
hallway entrance resulted in the gases traveling the length of the hallway as a ceiling jet 
(defined as hot flowing gases traveling along the ceiling). With a 0.20 m soffit at the 
hallway entrance, a buoyant jet of fire was seen leaving the compartment. The jet was 
observed to impinge upon the ceiling approximately 0.30 m downstream of the 
compartment, after which the gases flowed down the hallway as a ceiling jet. The 
additional oxygen entrainment induced by the buoyant jet resulted in better oxidation of 


CO and lower post-hallway yields in the 0.20 soffit case. 
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Ewens (1994 a) also performed sampling along the length of the hallway, and 
measured that unburned hydrocarbons (UHC) reduced in concentration more rapidly than 
CO. This result was consistent with the findings of Westbrook and Dryer (1984) which 
showed UHC (e.g. C,H,, CH,) having a higher oxidation rate of reaction compared with 
CO for temperatures less than 1100 K. The reduction in the levels of CO was also 
inhibited by the production of CO from the oxidation of the UHC. Ewens (1994a) also 
measured a significant reduction in CO oxidation when hallway temperatures fell below 
850 K. This slow CO oxidation is consistent with modeling results (Pitts, 1992 and 
Yetter, et al., 1991). 

In the experiments of Ewens (1994a), the fire gases exit the compartment into an 
overventilated hallway upper-layer. This situation resulted in the oxidation of CO to 
levels less than 0.50% by the exit of the hallway. These levels would be lethal to humans 
if exposed to them for approximately 20 minutes. Exposure to the levels of CO at remote 
locations in experiments performed by Morikawa et al. (1993) and Nelson (1988) (4.5%) 
would cause death in approximately 1 minute. The difference in the experiments appears 
to be the amount of oxygen allowed to entrain into the compartment fire gases entering 


the hallway. 
1.3 FOCUS OF RESEARCH 


The focus of the research discussed in this document is to investigate how and when 
high levels of CO are transported to locations remote from the fire origin. An 


experimental study was conducted to determine: 
1. the levels of CO transported to remote locations, 
2. the phenomena controlling the level of CO transported to remote locations, 
3. the movement of CO away from a room on the side of a hallway, and 


4. the effects of different fuels in the compartment. 


Chapter 1. Introduction 


The experimental findings will provide data which can be utilized to develop models, and 


to predict levels of CO transported to distant locations from the fire origin. 
The level of CO transported to remote locations is mainly dependent upon 
e the chemical composition of, 
e the oxygen entrainment into, and 
e the heat losses from 


the high temperature compartment fire gases entering the hallway. For a fixed chemical 
composition of compartment fire gases, the entrainment into the fire gases is 
hypothesized to be the parameter which controls the level of CO transported to remote 
locations. The entrainment affects both the amount of oxygen entrained into the jet of 
compartment fire gases, and the heat losses from the jet. The entrainment of O, into the 
fire gases is necessary to oxidize CO. However, rapid entrainment of cooler O, 
containing gases increases the heat losses. This results in fire gas temperatures reduced 
to levels where CO oxidation cannot occur. The focus of the research was to measure the 


oxygen entrainment into fire exhaust gases from compartment fires of various GER. 


The amount of oxygen entrained into the hot (>850 K) compartment fire gases 
entering and traveling within the hallway is the parameter believed to be controlling the 
level of CO transported away from the room. The oxygen entrainment into the gases was 
varied in the study by Ewens (1994a), but never quantified. The present research 
determined the physical parameters which most significantly affect the oxygen 
entrainment into the compartment fire gases entering the hallway. The effect of varying 


these parameters on the oxygen entrainment was then quantified. 


In actual building fires, it is hypothesized that the accumulation of an oxygen 
deficient layer of combustion gases in the hallway during the beginning stages of the fire 


inhibits the amount of oxygen allowed to be entrained into the compartment fire gases 
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entering the hallway. The result is the transport of fatally high levels of CO to distant 
locations in the building. Experiments were performed to quantify the effect of varying 
the depth of the upper-layer in the hallway on the CO level transported to remote 


locations. 


With the compartment on the side of the hallway, the fluid mechanics of the hallway 
gases are highly three-dimensional. The movement of toxic combustion gases within the 
hallway is believed to be non-uniform as evidence by the location of the victims in the 
Hillhaven nursing home fire pictured in Fig. 1.1 (Nelson and Tu, 1991). To verify this 
theory, the movement of toxic gases away from the compartment was experimentally 


measured. 


To investigate the effect of different fuels being burned inside the compartment, both 
liquid n-hexane fires and polyurethane foam fires were used in the study. The effect of 
the upper-layer depth in the hallway and the evolution of the hallway species 


concentrations were experimentally determined for the fires produced by each fuel. 


A methodology was developed to estimate the level of CO transported to remote 


locations. The method considers 

1. the uniformity of gas concentration in the adjacent space, 
2. the occurrence and/or absence of external burning, 

3. the hallway upper-layer depth, 

4. the stoichiometry of the compartment and hallway, and 
5. the entrainment into the compartment fires. 


A tool was also developed to provide an estimation of combustion product 
concentrations produced inside the compartment. The tool considers only the chemical 


kinetics of the gases inside the compartment. The results of the model were compared 


Chapter 1. Introduction 


with existing experimental data. A model for estimating the oxidation of compartment 


fire gases in the hallway is also given, but not tested. 


15 


Chapter 2. Experimental Methods 


CHAPTER 2. 
EXPERIMENTAL METHODS 


2.1 INTRODUCTION 


The following chapter is devoted to describing the experimental facility and the 
methodologies used to characterize the transport of high concentrations of CO away from 
a compartment on the side of a hallway. The chapter begins with a description of the 
experimental facility. A discussion of the experimental procedures used to acquire the 
data, and the treatment of the acquired data is then presented. The final portion of the 


chapter focuses on the types of experiments performed in the study. 


2.2 EXPERIMENTAL FACILITY 


The fire experiments were conducted in the fireproof structure which consists of a 
compartment placed on the side at the end of a hallway, see Fig. 2.1. The following 
sections describe the compartment, the hallway, the fume hood and exhaust system, and 
the gas analysis system. Details of the instrumentation in each of these areas of the 
experimental facility is also given. Lastly, the details of the data acquisition system are 


discussed. 
2.2.1 Compartment 


The compartment is seen in Fig. 2.2 to be divided into two levels to separate the flow 
of air into the compartment and the flow of combustion gases out of the compartment. 
The upper level of the compartment was 1.52 m wide, 1.22 m high, and 1.22 m deep, and 
was where the fire originated. This portion of the compartment was supported by a 6.35 
mm thick, 50.8 mm steel angle iron frame with 6.35 mm thick, 50.8 mm bar stock 
supports. The floor, walls and ceiling were lined with 25.4 mm thick Fire Master, UL 
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fire rated insulation board. The fire combustion products exit the compartment through a 
variable size opening with the maximum opening size, shown in Fig. 2.2, being 0.50 m 


wide, 0.75 m high (0.375 m’). The three opening sizes utilized in the study were: 
e 0.50 m wide, 0.25 m high (0.12 m’ in area), 
e 0.50 m wide, 0.16 m high (0.08 m’ in area), and 
e 0.25 m wide, 0.16 m high (0.04 m’ in area). 


The different opening sizes were achieved by blocking off the bottom portion of the 
opening. A 0.50 m wide, 0.63 m high hatch door was located on the wall opposite the 
wall containing the opening to allow easy and safe access to the inside of the 
compartment. To allow an upper-layer to form inside the compartment, there was a 
0.20 m soffit (defined as the distance between the ceiling and the top of the opening 


where gases exit the compartment) inside the compartment. 


The lower level of the structure was a 1.52 m wide, 0.37 m high, 1.22 m deep air 
distribution plenum constructed of 3.175 mm thick steel sheet metal. Air was naturally 
drawn into the plenum during a fire through a 0.30 m diameter, 3.96 m long circular duct 
attached to the plenum on the side of the compartment opposite the opening. The plenum 
was connected to the upper level of the compartment by two 1.22 m long, 0.13 m high 
thermally shielded openings located along the floor on either side of the compartment. 
The shields extended 0.28 m inward from the side walls. 


The mass flow of air naturally drawn through the duct and into the compartment was 
controlled by attaching different diameter orifices on the inlet of the air duct (see Fig. 
2.2). The orifice diameters utilized in the study were 0.30 m (no orifice), 0.25 m, 0.20 m, 
0.15 m. The mass flow of air entrained into the compartment was determined by 
measuring the velocity of the air in the duct using a Kurz model 415, 0-2 m/s linear 


velocity hot film probe which had an accuracy of +2.5%. The probe was placed 3.0 m 
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(10 diameters) downstream of the air inlet where the orifices were attached to the duct. 
The velocity profiles were measured to be turbulent, with Re=3,425 to 9,675, for each 
orifice diameter, and were found to be quite similar in shape, see Appendix A. The mean 
velocity was determined to be approximately 44.0 mm radially from the duct wall, and 
was at the same location for each orifice size considered. The probe could, therefore, 
remain in the same location for all experiments regardless of the attached orifice 
diameter. For further information on the procedure used to determine the velocity probe 
location see Appendix A. The mean temperature of the air flowing through the duct was 
also measured. The thermocouple was 3.35 m downstream of the duct air inlet. The air 
inlet of the duct was surrounded by a wooden cubical housing to minimize the external 


disturbances on measurement of air entrainment into the compartment. 


Fires were burned in circular fuel pan positioned in the center of the compartment. 
Fuel pan diameters of 0.12, 0.15, 0.20, 0.23, and 0.28 m were utilized in the experiments. 
The fuel pan was located on a platform resting on a 15 kg range A&D load cell (1 gram 
resolution) situated in the plenum. The time derivative of the temporal fuel weight was 


used to determine the mass loss rate. 


A spark ignition system was used to ignite the liquid pool fires. The system 
consisted of a voltage source and 14 gauge steel wire leads. The 120 V-AC powered 
source contained a transformer to generate the high voltages necessary for the spark. The 
two leads enter the compartment through the wall with the ends of the wires resting off to 
the side and just below the top of the fuel pan. After connecting the leads to the source, a 
spark was generated by flipping the power switch and turning the potentiometer. After 
the pool was ignited, the power source was turned off and disconnected from the leads. 
The development of the compartment fire was observed through three polished Vycor 
glass windows installed along the height of one wall. 
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An aspirated temperature rake was present in one corner of the compartment to 
measure the temperature profile of the upper-layer gases. The rake, constructed of 
6.35 mm stainless steel tubing, consisted of eight type K (chromel-alumel) thermocouples 
vertically spaced 0.10 m apart with the top thermocouple which was 0.10 m below the 
ceiling. The rake was placed 0.10 m from either wall to avoid wall jet effects. Gases 
were drawn into the rake by a Dayton Speedaire diaphragm vacuum pump, model 
#4Z024. A water trap in an ice bath and a Gelman glass fiber filter were upstream of the 
pump to prevent water and particulates from flowing through pump. A bare 
thermocouple was located at the side of the opening to measure the temperature of the 


gases exiting the compartment. 
2.2.2 Hallway 


The combustion gases exited the compartment and entered into a 1.22 m wide, 
5.18 m long hallway, see Fig. 2.3. The height of the hallway was dependent on the 
distance between the hallway ceiling and the top of the opening, termed the inlet soffit 
height. The hallway was 1.47 m high with an inlet soffit height of 0.0 m, and 1.67 m 
high with a 0.20 m inlet soffit. The hallway was constructed of 3 mm thick, 25.4 mm 
angle iron frame. The bottom 0.91 m of the hallway walls were constructed of 16.0 mm 
thick gypsum board lined with 1.5 mm thick Fiberfax ceramic fireproof paper. The upper 
portion of the hallway was fabricated of 25.4 mm thick Fire Master, UL rated, fire 
insulation board. Two Vycor glass access windows were placed in the upper portion of 
the wall at the dead end of the hallway where the compartment was located. The burning 
of the upper-layer gases within the hallway was observed through these windows. 


Experiments were conducted with two types of blockage at the open end of the 
hallway. The open end, termed the hallway exit, was positioned underneath of the fume 
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hood. The types of experiments were: 


1. upper portion of the hallway blocked, termed an exit soffit, from 0.0 m to 0.80 m 
(see Fig. 2.1), and 


2. bottom portion of the hallway blocked leaving a 1.22 m wide, 0.20 m high 
opening at the top of the hallway exit (see Fig. 2.3). 


The upper portion of the hallway was blocked to control the upper-layer depth inside the 
hallway. The bottom portion of the hallway was blocked to highly restrict the amount of 


air entrainment into the hallway. 


The automated sampling cart shown in Fig. 2.4 allowed the measurement of gas 
concentrations and temperature profiles at a variety of locations within the hallway during 
a single experiment. The cart contained a gas sampling probe, a vertical aspirated 
temperature rake with 9 type K thermocouples and two traversing assemblies. The cart 
was manually moved along the length of the hallway on a track. The two traversing 
assemblies, each consisting of a stepper motor, a lead screw and a driver nut, were used to 
move the gas sampling probe in either a vertical or horizontal direction via a computer. 
The vertical range of motion was 0.67 m while the horizontal range of motion is 0.76 m. 
The thermocouple rake was fixed in the vertical direction, and was only allowed to 
traverse in the horizontal direction. The wires from the motors, the Teflon tubing line 
necessary for the aspiration of the thermocouple rake, and the heated Teflon line used to 
transfer the gases from the hallway to the gas analyzers were all fed out of the hallway 
through a 0.10 m diameter flexible line feed-through duct. 


The schematic of the control system for the 4.8 Amp stepper motors is seen in Fig. 
2.5. A 486 DX66 personal computer containing an internal DT 2801-A Data Translation 
board controlled the movement of sample probe and thermocouple rake. An external DT 


707 screw terminal board interfaced the DT 2801 with the two Servo Systems CMD-50 
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Figure 2.5 A schematic of the electronics involved in computer controlling the motors 
on the automated sampling cart 
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controllers, each of which was connected to a stepper motor. A FORTRAN program was 
written to control the movement of the two stepper motors. The program communicated 
with the data acquisition board using PCLAB subroutines. The user entered the distance 
the sample probe should move and the time it should remain at that location into a data 
file. This enabled the user to map out an entire sampling grid before the beginning 


experiment. The program can be activated at any time during the experiment. 


The gas sampling probe was constructed of 6.35 mm diameter stainless steel tubing 
which had a 90 degree bend on one end forming an upside down L when mounted to the 


cart. Further information on the gas sampling system is given in section 2.2.4. 


The aspirated temperature rake contained 9 type K thermocouples mounted inside 
6.35 mm diameter stainless steel tubing. The top three thermocouples were 0.05 m apart 
with the top thermocouple 0.05 m below the ceiling to measure the large temperature 
variations near the ceiling. The remaining thermocouples were 0.10 m apart since more 
gradual temperature changes were expected this far below the ceiling. The thermocouple 
beads were recessed approximately 6.0 mm from the end of the aspirating tubes to shield 
the beads from radiation. To minimize the temperature effects of the aspirator tubes on 
the gases, the gas velocity in the tubes was greater than 5 m/s. With 5.5 m of 63.5 cm 
diameter Teflon tubing connecting the rake to the pump, a 1/3 hp rotary pump was used 
to achieve the desired velocities in the aspirator tubes. A water trap submerged in an ice 
bath and a Balston 915A, DX grade, glass fiber filter kept water and large particulates 
from flowing through the pump. 


Located approximately 1.0 m outside of the end of the hallway was a video camera 
which was used to make visual record of each experiment. The camera recorded the 
experiment from the view point of a person standing at the end of the hallway under the 


fume hood looking down (upstream fluid mechanically) the hallway. 
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2.2.3 Fume Hood and Exhaust System 


A 1.5 m by 1.5 m fume hood was situated directly over the exit of the hallway to 
collect the escaping combustion gas. As shown in Fig. 2.6, the hood was connected to a 
0.46 m diameter duct, termed the exhaust duct. The duct contained an orifice plate, a 
probe for downstream gas sampling, and a laser extinction system for the measurement of 
soot. A 142 m’/min (5000 cfm) blower connected to the exhaust duct draws the 


combustion gases into the hood. 


The 0.3048 m inside diameter, sharp edged orifice plate was designed in accordance 
with the A.S.T.M. standard for the measurement of the mass flow in a duct (Fluid Meters, 
1971). The flow rate of air within the duct was determined by accessing the D and D/2 


pressure taps. 


Approximately 5 diameters upstream of the orifice plate, a 0.15 m long, 6.35 mm 
diameter sampling probe was positioned in the exhaust duct, see Fig. 2.6. To obtain a 
sample of the average concentration of gases flowing within the duct, numerous 3.2 mm 
holes were drilled along the length of the probe. The transport of the sample gases to the 


analyzers and the analysis of the sample gases is discussed in section 2.2.4. 


A laser extinction system was present 3.7 diameters upstream of the orifice plate. 
The system consisted of a 5 mW, 670 nm visible red diode laser, an Oriel neutral density 
filter, an Oriel diffuser and a photodiode. Two 9.0 mm holes were made on opposite 
sides of the exhaust duct to allow the beam of light to pass through the gases flowing in 
the duct. After the light passed through the duct, the light was detected by the 
photodiode. The filter and diffuser were placed in front of the photodiode to generate a 


sensitive signal in the desired range. 
2.2.4 Gas Sampling System 


The piping system shown in Fig. 2.7 allowed gases to be drawn from the exhaust 
duct, the compartment, or the hallway. Gases were sampled from the facility through 
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Figure 2.7 A schematic of the piping system used to sample combustion gases from 
different locations in the facility. Also shown are the CO, CO, and O, 
analyzers. 
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6.35 mm tubing by a suction produced from two Thomas 2107 CA18 diaphragm vacuum 
pumps. All heated lines and filters were insulated with approximately 10 mm thick fiber 
glass insulation, and were maintained at approximately 120°C using either Thermolyne 
electric resistant heating tapes or Hotwatt 360 W single glasrope resistance heaters. Lines 


were heated to prevent higher order hydrocarbons and water from condensing in the lines. 


Gases sampled from the probe in the exhaust duct travel through a heated stainless 
steel line. A heated Gelman glass fiber filter was used to trap the soot particulates from 


entering the gas analysis equipment. 


When gases were sampled from either the hallway or the upper-layer of the 
compartment, they were drawn out of the facility through a heated 5.5 m long Teflon line. 
The flexible Teflon line used allowed the sampling cart to move freely along the length of 
the hallway. The line was thermally shielded from radiant heat given off by the hallway 
external burning by placing the line inside of a flexible metal duct. The heavy soot 
content of the gases sampled from the hallway and compartment required the gases to be 
filtered twice to ensure the removal of soot particles from the sample. The Teflon line 
was connected to a heated Balston model #915A filter which contained a DX rated, 
borosilica glass fiber filter rated for the retention of 93% of particulates larger than 0.1 
microns in diameter. The gases then proceeded through a heated stainless steel line 
which led to a heated Gelman glass fiber filter to retain any soot particulates remaining in 


the gases. 


2.2.5 Gas Analysis 


Once the gases had been filtered of soot, they branched off into two directions. A 
portion of the sampled gases remained wet and was drawn through a heated line by one of 
the vacuum pumps. The remaining gases were dried by being drawn through a water trap 
immersed in a Fisher model #910 refrigerated circulating bath held at -10°C. The 


unburned hydrocarbon (UHC) levels were measured using the wet gases with a flame 
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ionization detector (FID). The dry gases were passed through two non-dispersive infrared 
analyzers to measure CO and CO, concentrations, and through a paramagnetic analyzer to 


measure the O, concentration. 


2.2.5.1 Unbumed Hydrocarbon Measurements 


The portion of the sampled gases which remained wet was drawn into an oven 
(maintained at 105°C) to measure the UHC wet concentration within the sample, see 
Fig. 2.8. Gases remain wet to prevent the loss of higher order hydrocarbons during the 
water trapping process. The measurement was performed using a model #12-800 Gow- 
Mac FID connected to a model #40-900 Gow-Mac electrometer for conditioning of the 
FID signal. The FID contained a hydrogen flame which was produced using a 40% 
hydrogen, 60% helium mixture as the fuel and purified air as the oxidant. Using 
Matheson rotameters, 270 cc/min of fuel and 470 cc/min of oxidant were continuously 


flowing to the FID. 


As the sample gases entered the oven, they pass through a Fairchild 0-2 psig pressure 
regulator, model #10. For an accurate measurement of the sample, the pressure in the line 
was held at 1.5 psig using a flow meter and a pressure gauge downstream of the oven. 
The 0-10 SLPM Matheson rotameter controlled the bypass flow rate (typically ranging 
from 4.5-5.5 SLPM) to maintain a 1.5 psig reading on the Magnehelic 0-6 psig pressure 
gauge. The bypass flow rate also reduced the delay time of the sampling gas system to 
approximately 10 seconds. A portion of the gases flowing through the oven were 
directed through a capillary tube which led to the FID. The capillary tube was designed 
around the recommended flow rate of 20 to 40 cc/min. The products of the hydrogen 
flame were vented out of the oven through an opening in the top of the oven. Before the 
gases were vented to the ambient surroundings, the gases passed through a hydrogen 
detection system. If the flame was extinguished and high levels of hydrogen begin to be 
vented to the atmosphere, the hydrogen detector would automatically shut off the flow of 


29 


Chapter 2. Experimental Methods 


Cooling 


To Atmosphere Coils 


——$ 
-10°C 
Coolant 


40% H, 


AIVe 


60% He 


Rotameter PA FID 


Pressure Gauge Capillary 
Tube 


Rotameter 


Rotameter 


Lines Ye 
Pump Pressure Regulator 
Heated Sample 
Line Gas In Oven at 105°C 


Figure 2.8 A schematic of the equipment and piping in the UHC analysis system. 
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hydrogen using a solenoid valve. Cooling coils at -10°C prevented the temperature 
around the hydrogen detector from exceeding the maximum operating condition of 49°C. 
The coils were cooled using -10°C coolant from a Fisher model #910 refrigerated 


circulating bath. 


The FID signal was conditioned using the electrometer which contained four possible 
sensitivity ranges; three were commonly used in the experiments. Before the start of each 
experiment, the FID was calibrated by determining the zero point and the span point for 
the desired sensitivity range. At each of the ranges, the zero point was determined by 
passing nitrogen through the FID. The span point for a particular sensitivity range was 
generated by passing gases containing different concentrations of ethylene (C,H,) 
through the FID. Since the higher order hydrocarbons tend to pyrolyze to ethylene, using 
ethylene as a calibration gas for the measurement of UHC is a valid assumption 
(Westbrook and Dryer, 1984). Ethylene concentrations of 4.71%, 5456 ppm, 615 ppm 
were used to generate span points for each of the three sensitivity points. The lower 
ethylene concentration of 615 ppm was the range of UHC expected in the exhaust, and 
was used to calibrate the FID when sampling gases in the exhaust duct. The other two 


ranges were utilized when sampling gases in the compartment and hallway. 


2.2.5.2 CO, CO, and O, Measurements 


The portion of gases which was not drawn to the FID was dried to prevent the cells 
within the CO, CO, and O, analyzers from being damaged by condensation. After the 
gases were dried, they passed through the vacuum pump and were divided again. Part of 
the gases were directed to the analyzers while the remaining gases were vented to the 
surroundings (bypass flow). The bypass flow enabled the gases to be drawn into the 
analyzers at a faster rate allowing for sampling delay times of less than 10 seconds. The 
flow rate of gases into each of the three instruments was monitored at 1 1/min using 3 


separate 0-3.5 SLPM range Matheson rotameters. 
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The dry CO and CO, concentrations were measured using Rosemont Analytical 
model 880 NDIR analyzers. The 0-5 volt output signal was linearized with respect to the 
gas concentration using a linearizer contained within the analyzer. The linearizer 
required a zero point and a span point for calibration of the instrument. The zero point 
was set by passing pure nitrogen through the analyzers. The span point was determined 
by passing a calibration gas containing a concentration equivalent to approximately 90% 
of the upper bound of the desired range. The CO analyzer ranges typically used in the 
experiments were 1000 ppm, 1% and 10%. The CO, analyzer ranges used in the 
experiments were 5000 ppm, 2%, 15% and 20%. The range was selected such that the 
concentration of the CO or CO, within the sample of the same order of magnitude and 
never exceeding the upper bound of the range. When sampling diluted exhaust duct 
gases, the lower concentration ranges (0 to 1000 ppm and 1% for CO and 0 to 5000 ppm 
and 2% CO.) were used in the experiments. The higher ranges were used when sampling 


inside the hallway and compartment. 


The dry concentration of O, was measured using a Siemens paramagnetic Oxymat 
SE analyzer. The output from the instrument was 4-20 mA. This was converted to 1-5 
volts using a 250 ohm resistor to enable the signal to be recorded in the computer. The 
analyzer operated in the 0 to 22% range, since concentrations of O, less than levels 
observed in the ambient surroundings (approximately 21%) were expected. The 
instrument was calibrated by first determining the zero point, and then setting the span 
point. Purified nitrogen was passed through the cell to determine the zero point. The 
span point was set using a gas mixture containing 4.75% O,. Since the majority of 
sample gases where accurate O, measurements were necessary contained low 
concentrations of O,, a low concentration calibration gas was used to calibrate the 


analyzer. 
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2.2.6 Data Acquisition 


A 386 computer containing various Data Translation data acquisition boards was 
used to collect the data. The computer was equipped with three internal DT 2801-A 
digital-to-analog data acquisition boards containing 12-bit resolution. The internal boards 
were connected to three external boards, one DT 707 screw terminal board and two DT 
756-Y amplified boards. The DT 707 was utilized to measure the differential outputs 
from the UHC, CO, CO,, and O, analyzers, the load cell, the air duct velocity probe and 
the laser extinction system. One of the DT 756-Y boards was used to process the signals 
of the thermocouples in the compartment rake, the air inlet duct, the exhaust duct and the 
opening between the compartment and hallway. The other DT 756-Y board was used to 
process signals of the thermocouples in the hallway rake. The two DT 756-Y boards 
were amplified, multiplexing terminal boards containing cold junction compensation for 
thermocouple measurements. The boards were used to amplify the low voltage signals 


given by the thermocouples allowing the thermocouple signals to be read directly. 


The data acquisition was controlled using a BASIC program. The program accessed 
PCLAB subroutines to acquire the data from the different measuring devices listed above. 
Before the acquisition began, the user was required to enter a file name where the 
computer will store the raw data. The data acquisition was started approximately 30 
seconds before the ignition of the fire in the compartment. After the data acquisition 
program was activated, the program began to store data in the raw data file every 2 
seconds. The raw data was determined by averaging the data samples sent to the 
computer 1 second before and 1 second after the recording time. There were typically 
60+2 samples of data which were averaged to generate a single raw data point for each of 
the 28 signals. The thermocouple signals were converted to degrees Celsius using a 
PCLAB subroutine and then stored in the raw data file. The data acquisition program 
was permitted to run 30 seconds after the compartment fire had extinguished to allow the 
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gases sampled at the end of the experiment to be analyzed. The reduction of the raw data 


file is described in section 2.3. 


2.3 DATA REDUCTION 


Once the numerical and visual data were collected, they needed to be reduced to a 
useful form. The basic data reduction converted the raw data to a useful form, in addition 
to performing other calculations. After the basic reduction, the temporal data was 
averaged over one or several time periods. The output from the basic reduction program 
was also used in other programs to calculate parameters such as the ignition index. 
Finally, visual data was analyzed in certain cases where external burning was particularly 
interesting. The remainder of the section describes in detail the basic data reduction, the 


data averaging, the calculation of the ignition index and the analysis of visual data. 


2.3.1 Basic Data Reduction 


A FORTRAN program titled FJIRERED3.FOR, see Appendix B, was written to 
convert the raw data file into seven data files. Six of the output files contain columns of 
the temporal data for the calculated and measured values, while the seventh file contains 
descriptive information about the testing conditions. The treatment of the raw data and 


the methods used for calculating various parameters are described below. 


2.3.1.1 Temperature Data 


The temperature data was stored in the raw data file in units of degrees Celsius. The 
reduction program converted the raw temperature data to Kelvin. The temporal 
temperature data from the compartment, air inlet duct and exhaust duct were stored in one 


of the output files, while the hallway temperature data was stored in a separate output file. 
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2.3.1.2 Species Concentrations 


The temporal CO, CO,, O, and UHC concentrations were stored as voltages within 
the raw data file. The voltages associated with CO, CO, and O, were converted to dry 
concentrations by entering the sampling range used for each species during the 
experiment. The UHC voltages were converted to wet concentrations by entering the 
sensitivity range of the electrometer (either 1 x 10°, 107°, or 107! with 107” being the 
most sensitive range) used during the experiment and the voltage output of the FID with 


the span gas flowing through it. 


The wet concentrations of CO, CO,, and O, were calculated in the program using the 
dry concentrations and an assumption about the amount of H,O within the sample gases. 
It was assumed that all of the H,O in the sample gas was produced from the combustion, 
with a stoichiometric proportion of H,O and CO, being produced. With hexane as the 
fuel, 6 moles of CO, and 7 moles of H,O are produced when stoichiometric burning 
occurs. The dry concentration of a species 7 were converted to the wet concentration of 


species j, C;,. using the equation 


De 
Ge 00 eel hl (2.1) 
: |, | 


having knowledge of the dry mole fraction of species j (Xz) and the dry mole fraction of 


CO, (X¢o,, ). Equation 2.1 becomes less accurate in very fuel rich environments, global 


equivalence ratio greater than 3.0. 


A lag time exists between the data measured instantaneously (i.e. temperature) and 
the species concentration measurements. The lag time of the concentration measurements 
depended on the location of the sampling (hallway, compartment or exhaust duct) and the 
instrument measuring time. The lag time associated with sampling in the hallway and the 


compartment was measured to be 8-10 seconds by injecting CO, directly into the 
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sampling probe. Knowing the time in which the CO, injection occurred, the lag time was 
determined by subtracting the time which the CO, was sensed by the analyzer from the 
injection time. The hallway and compartment lag time was the same since the same 


sampling line was used to sample from these locations. 


The lag time of the species concentrations sampled from the exhaust duct was 
defined slightly differently. When sampling in the exhaust duct, the measured species 
concentrations were aligned with the time of the fuel mass loss rate which produced the 
species. The concentrations and fuel vaporization were aligned so that the post-hallway 
yield calculations were performed accurately. The lag time was estimated by subtracting 
the ignition time of the fire from the time when CO, levels produced by the fire were first 
measured in the exhaust duct. The ignition time of the fire was estimated as the time in 
which a jump in the temperature of the top thermocouple in the compartment rake 
occurred. This lag time was sensitive to the residence time of the gases in the hallway 
which was governed by the blockage at the end of the hallway. The lag time ranged from 
20-22 seconds with no blockage at the end of the hallway to 30-32 seconds with a 0.80 m 
high exit soffit. 


2.3.1.3 Species Yields 


In the past, Beyler (1986a,b) and Tewarson (1984) used mass based species yields to 
quantify the production and consumption of species in an upper-layer fire environment. 


In general, the mass based species yield, 


Yin (2.2) 


is defined as the mass flow rate of species j (m,;) divided by the fuel mass loss rate (m,,,7). 
Except for O., all species yields calculated represented the production of species relative 
to the fuel mass loss rate. In the case of O,, the yield represented the consumption of O, 


relative to the fuel mass loss rate. Species yields were calculated when sampling was 
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performed in the compartment and in the exhaust duct where the mass flow of the gases 


in and out of the system was well defined. 
The species yields in the compartment upper-layer, 


KX jy MW, (tt pyet + Mair) 
TM fei MW, 


Ve 


Jc 


(2.3) 


were calculated using the wet species concentration of interest, the air entrainment rate 
into the compartment and the fuel mass loss rate. The molecular weight of the upper- 
layer gases, MWy,, was assumed to be equivalent to that of air. This assumption is 


accurate to within +5%. 
The post-hallway species yields (the species yield of the gases exiting the hallway), 


y. ea X ; wMW jh, 
J,ex 3 2 
M fel 


(2.4) 
were calculated using measurements of species concentrations in the exhaust duct and 
with the knowledge of the molar flow rate in the exhaust duct, 7,, . The molar flow rate 


through the duct, 


ex 


Ss 0 


P.MW,, } [mol/s] (2.5) 


was determined by applying the ideal gas law, inserting the measured exhaust duct 
temperature, 7,,, and assuming the recorded ambient pressure, P,, was equivalent to the 


exhaust duct pressure, P,.. The general equation to calculate the volumetric flow rate 


through a duct, Q,,, 


2 
2, =0.550 os 5 We [tes] (2.6) 


J1-p* 
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was taken from 1 : ir icati (1971). After determining 


the flow constants (c, y and F,) and the ratio between the orifice diameter and the duct 


diameter, the volumetric flow rate through the duct (in SI units) was calculated from the 


Q,, = 1.2099 pe [m?/s] (2.7) 


With the pressure drop across the orifice plate, Ap, measured to be 14.3 mm Hg, the 


following equation: 


volumetric flow rate through the duct was determined using the temperature of the gases 
flowing through the exhaust duct, 7,,, and the ambient pressure, P,. Since the exhaust 
duct gases were highly diluted (approximately 30:1), the assumption of the properties 
(MW, p, and C) of the gases flowing through the exhaust duct being equivalent to that of 


air is well supported. 
2.3.1.4 Smoke Extinction Coefficient and Yield 


Downstream of the hallway in the exhaust duct, the smoke level was determined with 
the use of a laser extinction system previously described in section 2.2.3. The reference 
intensity, J,, of light was measured by averaging the initial 10 seconds of the photodiode 
output when no combustion gases were present inside the exhaust duct. The signal 
produced by the photodiode, J, is proportional to the amount of attenuated laser light 
which impinged upon it. 


The extinction coefficient, 


o =(1/1,)in(1/I,), [1/m] (2.8) 


is defined as the natural log of the attenuated intensity, J, divided by the reference 
intensity, J,, all divided by the optical path length, /,. The optical path length is the 
distance between the laser and the photodiode and was taken to be 0.457 m. The mass 
optical density, 
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1 7 
Daas = —(%, 303) [m/e] 2.9) 


was calculated using the fuel mass loss rate, m,,,,, the volumetric flow rate through the 
duct, Q,,, and the extinction coefficient, o (Tewarson, 1988). 


To account for the size of the smoke particles attenuating the laser light passing 


through the duct, the specific extinction coefficient, 


- BE 
MP sk 


é [g/m] (2.10) 


was estimated using the correlation for overventilated fires. This was used because no 
correlation for underventilated fires was available (Newman and Steciak, 1987). The 
wave length of laser beam, A, was 0.670 pm, while the smoke density, p,.,, was 


estimated from experimental data to be 1.1 g/cm” (Newman and Steciak, 1987). 


Therefore, the smoke yield, 


ee Dns : (2.11) 


was calculated by dividing the optical mass density by the specific extinction coefficient 
(Tewarson, 1988). 


2.3.1.5 Fire Size 


The size of the fire was estimated as the ideal heat release of the mass loss rate of 


fuel. The ideal heat release, 
O= MyAH., [kW] (2.12) 


was determined by multiplying the mass loss rate of fuel by the lower heating value of 


the heat of combustion, AH., of the fuel. 
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The fuel mass loss rate was determined using the temporal fuel weight data. The fuel 
weight data was stored in the raw data file as 1-5 volt signal (which corresponded to a 0- 
10 kg weight), and was converted to units of mass, [kg], using the data reduction 
program. The temporal derivative of the fuel weight data was used to determine the fuel 
mass. loss rate. The temporal derivative was determined numerically over points +10 


seconds around the point of interest. 


2.3.1.6 Compartment Global Equivalence Ratio 


The degree to which the fire inside the compartment was ventilated was defined by 
the equivalence ratio. The global equivalence ratio, $,, is defined as the mass of gas in 
the upper-layer from the fuel divided by the mass of gas in the upper-layer from the air all 
normalized by the stoichiometric fuel-to-air ratio. The stoichiometric fuel-to-air ratio is 


defined as the ratio of fuel-to-air which results in only CO, and H,O as products. The 


fev | 
a (mie) 


st 


plume equivalence ratio, 


(2.13) 


is the fuel mass loss rate divided by the air entrainment into the compartment, again, 


normalized by the stoichiometric fuel-to-air ratio. 


Species production inside the compartment has previously been correlated to the 
global equivalence ratio, but this parameter is difficult to directly measure with any 
accuracy. The plume equivalence ratio of the compartment fire was easily calculated 
using the air entrainment rate into the compartment and the fuel mass loss rate. When the 
fire was fully developed (during the post-flashover period), a steady-state was reached 
inside the compartment. During the steady-state period, the fuel and the air entrained into 
the compartment was all assumed to enter the upper-layer through the plume. The values 
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of $, and $,, therefore, were equivalent to one another during this period. The value of $, 
is known during the steady-state time of the fire, and will be referred to herein as simply 
o. 

The air entrainment rate into the compartment was determined by measuring the 
mean velocity (see section 2.2.1) and temperature of the air flowing through the air inlet 
duct. The density of the air flowing through the duct was determined using the ideal gas 
law. The air entrainment into the compartment was, therefore, calculated using the mean 
velocity, the density of the gases and the cross sectional area of the duct. Refer to the 


previous section, 2.3.1.5, for details on how the fuel mass loss rate was determined. 


2.3.1.7 Residence Times 
The residence time is defined as the amount of time a unit volume of gases remains 
within the area of interest. The residence time inside the compartment, 


V, 
t UL 


Sage tle ne 2.14 
bay O\ Ten | fa) Bs 


calculates the amount of time a unit volume of air (temperature corrected) remained in the 
upper-layer of the compartment (which was assumed to be 0.60 m deep) (Gottuk et al., 
1992). 


2.3.1.8 Compartment Fire Growth and Steady-State 


The growth of a compartment fire has previously been defined using a number of 
measured variables (e.g. upper-layer temperature, fuel mass loss rate, equivalence ratio 
etc.). Since the growth of the fire is mainly governed by the fuel mass loss rate, the fire 
growth was defined by Gottuk et al. (1992c) as, 


wie 
dt 
G=t 


res,c : 2 


(2.15) 
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the residence time of the gases in the compartment multiplied by temporal change in the 
fuel mass loss rate divided by the mass loss rate. The parameter G is termed the fire 
growth parameter. The term multiplying the residence time of the gases in the 
compartment defines the rate of growth of the fire. When the growth rate of the fire is 
small compared with the residence time of the compartment gases, G goes to zero and the 


compartment is considered to be in a steady-state. 


A plot of the fire growth parameter with time, Fig. 2.9, shows the four stages of a 
typical compartment fire. The compartment temperature is often used to characterize fire 
growth inside a compartment. For comparison purposes, the temperature 0.05 m below 
the compartment ceiling is also shown in Fig. 2.9 with the fire growth parameter since. 
After the fire was ignited, a gradual growth of the fire (termed the pre-flashover period) 
was observed to occur within the compartment. As the fire became oxygen deficient 
(o>1) and the fire growth parameter reached a maximum, flashover occurred inside the 
compartment. Flashover is defined as the point at which the entire upper-layer in the 
compartment ignites. After flashover, the fire growth parameter is near zero. The fire 
was considered to be in a steady-state (termed the post-flashover period). The fire enters 
the decay period of its life when it had depleted almost all of the fuel in the compartment. 
Extinguishment occurred when all of the fuel in the compartment had been depleted. 


2.3.2 Temporal Averaging of Data 


Two types of data averaging were performed, with the type of averaging depending 
on the sampling location. When gases were sampled in the exhaust duct, data was 
averaged over a single window of time which encompassed a large portion of the post- 
flashover period. Data was averaged over several windows of time when gases were 
sampled from the compartment and hallway. Details are given below for the procedures 


of both types of averaging. 
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Figure 2.9 The utilization of the fire growth parameter(O) and the compartment upper- 
layer temperature (—) to describe the different stages of growth occurring 
during a compartment fire. 
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2.3.2.1 Averaging over a Single Window of Time 


Post-hallway sampling was performed in the exhaust duct to determine the overall 
oxidation of the gases in the hallway. Since the highest levels of CO were formed in the 
post-flashover period of the fire (typically 140 seconds in length), the averaging was 
performed during this period of the fire. The species levels were quite uniform over the 
majority of the post-flashover period, thus averaging was performed near the middle of 
the period. If external burning occurred in the hallway during the post-flashover period, 
the data was averaged during this time window. The averaging window of time was 
typically 60 seconds in length. The data was averaged during a window of time where 
the fire growth parameter was approximately zero. The FORTRAN program 
FIREAVG3.FOR, see Appendix B, was used to average the data. 


2.3.2.2 Averaging over Several Windows of Time 


Sampling was performed in the compartment and the hallway to investigate the 
evolution of the gases within the facility. To study this transient phenomenon, the data 


was averaged over 4 second time periods before and after flashover. 


Since sampling was performed at one location during an experiment, 25 separate 
experiments were performed to map the species concentrations in the hallway and 
compartment. To view the distribution of the gas concentrations within the hallway and 
compartment at a certain point in the fire, all of the experiments had to be temporally 
aligned with one another. Unfortunately, the time of ignition was not a good reference 
because the growth time of the fire to flashover was variable. The steady state time of 
each experiment, estimated as the time at which the fire growth parameter goes to zero, 


was instead used to align the experiments with one another. 


A FORTRAN program was used to perform the multiple window averaging. Each of 
the output files produced from the program contained data from a single averaging 


window for all of the experiments of interest. 
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2.3.3 Ignition Index 


From a fire safety point of view, it is important to be able to predict the ignition of 
gases in the hallway. A relationship based on classical empirical relations for lean 
flammability limits, termed the ignition index, was derived by Beyler (1984) to predict 


the occurrence of ignition of an upper-layer inside a compartment. The ignition index, 


xX AHS, 
iy eee, (2.16) 
j 7 Nyrodc pal 
where, 


j - fuel species of interest (CO, UHC and H, are considered) 
X; - wet mole fraction of the species j 


AH, ; - heat of combustion of the species j, [kJ/g-mol] (1411 kJ/g-mol for C,H,, 
238 kJ/g-mol for CO, and 242 kJ/g-mol for H) 


Ts,; - adiabatic flame temperature of the stoichiometric mixture for fuel species 
j at the lower flammability limit, [K] (1700 K for UHC, 1450 K for CO 
and 1080 K for H,) 


To - temperature of the stoichiometric mixture of fuel species prior to reaction, 


[K] 


Nproq ~ Number of g-moles of products of complete combustion per g-mole of 
reactants (fuel species mixture) 


C,; ~ average heat capacity of products of complete combustion, [kJ/g-mol K], 


was calculated using a FORTRAN program, and was used in this study to attempt to 
predict the occurrence of external burning in the hallway. The adiabatic flame 
temperature for a stoichiometric mixture at the lower flammability limit is approximately 
100 K higher than the adiabatic flame temperature at the lower flammability limit. Note 
that the mixture at the lower limit is stoichiometric. The emphasis on having a 
stoichiometric mixture is that the ignition of a diffusion flame will locally occur where 


the mixture is at approximately stoichiometric conditions. 
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The ignition index is the ratio of the potential energy of the gases divided by the 
energy necessary for ignition of the gases. Once the potential energy is greater than the 
energy necessary for ignition, /. rises above unity and ignition is expected to occur. 


This is of course assuming that some source of ignition, a pilot, is present. 


The ignition index was calculated for gases sampled from the center of the opening 
connecting the compartment and the hallway. A few assumptions were applied in order 
to determine the JJ. The initial temperature of the reactants, Tg in Eqn. (2.16), was 
assumed to be equivalent to the average temperature of the top three thermocouples in the 
compartment upper-layer. The UHC were assumed to only consist of ethylene (C,H,). 
As previously mentioned, wet concentrations of the species were based on the fact that 
CO, and H,O are present in stoichiometric proportions. Finally, the concentration of H, 
had to be estimated. In hood experiments performed by Beyler (1983) with liquid 
hexane, the concentration ratio of H, to CO varied from 0.26 to 0.67. Gottuk (1992) 
determined that a ratio of H, to CO of 0.5 yielded ignition index results which were 
+10% of the expected value of J. The ratio of H, and CO equal to 0.5 was also utilized 
in this study. 

The thermodynamic data was taken from several sources. The values of C, for the 
products of complete combustion was found in Van Wylen and Sonntag (1986), while the 
values of AH, were obtained form Beyler (1984) and Drysdale (1985). The adiabatic 
flame temperatures of a stoichiometric mixture at the lower flammability limit, Ts, ;, were 


also obtained from Beyler (1984). 
2.3.4 Video Record of Experiments 


Each experiment was recorded using a VHS video camera which was located at the 
end of the hallway where the effluent enters the fume hood. Video taping provided a 
permanent record of each experiment, and allowed the opportunity to observe 


experiments again to document the occurrence of important events during the test. 
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The external burning in the hallway was one of the most important visual 
observations in the experiments. Video records of the experiments were viewed before 
data averaging was conducted to verify the occurrence(s) of external burning in the 
hallway. Accurately determining the times when the external burning begins and ends in 
the hallway was also important when attempting to predict the ignition of gases in the 


hallway, and verifying trends in the measured gas concentrations. 


Video record of the experiments also played an integral role in visualizing the flow 
of gases in the hallway. Since laboratory techniques for measuring the flow field in the 
hallway were not practical for the reported experiments, the flame emissions in the 
hallway (external burning) were used to reveal how the gases were transported away from 
the compartment and down the hallway. The path of the combustion gases was most 
easily visualized when the most significant burning occurred in the hallway. This 
occurred when the gases exiting the compartment encountered concentrations of oxygen 
near those measured in the ambient surroundings (i.e. experiments with an exit soffit of 


0.10 m or less). 


2.4 TYPES OF EXPERIMENTS 


The remainder of this chapter is dedicated to discussing the different types of 
experiments performed to measure the transport of CO away from a compartment on the 
side of a hallway. Initially, the goals of the experimental execution are explicitly defined. 


The experiments conducted to accomplish these goals are then described. 
2.4.1 Experimental Goals 


There were two main goals in the study. The first goal was to quantify the 
conditions where high concentrations of CO will be transported to locations downstream 


of the compartment where gas temperatures have decreased to levels where CO becomes 
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non-reactive, termed remote locations. The controlling factor in the transport of high 
levels of CO to remote locations was hypothesized to be the air entrainment into the 
plume of combustion gases exiting the compartment. Secondly, the study investigated 
the evolution of CO within a hallway when the burning compartment was located on the 
side. The evolution of CO in the hallway was hypothesized to be a function of the time 
necessary for the hallway gas concentrations to reach a steady-state, the occurrence of 


external burning in the hallway and possibly the type of fuel in the compartment. 


2.4.2 General Experimental Procedure 


Prior to the ignition of the compartment fire, several procedures were conducted to 
ensure the accuracy of the acquired data and to limit the number of experimental 


variables. 


To be confident that the data acquired was accurate, the gas analyzers were calibrated 
for the desired range before each experiment. The load cell was calibrated at the 
beginning of the day of experimentation. After every 5 experiments, sampling lines and 
aspirated thermocouple lines were cleared of soot which had collected on the walls of the 


lines by blowing high pressure air through the lines. 


The variations in the initial gas temperature in the compartment and hallway have 
previously been noted to have effects on CO formation and oxidation. For this reason, all 
experiments were conducted with a compartment temperature of 20-30°C, and a hallway 
temperature of 15-25°C. Wall temperatures also have an effect on the heat losses to the 
surroundings. To ensure approximately the same losses, the facility was allowed to cool 
(for 45 to 75 minutes depending on the ambient temperature) until the air temperatures 
within the facility were within the ranges state above. Lastly, the initial weight of the 


fuel inside the compartment was kept as constant as possible. 
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2.4.3 Fuel Samples 


Three types of fuels were used in the investigation, liquid hexane, polyurethane 
foam, and Douglas fir plywood,. The majority of the experiments being conducted using 
liquid hexane. The composition of each fuel in addition to some of the fuel properties are 
shown in Table 2.1. The type of fuel was chosen based on the existence of experimental 


compartment data, and the levels of combustion gases produced by the compartment fire. 


A number of experimental studies have utilized liquid hexane, CgH,,, as a fuel since 
its composition is well characterized. Experiments performed by Beyler (1983) and 
Gottuk (1992a) were used as a basis for the species production in liquid n-hexane 
compartment fires. Liquid n-hexane was also used in the study by Ewens (1994a) where 
the oxidation of combustion gases down a hallway. To expand on this data base of 
information, liquid n-hexane was used in the majority of the experiments in this study. 
The liquid n-hexane was burned in 0.064 m deep steel fuel pans which were 0.15 m to 
0.28 m in diameter. The pans were completely filled with liquid n-hexane just prior to 


the ignition of the fuel. 


To validate that the experimental results were not fuel specific, polyurethane foam 
was also burned inside the compartment. The polyurethane foam had a density of 44.1 
kg/m° and an elemental composition of CH, 7409 323No.79 (“Special Reference”, 1976). 
An elemental analysis performed by Galbraith Labs confirmed these results (Gottuk, 
1992a). The foam contained a filler which accounted for 45% of the weight. After the 


Table 2.1 The composition and some properties of the fuels utilized in the study. 


Volatiles (F/A), Heat of Combustion 
Fuel Composition (by mass) (LHV), AH,, [kJ/kg] 


Liquid n-Hexane 0.0658 44,735 
Polyurethane Foam CH, 7400 323No.608* 0.113 26,570 
Douglas fir Plywood 0.211 26,657 


* Not including nonvolatile filler 
** Assuming 15% char 
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foam compartment fire had extinguished, the filler remained as a powder and solid crust. 
A 0.457 m square, 0.30 m thick slab of polyurethane foam was burned in a stainless steel 
pan in the center of the compartment. Ignition of the foam block was instigated by 


placing two lit wooden matches on top of the block. 


In a number of liquid hexane experiments, Douglas fir plywood was hung 0.05 m 
below the ceiling of the compartment to investigate the effects of an oxygen containing 
fuel in the upper-layer on species formation in the compartment upper-layer. Three 
0.30 m wide, 0.91 m long 0.654 cm thick sections of Douglas Fir plywood were 
suspended below the ceiling so the mass loss rate of the plywood could be measured. 
The elemental composition of the plywood was determined to be CH) 590, 97; (assuming 
15% of the wood mass formed char) from an elemental analysis performed by Galbraith 


Labs. These experimental results are presented in Appendix C. 
2.4.5 Estimation of Compartment Fire Gas Jet Air Entrainment 


The oxygen entrainment into the compartment fire gas jet issuing into the hallway 
was hypothesized as the parameter which primarily governed the degree of CO oxidation. 
Overventilated compartment fires (<1) were performed to quantify the air entrainment 
into the compartment fire gases for different opening sizes (0.12, 0.08 and 0.04 m’) and 
for cases with no inlet soffit and with a 0.20 m inlet soffit. The automated cart was 
utilized in these experiments to generate species concentration profiles of the combustion 
gas layer present in the hallway during the steady-state time of the compartment fire. The 
entrainment was determined by the increase in O, concentration within the gases as they 


flowed away from the compartment. 
2.4.5 Study of Conditions Necessary for the Transport of High CO Levels 


Controlled experiments were conducted to determine the conditions under which 


high levels of CO were transported to remote locations. To quantify the overall effect of 
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each experimental variable on CO, CO, and UHC levels, gas sampling was performed in 


the exhaust duct. 


The air entrainment into the combustion gases exiting the compartment was 
hypothesized as the primary variable in determining the yield of CO transported. The air 


entrainment was experimentally varied by changing: 
1. the height of the inlet soffit, 

2. the opening size, and 

3. the depth of the upper-layer in the hallway. 


The majority of the experiments were performed with a 0.20 m soffit, but some 
experiments were performed with no inlet soffit for comparison purposes. Opening sizes 
of 0.12 m’, 0.08 m? and 0.04 m? were utilized in the study to vary the size of the jet 
entering the hallway. The depth of the low oxygen containing upper-layer in the hallway 
was experimentally controlled by the exit soffit height which was varied from 0.0 to 0.80 


m in 0.10 m increments. 


In addition to the air entrainment parameters, the investigation also included the 
effects of the compartment global equivalence ratio, the absence or presence of external 


burning and the type of fuel on the species levels transported downstream. 


When a door connects the compartment and the hallway (instead of a window), air 
enters the compartment through the lower portion of the door and combustion gases exit 
the compartment through the upper portion. The fluid mechanics at the compartment and 
hallway interface are different than those in the case with a window containing 
unidirectional flow. The effects of having a 0.50 m wide, 0.75 m high door connecting 
the compartment and hallway on the species production inside of the compartment, the 
oxidation of gases in the hallway and the transport of gases to remote locations is 


discussed in Appendix C. 
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2.4.6 Study of Species Concentration Evolution within the Facility 


The species evolution in the hallway was investigated under conditions where high 
post-hallway yields of CO (the yield of CO exiting the hallway) have been measured. 
The effects of the occurrence or absence external burning, the time necessary for hallway 
gas concentrations to reach a steady-state, and the fuel type in the compartment on the 
species levels in the hallway were investigated. The evolution of the species in these 
cases was measured by performing 25 different experiments at a variety of locations 
within the hallway and compartment, see Fig. 2.10. All samples were taken 0.05 m 
babar the ceiling of the hallway and compartment except for the sample location at the 
opening (location b). At location b, the sampling was performed in the center of the 
opening, approximately 0.33 m below the ceiling. The sampling locations are shown on 


each plot of the species variation in the facility. 


The upper-layer height in the hallway was observed to have an effect on the 
occurrence or absence of external burning in the hallway. Liquid n-hexane pool fire 


experiments were conducted with a 0.60 m exit soffit to study this effect. 


The effect of the time necessary for the species concentrations in the hallway to reach 
a steady-state on the species evolution was investigated by varying the conditions at the 
hallway exit. In experiments with the bottom portion of the hallway exit blocked, nearly 
the entire post-flashover period (of the compartment fire) was necessary for the gas 
concentrations to reach a steady-state in the hallway. These experiments were termed as 
having a long transient time. Experiments containing exit soffit heights of 0.60 m 
resulted in the hallway gas concentrations reaching a steady-state approximately 50-60 
seconds after flashover. These experiments were termed as having a short transient time 


experiments. 


The effect of the compartment fire fuel on the evolution of species in the hallway 


was investigated by comparing liquid n-hexane fire experiments with polyurethane foam 
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fire experimental results. A coarse mapping of the hallway was performed in the 
polyurethane foam experiments, see Fig. 2.11. The mapping was designed to capture the 


non-uniform concentration of CO across the hallway which was seen to exist in the n- 


hexane experiments. 
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CHAPTER 3. 


RESULTS: TRANSPORT OF CHEMICAL SPECIES 
TO REMOTE LOCATIONS 


3.1 INTRODUCTION 


The oxygen entrainment into compartment fire gases entering the hallway was 
hypothesized in section 1.3 as being the parameter which controls the level of species, 
particularly CO, transported to remote locations. Overventilated n-hexane compartment 
fire experiments were performed to measure air entrainment into the gases entering the 
hallway. The entrainment was measured for different opening sizes and different soffit 
heights above the opening. Underventilated fire tests were performed to measure the 
effect the experimental variables which alter the air entrainment on the post-hallway yield 
levels. The results of a parametric study quantifying the effects of the opening size, soffit 
height and upper-layer depth on the post-hallway CO, UHC, and CO), yields are 
presented. The effect of varying the compartment global equivalence ratio is also given. 


3.2 AIR ENTRAINMENT INTO COMPARTMENT FIRE EXHAUST GASES 


Oxygen entrainment is hypothesized as being the most significant parameter which 
determines the species levels transported to remote locations. Mullholland et al. (1991), 
Morehart et al. (1991), and Morehart et al. (1992) investigated limited air entrainment 
into compartment fire plumes, but the effects of oxygen entrainment into gases exiting 
the compartment has been addressed in only one study (Ewens, 1994a). In this section, a 
simplified model for determining the air entrainment into a non-reacting jet of 


combustion gases entering the hallway is derived from the conservation of mass. The 
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model is then applied to data produced from overventilated compartment fires to quantify 


the effects of different window and soffit sizes on the air entrainment into the jet. 
3.2.1. Background on Air Entrainment into Ceiling Jets and Plumes 


Air entrainment into a jet of gases is defined here, after Ellison and Turner (1959), in 


terms of the entrainment function, 
V, 
E=+, 3.1 
7 (3.1) 


which is defined as the ratio of the average velocity of entrained gas V, with the average 
velocity of the jet, V; (Ellison and Turner, 1959). Air entrainment into ceiling jets and 
plumes is a function of the distance the gases have traveled from the source, the viscous 
effects between the jet gases and the surrounding gases (the Reynolds number), and the 


gravitational effects (the Richardson number). 


In experiments with turbulent, unheated jets, Bakke (1957) and Poreh (1967) 
determined that the value of the entrainment function to be nearly independent of the 


distance traveled and the Reynolds number, 


Re = PPh (3.2) 
Hi; 


making the entrainment mainly a function of the Richardson number, 


pj = 22» (Pe Pi) (3.3a) 
beds Fe 
where 
4A 
Die —- (3.3b) 


= | 
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These findings were supported by the salt water floor jet experiments performed by 
Ellison and Turner (1959). As a result of their work, Ellison and Turner (1959) 


determined that, 
E = 0.075 exp(-3.9 Ri), (3.4) 
the entrainment function decayed exponentially in magnitude with an increase in 
Richardson number. 
3.2.2 Air Entrainment Model for Compartment Fire Exhaust Gases 


To estimate the air entrainment into a non-reacting jet of compartment fire gases, a 
simplified model was derived using a control volume analysis. Shown in Fig. 3.1 is a 


control volume (CV) drawn around a portion of upper-layer gases in the hallway. 


Hallway Ceiling 


Figure 3.1 The simplified control volume analysis of the hallway upper-layer gases. 


Assuming one-dimensional flow, the gases flowing into the CV, m,, are less diluted than 


those flowing out of the CV, m,, due to the entrainment of air, 7m,. 
Applying the conservation of mass to this CV yields 
m, +m, = m,. (3.5) 


Since the molecular weight of all the gases considered is nearly that of air, 28.84 


kg/kmol, the above relation can be transformed to molar flow rate form; 
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i,+n, =h,.- (3.6) 
Performing a similar analysis for the flow of O, within the system, the following relation 
is produced 


Kio, aN Ne, oa No.0, (3.7) 


If equation (3.7) is divided by 7, and equation (3.6) is applied, then 


ey 
Ni; +n, Ni; it No 


1:0, NO, lel No.0, 3 8) 


Dividing the top and the bottom of the first term on the left hand side by 7, and dividing 
the top and the bottom of the second term on the left hand side by 7,, equation (3.8) is 


transformed to 


I 


Xo hogs 
: rf “fF : z =X,0.> 
Ca ee od CLD « 


(3.9) 
where X denotes the oxygen volume fraction at the location denoted by the subscript. 
Through some algebra, a quadratic equation in terms of 7,, | 
(pa ; Xoen) 12 (Xap + X.0, -2X,0) 
any pea E + PE Ee SEAT LSS SE ST Re I a 


N; Ni; 


hn, +(Xio, -Xp0,)=0 (3.10) 


is formed. With the restriction of 7,>0, the quadratic formula yields 


sous —b+Vb’ —4ac 


3.11 
e ae (3.11) 


(orgs 1X, Oo 20a) 
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The mass flow of the air entrainment into the CV, 
Mm, x MW ile (3.12) 
is determined by multiplying by the molecular weight of air, MW,,,,. 


In summary, to determine the air entrainment mass flow rate into a non-reacting jet 


of gases between two points, the following must be known: 
e molar flow rate of gases into the CV, 
e volume fraction of O, flowing into CV 
e volume fraction of O, entrained into CV, and 
e volume fraction of O, flowing out of CV. 


Since the mass flow rate out of the compartment was known, the molar flow rate into 
the control volume was estimated by dividing the mass flow rate by the molecular weight 
of the gases (assumed to be that of air). By measuring the concentration of the gases 
exiting the compartment and at some location downstream of the compartment, the air 
entrainment which has occurred between these two points can be determined assuming 


that the O, concentration in the gases being entrained is known. 
3.2.3 Air Entrainment into Ceiling Jets 


A ceiling jet is defined as hot gases flowing along the ceiling of the hallway. The 
compartment fire gases entered the hallway as a ceiling jet when no inlet soffit was 
present. These tests are termed ceiling jet experiments. The air entrainment into 
overventilated compartment fire gases entering the hallway as a ceiling jet was estimated 
for three different window sizes, 0.12, 0.08 and 0.04 m’. A 0.10 m diameter fuel pan 
filled with liquid n-hexane was utilized in the experiments to produce global equivalence 


ratios ranging from $=0.1-0.4. Experiments were performed with no exit soffit or 


blockage at the exit of the hallway. 
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Since the species concentrations in the upper-layer were not constant with height, 
species concentration profiles at a particular location in the hallway were experimentally 
determined by moving the sampling probe during the steady-state time of the fire. The 
sampling probe was moved from 0.05 m below the ceiling to 0.30 m below the ceiling in 
0.05 m increments. With the steady-state lasting approximately 300 seconds in these 
experiments, the probe remained at each location for 20 seconds before being moved 
down another 0.05 m. In the underventilated fire experiments, the majority of species 
oxidation was measured to occur as the combustion gases exit the compartment and move 
across the hallway. The measurement of species profiles at the opening and 0.60 m 
across the hallway was, therefore, deemed sufficient for quantifying the difference in the 


entrainment between the three opening sizes. 


To apply Egn. (3.11), the volume fraction of O, at each location must be known. An 
integrated average of the measured O, profile was used as the volume fraction of O, in 
the upper-layer at each sampling location. Air was assumed to be entrained into the 
upper-layer gases, so X,9 was assumed to be 0.21. Note that the entraining gases may 
not always contain 21% QO, (i.e. when deep upper-layers are accumulated in the hallway). 
This is an important point to considered when determining the oxygen entrainment into 
the compartment fire exhaust gas jet. This is addressed in Chapter 5 when correlations to 
predict CO and UHC post-hallway yields are developed based on the oxygen 


entrainment. 


Since the mass flow rate of the gases exiting the compartment was known, an oxygen 
profile was measured at the opening where the compartment fire gases entered the 
hallway. The other O, profile was measured at a location 0.60 m downstream of the 
compartment opening. The mass flow rate of gases entering the hallway was converted 
to molar flow rate by dividing by the molecular weight of the gases (assumed to be 


equivalent to air). 
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Table 3.1 The experimental conditions in the overventilated fire experiments used to 
determine the air entrainment into the ceiling jet. 


Upper- Entrained 
Layer Gas Gas 

Opening, | Density (as Density (as 
m/s air), kg/m* air), kg/m* 


Mass Flow | Velocity of 
Gases at 


Opening 
moa 


Table 3.2 The air entrainment into the ceiling jet using the overventilated compartment 
fire experimental results. 


Oxygen Oxygen 
In, mole Out, mole 
fraction fraction 


The conditions and the results of the analysis are tabulated in Table 3.1 and Table 
3.2, respectively. To compare the air entrainment results with those of Ellison and Turner 
(1959), the entrainment function, see Eqn. (3.1), must be determined for each case. The 
average velocity entering the hallway, V; , was determined by dividing the mass flow of 
gases through the opening by the density (taken to be the density of air at the gas 


temperature) and the opening area. An average entrainment velocity, 


ae (3.13) 
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was determined by dividing the mass flow of entrainment by the density of the entrained 
air (density of air at the gas temperature below the upper-layer) and the wetted area of the 
jet. The wetted area of the jet, 

Aye =U(2h+w), (3.14) 
was defined as the product of the length, /, between the two sampling points (0.60 m) and 
twice the height, h, plus the width, w, of the opening. For the ceiling jet experiments, the 


Richardson number, 


4p; gA2 = 
Ri= 40: 84open (Pe= Ps). (3.15a) 
m(2h+w)  p, 
where; 
meets 
D, =——"— (3.15b) 
(2h+w) 


was determined using the conditions at the opening. 

The air entrainment function is observed in Fig. 3.2 to exponentially decrease with 
respect to the Richardson Number , and is represented by the curve: 

E = 0.07 exp(-1.1Ri)+0.006 . (3.16) 

The variation of the entrainment function with the Richardson number is similar to that 
found by Ellison and Turner (1959), see Eqn. (3.4), in their salt-water experiments, but is 
slightly shifted upward. The upward shift in the value of E£ is attributed to both not 
considering the expansion of the ceiling-jet as it entrains air and the wall effects. Not 
accounting for the jet expansion results in a smaller surface area, A,,.,. over which the 
gases are entrained into the jet. This results in a calculation of a higher average 
entrainment velocity, see Eqn. (3.13). The wall effects were noted by Ellison and Turner 
(1959) to induce further entrainment into the jet by forcing mixing near the wall. The 
model does, however, provide a reasonable estimation of the air entrainment into the 


compartment fire gases. 
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Figure 3.2 The entrainment into a ceiling jet produced from an overventilated 
compartment fire. The data is from experiments having opening sizes of O- 
0.12 m’, 01-0.08 m’, and A-0.04 m? with the solid line being the curve fit to 
the data. The dashed line is the ceiling jet curve from Ellison and Turner 


(1959). 
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The effect of the opening size on the mass flow of air entrainment is also seen in 
Fig. 3.2. As the opening size was decreased, the air entrainment was measured to 
increase. 


3.2.4 Air Entrainment into Plumes 


A plume is defined as a jet of hot gases which are permitted to rise vertically due to 
buoyancy. In experiments with a 0.20 m soffit, the gases exiting the compartment were 
allowed to rise 0.20 m until they impinged upon the ceiling approximately 0.30 m 
downstream of the opening. After ceiling impingement, the gases traveled across the 
hallway as a ceiling jet. Air entrainment into plumes has been noted by Ellison and 
Turner (1959), Alpert (1972, 1980) and Hinkely et al. (1984) to be considerably larger 
than that observed in ceiling jets. The turning region, where the jet impinges on the 
ceiling, was noted as having the next highest entrainment rate, while the ceiling jet had 
the least efficient air entrainment. The air entrainment into the fire gases as they flowed 


across the hallway first as a plume and then a ceiling jet is quantified in this section. 


Overventilated compartment fire experiments were performed with a 0.15 m 
diameter fuel pan and no orifice on the air inlet duct. The three different opening sizes 
(0.12, 0.08 and 0.04) used in the ceiling jet study were also used in these experiments. 
To capture the three regions of the flow (the plume, the turning region and the ceiling jet) 
present as the gases moved across the hallway, species concentration profiles were 
measured at the opening (0.0 m), 0.15 m, 0.30 m, 0.45 m, 0.60 m and 0.90 m across the 
hallway. The profiles were produced by varying the sampling location from 0.05 m to at 
least 0.35 m below the ceiling in 0.05 m increments (remaining at each location for 20 
seconds) during the steady-state time of the experiment. The molar flow rate into the 
initial control volume, 7,, was determined from the mass flow of gases out of the 
compartment with each subsequent control volume using the value of 7, for its value of 
nh, 


i? 
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Figure 3.3 
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The mass flow rate of air entrainment into the overventilated compartment 
fire gases in experiments with a 0.20 m inlet soffit where the plume, the 
turning region and the ceiling jet were present in the initial 0.60 m of travel. 
Experiments utilized O-0.12 m’, 01-0.08 m’, and A-0.04 m? opening sizes. 
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Table 3.3 The experimental conditions in the overventilated fire experiments with a 0.20 
m inlet soffit where the air entrainment into compartment fire gases in the 
initial 0.60 m of the hallway is determined. 


Entrained 


Mass Flow | Velocity of Upper- 


Opening Opening Gases at Layer Gas Gas 
Area, m? Width, m Opening, | Density (as Density (as 
m/s air), kg/m° air), kg/m* 


[008 | 0.16 | 050 | 0.0613 | 1.08 | 0.7105 | 0.8905 _ 


Table 3.4. Using overventilated compartment fire experimental results with a 0.20 m 
inlet soffit, the calculated air entrainment into the compartment fire gases in 
the initial 0.60 m where the plume, the area of impingement and the ceiling jet 
are all present. 


Oxygen 
In, mole 
fraction 


[008 [0.175 | 0.193 | 0.060 | 0.10 | o41__| 0.097 _] 


Oxygen 
Out, mole 
fraction 


For all three opening sizes, the air entrainment is seen in Fig. 3.3 to be the most 
significant in the first 0.30 m of the hallway where the plume was present. As noted in 
previous experiments performed by Hinkely, et al. (1984), less significant was measured 
to occur in the impingement area (0.30<y<0.45 m) and in the area of the ceiling jet 
(0.45<y<0.90 m), see Fig. 3.3. Also, the entrainment into the compartment fire gases was 


measured to increase as the area of the opening was decreased. 


The entrainment function for the first 0.60 m of the hallway was calculated, see 
Table 3.3 and 3.4, so quantities could be compared with the ceiling jet data. Assuming 


the plume impinges on the ceiling 0.3 m across the hallway, the wetted area of the jet is 
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Ane, = 0.3(4h+3w). | (3.17) 


The Richardson number of the gases entering the hallway was defined as 


49° gA? —p, 
pj = 1: B4oven (Pe p;) (3.18a) 
n¢(2h+2w) pp, 
where; 
4A 
Das =——— - (83.18b 
" (2h+2w) Se 


Utilizing the new definitions for the wetted area and the Richardson number, the 
entrainment function was calculated to be significantly higher than that measured in the 


ceiling jet experiments, see Fig. 3.4 and Table 3.4. The entrainment function 
E =0.14exp(—4.0Ri)+0.084 . (3.19) 


still, however, decreases exponentially with an increase in the Richardson number. As 
before, the entrainment functions for the 0.12 m” and the 0.08 m” openings are shown in 
Fig. 3.4 to be nearly equivalent, but the magnitude of E has risen to 0.10. The presence 
of a 0.04 m? opening resulted in a 50% increase in the entrainment function to a value of 


Os153 
3.2.5 Additional Parameters Affecting Air Entrainment 


Average velocities were used in the air entrainment model since the measurement of 
the flow field was not possible. It should be mentioned that combustion instabilities and 
higher gas temperatures play a role in the actual air entrainment into compartment fire 
gases in experiments where the compartment fire is underventilated. The air entrainment 
model only considers a one-dimensional flow field, thus decreasing the accuracy of the 


model in areas where two and_ three-dimensional flow field’ effects 
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Figure 3.4 The air entrainment function for the overventilated compartment fire gases in 
the first 0.60 m of the hallway in experiments with a 0.20 m inlet soffit where 
both a plume and a ceiling jet are present. The data is from experiments 
having opening sizes of O-0.12 m?, 1-0.08 m’, and A-0.04 m? with the solid 
being the curve fit to this data. The dash-dot line is the ceiling jet curve from 
this study while the dotted line is the ceiling jet curve from Ellison and 


Turner (1959). 
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are significant. For this reason, the measurement of the entrainment was confined to 
regions close to the compartment where wall effects on the flow field were not as 
significant. These experiments, however, do provide a basis for comparing the effects of 


different opening sizes and the presence of a 0.20 m inlet soffit on the air entrainment. 


3.3 EFFECT OF ENTRAINMENT ON POST-HALLWAY SPECIES YIELDS 


The oxygen entrainment into the compartment fire gases entering the hallway was 


experimentally altered in a number of ways: 
e changing the opening size, 
e changing the size of the fire, 
e changing the air inlet orifice diameter, 
e changing the inlet soffit height, and 
e changing the exit soffit height (or exit blockage). 


The first three experimental variables changed the entrainment into the compartment fire 
gases by altering the velocity of the gases entering the hallway, V;. The inlet soffit height 
at the opening determined whether the gases entered the hallway as a ceiling jet or as a 
plume. The inlet soffit, therefore, governed the entrainment into the compartment fire 
gases at the entrance of the hallway where the highest temperatures were present and 
chemical kinetics rates were fastest. The height of the exit soffit determined the O, 
concentration of the gases being entrained into the compartment fire products by defining 
the depth of the oxygen deficient combustion gas layer in the hallway. The effects of 
each of these variables on the post-hallway species yields (the yields of the species in the 


gases exiting the hallway) is discussed in the following sections. 
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3.3.1 Effect of Opening Size 


Opening sizes of 0.12, 0.08 and 0.04 m’ were utilized in the study to alter the air 
entrainment into the compartment fire gases entering the hallway. As seen in Fig. 3.2 and 
Fig. 3.4, a decrease in the opening size results in an increase in the entrainment. The 
oxidation of incomplete compartment fire products entering the hallway is dependent on 
the amount of oxygen available. A decrease in the opening size is, therefore, expected to 
result in lower post-hallway yields of CO and UHC (with the four remaining effects 
mentioned in the beginning of the section held relatively constant). It should be noted 
that the reduction in the yield is related to the oxygen entrainment which occurred in the 
region of the hallway where temperatures were high enough (>850 K for CO and 800 K 
for UHC) for oxidation to occur. After this point, gases were merely diluted by the 


process of entrainment since the chemical kinetics were frozen. 


As stated in section 1.2.2., Ewens (1994a) investigated the oxidation of compartment 
fire gases down a hallway with near ambient levels of oxygen being entrained into the 
compartment fire gas jet entering the hallway. Reevaluating the data taken by Ewens 
(1994b), the effects of the opening size on the post-hallway species yields are shown in 
Fig. 3.5. As anticipated from the air entrainment results seen in Fig. 3.4, the CO and © 
UHC yields are seen in Fig. 3.5 to decrease in magnitude with a decrease in the opening 
size. The increase in the oxidation of CO and UHC is also observed in Fig. 3.5 by an 
increase in the CO, yield. It should be noted that these trends become invalid when the 
compartment fire gas jet entering the hallway is entraining gases containing oxygen levels 
significantly lower than those observed in the ambient surroundings. This topic is 


addressed specifically in section 3.3.3. 
3.3.2 Effects of Air Inlet Orifice Diameter and Fire Size 
The size of the fire and the diameter of the orifice attached to the air inlet duct also 


control the mass flow rate of the compartment fire gases entering the hallway. As 
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Figure 3.5 The effect of the opening size on the post-hallway @-CO, HI-UHC and A- 
CO, yields for the cases with Q=448+154 kW, 6=2.57+0.06, a 0.20 m inlet 
soffit, and no exit soffit. Experimental data of Ewens (1994b). 


rp: 


Chapter 3. Results: Transport of Chemical Species to Remote Locations 


mentioned in section 2.2.1, the size of the fire is controlled by the diameter of the fuel 
pan. The larger fuel pan diameters resulted in a larger fire size. The air mass flow rate 
into the compartment was controlled by the diameter of the orifice placed on the end of 
the air inlet duct. A decrease in the orifice diameter resulted in a decrease in the air mass 
flow rate. These two experimental variables, pan diameter and orifice diameter, were 
used in the study to control the compartment GER and the mass flow rate of the gases 
entering the hallway through a particular opening size. 


With the diameter of the fuel pan controlling the fuel mass loss rate and the orifice 
Remcien controlling the air mass flow rate, the compartment GER was affected by 
altering these two parameters. A decrease in the orifice diameter reduced the air mass 
flow rate into the compartment causing an increase in the compartment GER for a given 
fire size. The decrease in the orifice diameter resulted in a deeper upper-layer within the 
compartment. The compartment upper-layer temperature fell by 125 K when the orifice 
diameter was decreased from 0.30 m (no orifice) to 0.10 m, see Fig. 3.6. The temperature 
drop was attributed to less complete combustion of the volatilized fuel. With the fuel 
volatilization rate dependent on the radiation incident upon the exposed fuel surface, the 
drop in upper-layer temperature could drastically affect the mass loss rate of fuel. The 
fuel mass loss rate, however, is shown in Fig. 3.7 to remain relatively unchanged with a 
decrease in orifice diameter. The decrease in the incident radiation from the fall in the 
upper-layer temperature was obviously compensated by a decrease in the distance 
between the upper-layer and the fuel surface (an increase in the upper-layer depth). For a 
particular fire size, the compartment GER could, therefore, be increased by decreasing the 


diameter of the orifice connected to the air inlet duct. 


The effects of fire size on the post-hallway species yields with no orifice attached to 
the air inlet duct (0.30 m diameter) were determined through further analysis of the data 
of Ewens (1994b). An increase in the fire size resulted in an increase in the mass flow of 


air drawn into the compartment. This resulted in an increase in the mass flow rate of 
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Figure 3.6 The variation in the compartment upper-layer temperature profile with inlet 
orifice diameter of O-0.30m (no orifice), (J-0.25 m, A-0.20, and V-0.15 m 


74 


Fuel Volatilization Rate, m_, [kg/s] 


Figure 3.7 


Chapter 3. Results: Transport of Chemical Species to Remote Locations 


0.016 


0.014 


0.012 


0.010 


0.008 


0.006 


0.004 


0.002 


0.000 
0 100 200 300 400 500 600 700 800 900 


Adjusted Time Using the Fire Growth Parameter, t, [sec] 


The variation in the fuel volatilization rate in n-hexane fires with inlet orifice 


diameter of ----0.30m (no orifice), * * *0.25 m, —*—0.20, and —* e«—0.15 m. 


ip 


Chapter 3. Results: Transport of Chemical Species to Remote Locations 


gases exiting the compartment into the hallway. For a fixed opening size, a large increase 
in fire size is, however, seen in Fig. 3.8 to only slightly increase the momentum of the jet 
entering the hallway. An increase in the jet momentum (lower Richardson number) of 
compartment fire gases entering the hallway was seen in Fig. 3.2 and Fig. 3.4 to result in 
a significant increase in the entrainment of surrounding gases into the jet. The effect of 
the opening size on the momentum of the jet entering the hallway is seen in Fig. 3.8 to be 


much more significant than the effect of the fire size. 


The effect of fire size on the post-hallway species yields was determined from further 
evaluation of Ewens (1994b) data. A broad increase in fire size (for a fixed opening size 
and compartment GER) is shown in Fig. 3.9 to result in only a slight increase in the CO 
and UHC yields with a corresponding decrease in CO, yield. 


3.3.3 Effect of Inlet Soffit Height 


The height of the soffit above the opening connecting the compartment and hallway 
determines whether the gases enter the hallway as a ceiling jet (0.0 m soffit) or as a plume 
(0.20 m soffit). The entrainment into a plume was seen in Fig. 3.4 to be significantly 
higher compared to the entrainment into a ceiling jet (see Fig. 3.2). The effects of the 


ceiling jet and the plume on the post-hallway yields is described in this section. 


The effect of a 0.0 m inlet soffit versus a 0.20 m soffit was determined by 
reevaluating the data of Ewens (1994b). The experiments of Ewens (1994a) were with a 
well-ventilated hallway allowing the effect of the inlet conditions on species yields to be 
studied. The post-hallway CO yields in experiments with a 0.20 m are shown in Fig. 3.10 
to be significantly lower compared with experiments having a 0.0 m soffit. Lower post- 
hallway UHC yields were also measured for experiments with a 0.20 m inlet soffit 
compared to levels measured in 0.0 m inlet soffit experiments, see Fig. 3.11. Conversely, 
an increase in the CO, post-hallway yields is observed in Fig. 3.12 for experiments with a 


0.20 m inlet soffit compared with experiments having 0.0 m inlet soffit. 
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Figure 3.8 The effect of the fire size on the momentum of the compartment fire gas jet 
entering the hallway for opening sizes of O-0.12 m’,O -0.08 m’ and A-0.04 
m” is demonstrated. Experimental data of Ewens (1994b). 
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Figure 3.9 The effect of the fire size on the post-hallway @-CO, HI-UHC and A-CO, 
yields with a 0.12 m opening, ¢=2.35+0.41, no inlet soffit and 0.20 m exit 
soffit is shown. Experimental data of Ewens (1994b). 
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Figure 3.10 The variation in the post-hallway CO yield with no inlet soffit (@) and a 
0.20m inlet soffit (Ml) is shown. The experimental conditions were 
Q=431+150 kW, $=2.64+0.17, and no exit soffit. Experimental data of 


Ewens (1994b). 
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Figure 3.11 The variation in the post-hallway UHC yield with no inlet soffit (@) and a 
0.20m inlet soffit (@) is shown. The experimental conditions were 
Q=431+150 kW, $=2.64+0.17, and no exit soffit. Experimental data of 


Ewens (1994b). 
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Figure 3.12 The variation in the post-hallway CO, yield with no inlet soffit (@) and a 
0.20m inlet soffit (Ml) is shown. The experimental conditions were 
Q=431+150 kW, 6=2.64+0.17, and no exit soffit. Experimental data of 
Ewens (1994b). 
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The significant difference in the post-hallway yields was attributed to the difference 
in entrainment between experiments with a 0.20 m soffit and with a 0.0 m soffit. The 
higher entrainment rate in the 0.20 m soffit experiments (compared with the 0.0 m soffit 


experiments) resulted in lower CO and UHC yields and higher CO, yields. 


3.3.4 Effect of Exit Soffit Height 


The exit soffit height was utilized in the experiments to control the depth of 
combustion gas layer which accumulated in the hallway during the pre-flashover period 
of the compartment fire. In actual building fires, it is hypothesized that the accumulation 
of such a layer in the space adjacent to the burning room results in lower oxygen 
concentration in the gases iashin entrained into the compartment fire gases entering that 
space. The effects of deepening the hallway upper-layer on the downstream species 
levels is explored in this section. All experiments were performed with a 0.20 m inlet 


soffit unless otherwise noted. 


The exit soffit height was varied in the experiments from 0.0 m to 0.80 m in 0.10 m 
increments. The upper-layer depth was visually observed to be approximately equivalent 
to the exit soffit height except for experiments with an exit soffit equal to or less than 
0.20 m. In the cases with exit soffits of 0.0, 0.10 and 0.20 m, the upper-layer depth was 
determined from the hallway vertical temperature profiles, and was found to be 


equivalent to approximately 0.20 m in all three instances. 


The post-hallway yields of CO, UHC and CO, are plotted versus a nondimensional 
upper-layer depth, 
y=6/Z, (3.20) 


which is defined as the upper-layer depth, 5, divided by the distance between the ceiling 
and the bottom of the opening, z, connecting the compartment and the hallway. 
Physically, y represents the degree to which the plume of compartment fire gases entering 
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the hallway was surrounded by oxygen deficient combustion gases. Opening sizes of 
0.12, 0.08 and 0.04 m? corresponded to values of z of 0.44, 0.36 and 0.36 m, respectively. 
For each opening size, the effect of varying the upper-layer depth on the post-hallway 
CO, UHC and CO, yields was measured. The data points were generated by averaging 
60 seconds of temporal data during the compartment post-flashover period. 


With an opening of 0.12 m’ and an upper-layer depth no greater than the inlet soffit 
height, y=0.45, the post-hallway species yields are shown in Fig. 3.13 to be relatively 
constant. The average yield of 0.098 for CO and 0.046 for UHC were the lowest 
measured yields, while the CO, was at its highest level of 2.58. 


With a hallway upper-layer depth below the top of the 0.12 m” opening, the sides of 
the jet started to become surrounded by oxygen deficient combustion gases. This resulted 
in the entrainment of oxygen deficient combustion gases into the jet. The limited oxygen 
entrainment is evident in Fig. 3.13 with the gradual increase in the CO and UHC yields 
and decrease in CO, yield. From the work of Westbrook and Dryer (1984), the reduction 
in available O, reduces the pool of radicals (H, HO,, OH, O, HCO, etc.) available for the 
oxidation of CO and UHC. The rapid generation of radicals is critical to the oxidation of 
CO. Carbon monoxide will not begin to oxidize until the UHC have been reduced to 
levels where the CO oxidation reactions are competitive with the UHC oxidation 
reactions (Westbrook and Dryer, 1984). If the generation of radicals is inhibited, less CO 
will be oxidized before the gas temperature falls below 850 K where the CO oxidation 
rate decreases drastically (Gottuk, 1992a, Pitts, 1992, Yetter et al., 1991). The gradual 


increase in CO and UHC yields is seen in Fig. 3.13 to continue until y>1.4. 


When the hallway upper-layer depth was increased such that y>1.4, the post-hallway 
species yields were measured to be less dependent on the upper-layer depth and more 
dependent on the presence or absence of external burning in the hallway. With external 


burning occurring in the hallway and y>1.4, the post-hallway CO yields were 
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Figure 3.13 The effect the non-dimensional upper-layer height on the post-hallway @- 
CO, M-UHC, and A-CO, yields for experiments with a 0.12 m* opening, 
Q=291+36 kW, and 6=3.20+0.54. Filled symbols represent experiments 
with external burning while open symbols were experiments with no 
external burning. 
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measured to be 0.27 (on average), see Fig. 3.13. These post-hallway yields are 
approximately 23% higher than yields measured inside the compartment by Gottuk 
(1992a) (on average 0.22 for o>1.8). The post-hallway UHC yields were measured to be 
0.24 which is approximately 27% lower than the yields measured inside the compartment 
(on average 0.33 for o>1.8) (Gottuk, 1992b). Again, the lack of oxidation of CO and 
UHC is attributed to the limited radical pool generated in the oxygen deficient 
environment surrounding the jet of hot compartment fire gases entering the hallway 
(Westbrook and Dryer, 1984). From the data, it appears the limited oxidation of CO is 
less than the production of CO from the UHC oxidation. The increase in the post- 
hallway CO yields was, therefore, attributed to the oxidation of UHC which was visually 


evident by the external burning. 


If the upper-layer depth was increased to levels where y>1.7, the jet of compartment 
fire gases issued into a hallway upper-layer environment with low oxygen levels which 
inhibited the ignition of the combustion gases. The post-hallway yields of CO, UHC, and 
CO, were all measured to be equivalent in magnitude to levels measured inside the 


compartment, see Fig. 3.13. 


With a 0.08 m” opening connecting the compartment to the hallway, the distance 
between the ceiling and the bottom of the opening was decreased to 0.36 m. The jet 
issuing from the 0.08 m° opening was also measured to have a slightly higher entrainment 
compared with the jet from a 0.12 m° opening, see Fig. 3.4. The post-hallway yields are 
seen in Fig. 3.14 to show a trend similar to that observed with a 0.12 m* opening (see Fig. 
3.13), but magnitudes are shifted downward when y<1.7. The downward shift is 
attributed to the slightly better O, entrainment in the 0.08 m* opening case. With y 
slightly less than 1.7 and external burning occurring, post-hallway CO yields are 
observed in Fig. 3.14 to be 5% greater than compartment levels. The post-hallway UHC 
yields are, however, seen in Fig. 3.14 to decrease to levels 65% less than in-compartment 


yields. Again, the slight increase in CO and UHC oxidation compared with the 0.12 m* 
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Figure 3.14 The effect the non-dimensional upper-layer height on the post-hallway @- 
CO, M-UHC, and A-CO, yields for experiments with a 0.08 m* opening, 
Q=322+64 kW, and o=2.56+0.41. Filled symbols represent experiments 
with external burning while open symbols were experiments with no 
external burning. 
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opening case was attributed to the increase in oxygen entrainment. The post-hallway CO 


yields were measured to be greater than compartment levels only at a y=1.68. 


A 0.04 m’ opening was utilized to investigate the effect of significantly increasing 
the O, entrainment on the post-hallway yields. The distance between the ceiling and the 
bottom of the opening was 0.36 m (the same as the 0.08 m* opening case). The CO and 
UHC post-hallway yields are seen in Fig. 3.15 to remain nearly equivalent at a magnitude 
of approximately 0.025 for y<1.2. For a 1.2<y<1.7, the CO and UHC yields increase to 
levels observed inside the compartment. Unlike the other two cases, the post-hallway CO 
yields were never measured to be greater than those measured inside the compartment. 


As before, the absence of external burning occurred at a value of y equal to 1.7. 
3.3.5 Effect of All Parameters with a Deep Hallway Upper-Layer 


From the results presented in section 3.3.4, it appears that the accumulation of a 
deep, y21.0, oxygen deficient upper-layer in the hallway during pre-flashover limits the 
oxygen entrained into the compartment fire gases entering the hallway by reducing the 
oxygen concentration of the gases surrounding the jet. The post-hallway species yields, 
therefore, become less sensitive to the entrainment induced by the size of the opening, the 
fire size, and the inlet soffit height. To demonstrate this point, a series of experiments 
were conducted with an upper-layer depth of 6~0.50. The effects of varying the opening 
size, the fire size, the inlet soffit height and the compartment global equivalence ratio on 
the post-hallway species yields were measured. The value of y for these experiments 
varied depending upon the opening size and the inlet soffit height (2.1 for 0.12 m” and 3.1 
for 0.08 m* and 0.04 m” with no inlet soffit, and 1.1 for 0.12 m’ and 1.4 for 0.08 m* and 
0.04 m? with a 0.20 m inlet soffit). 


The post-hallway CO yield data shown in Fig. 3.16 was generated from experiments 
where the opening size, the fire size, the inlet soffit height and the compartment global 


equivalence ratio were all varied. The shape of the curve is similar to that of the 
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Figure 3.15 The effect the non-dimensional upper-layer height on the post-hallway @- 
CO, M-UHC, and A-CO, yields for experiments with a 0.04 m” opening, 
Q=259+41 kW, and o=2.77+0.57. Filled symbols represent experiments 
with external burning while open symbols were experiments with no 
external burning. 
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compartment CO yield versus global equivalence ratio curve determined by 
Gottuk (1992a). In experiments with no external burning, downstream CO yields were 
similar in shape and magnitude compared with yields measured inside the compartment. 
With external burning in the hallway, the downstream CO yields were equivalent to the 
no external burning yields up to a 6 of 2.0. With o>2.0, the yield in experiments with 
external burning was approximately 0.16 compared with 0.22 in experiments without 
external burning. The post-hallway CO yield appeared to be relatively insensitive to the 
other variables (i.e. inlet soffit height, fire size, and opening size) with a deep upper-layer 


present in the hallway. 


The plot of CO yield versus global equivalence ratio can be divided into three 
regions: low, transitional and high toxicity. In the low toxicity region, 0<o<0.6, an 
undetectable amount of CO is transported downstream of the compartment. This is 
expected since there is virtually no CO being produced inside the compartment in this o 
range. Characteristic of the transitional toxicity region, 0.6<<2.0, is a steep rise in the 
CO yield with an increasing . The slope of the rise depends on whether or not external 
burning occurs within the hallway. The slope is 0.16 with no external burning and 0.11 
with external burning. The region is broader and the rise in the CO yield is less steep 
compared to the in-compartment data taken by Gottuk (1992a) where the transitional 
toxicity region was 0.6<o<1.5 and the slope was 0.24 (see Fig. 1.2). The high toxicity 
region in Fig. 3.16 encompasses all @>2.0. The CO yields plateau in this region and 
remain fairly constant. The level at which the CO plateaus depends on whether or not 
external burning occurred inside the hallway. The CO yield averages 0.22 with no 
external burning and 0.16 with external burning. The plateau level with no external 
burning is consistent with the in-compartment results reported by Gottuk (1992a). 
Overall, the CO yield in experiments with external burning decreases by 0.05-0.06 from 
the no external burning case for all @ greater than 1.2, see Fig. 3.16. 
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Figure 3.16 The effect of varying all parameters on the post-hallway CO yield when a 
deep upper-layer, y>1.0, was present in the hallway. Filled symbols 
represent experiments with external burning while open symbols were 
experiments with no external burning. 
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The post-hallway UHC yields are seen in Fig. 3.17 to also have a strong dependence 
on the occurrence of external burning in the hallway. The UHC yields are slightly 
scattered for the experiments with no external burning. The post-hallway UHC yields in 
experiments with external burning are seen in Fig. 3.17 to rise vary with an increase in 9. 
It is interesting to notice that the post-hallway UHC yields were not at a measurable level 
until @ was greater than 1.5. External burning in the hallway prevented measurable levels 


of UHC downstream until a o>1.5. 
3.3.6 Major Parameters Controlling Post-Hallway Species Yields 


The parameters having the most significant effect on the post-hallway species yields 


were the: 
e opening size, 
e inlet soffit height, 
e depth of upper-layer in hallway (hallway exit height), 
e occurrence of external burning, and 
e global equivalence ratio of the compartment fire. 


The first two of the parameters are the variables which affected the entrainment into the 
compartment fire gas jet entering the hallway, while the third parameter controlled the 
oxygen concentration entrained into the jet. The last two parameters govern the 


chemistry which occurred in the hallway and compartment. 
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Figure 3.17 The effect of varying all parameters on the post-hallway UHC yield when a 
deep upper-layer, y>1.0, was present in the hallway. Filled symbols 
represent experiments with external burning while open symbols were 
experiments with no external burning. 
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3.4 POST-HALLWAY SPECIES YIELDS FROM 
POLYURETHANE FOAM FIRES 


A series of experiments were conducted using polyurethane foam as the fuel inside 
of the compartment. The experiments exhibited an upper-layer temperature of 960 K, on 
average, at d=2.1 which is approximately 25 K lower than that observed in n-hexane 
experiments (985 K on average). The post-hallway yields from polyurethane fires are 
presented in this section, and compared to the post-hallway yields of the n-hexane fire 


experiments. 


The experiments were conducted using a 0.20 m inlet soffit, a 0.04 m” opening and a 
compartment global equivalence ratio of 2.1. No sustained external burning occurred in 
the hallway regardless of the exit soffit height. The result of the absence of external 
burning in the hallway on the post-hallway CO yields was a yield greater than 0.15 at all 
hallway upper-layer depths tested, see Fig. 3.18. The post-hallway CO yields were 
measured to be approximately 0.15 at a y=0.55, and increased to a level of 0.19 with 
y>1.0. The post-hallway CO yields for y>1.0 were equal to the yields measured inside the 
compartment by Gottuk (1992a). The low post-hallway CO, yields shown in Fig. 3.18 
further demonstrates the lack of oxidation of CO in the hallway. The UHC yields plotted 
in Fig. 3.18 are approximately 4 times lower than yields measured in n-hexane fires. The 
low UHC yields were attributed to polyurethane fires producing high levels of 
oxygenated hydrocarbons (C,H,O), unsaturated hydrocarbons (C,H,) and higher order 
hydrocarbons. Oxygenated hydrocarbons (i.e. aldehydes and ketones) and unsaturated 
hydrocarbons were reported by Colket et al. (1974) to produce FID responses 
significantly lower than expected from the calibration. For example, C,H,O and C,H, 
produces a signal 60% and 65%, respectively, lower than that of C,H, (Colket, et al., 
1974). The other source of possible error in the UHC measurements was the higher order 


hydrocarbons condensing in the sampling lines. 
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Figure 3.18 The effect of the nondimensional upper-layer depth on the post-hallway O- 
CO, O-UHC and A-CO, yield in experiments with polyurethane foam as the 
fuel inside the compartment. The experiments utilized a 0.04 m? opening, a 
0.20 m inlet soffit, Q=164+11 kW and o=2.08+0.01. Filled symbols 
represent experiments with external burning while open symbols were 
experiments with no external burning. 
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Compared with the n-hexane experiments, the polyurethane foam experiments 
contained higher post-hallway CO yields at lower values of y due to the absence of 
external burning. The post-hallway CO yields in the foam experiments were measured to 
be approximately 14% less than the yields in the n-hexane experiments with no external 
burning. The low post-hallway CO, yields of 1.3 were consistent with the post-hallway 


yields measured in the n-hexane experiments with no external burning. 
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CHAPTER 4. 


RESULTS: EVOLUTION OF SPECIES 
CONCENTRATIONS IN THE HALLWAY 


4.1 INTRODUCTION 


The effect of experimental conditions on temporal and spatial distributions of species 
concentration and temperature in the hallway cannot be measured post-hallway. In- 
hallway sampling was performed to measure the distribution of CO, UHC, CO, O, and 
temperature in the hallway. The data is used to develop a more detailed understanding of 
the parameters controlling the oxidation of CO. The experimental conditions for the in- 
hallway sampling experiments were consistent with those where high CO yields were 
measured post-hallway. The time necessary for hallway species concentrations to reach a 
steady-state (termed transient time), the absence or presence of external burning, and the 
type of fuel burning inside the compartment were experimentally varied. The effects of 
these parameters on the species concentrations and temperatures in the hallway were 


measured and are presented in this chapter. 


4.2 MOVEMENT OF COMBUSTION GASES IN THE HALLWAY 


With the compartment-hallway orientation shown in Fig. 2.1, the flow field in the 
hallway is highly three-dimensional. Since laboratory techniques for measuring the flow 
field in the hallway were not practical for the reported experiments, flame luminosity in 
the hallway (external burning) was used as a flow visualization tool to reveal the flow 
field present in the hallway. The path of the combustion gases was most easily visualized 


with the most intense burning in the hallway. This occurred when the gases exiting the 


96 


Chapter 4. Results: Evolution of Species Concentrations in the Hallway 


compartment encountered concentrations of oxygen near those observed in the ambient 
surroundings (i.e. experiments with an exit soffit of 0.10 m or less). The visualization 
was performed using a video camera placed at the hallway exit, and is described from this 
point of view (as if one were standing at the exit looking down the hallway toward the 


compartment). 


An image of the typical external during the steady-state period (termed the post- 
flashover period) of the compartment fire is shown in Fig. 4.1. As the combustion gases 
exited the compartment, a buoyant plume was formed at the hallway entrance due to the 
presence of a 0.20 m soffit. The plume impinged upon the ceiling resulting in three 
major paths for the combustion gases transport down the hallway. Since there was an 
inlet soffit, the first means of transport was a counterclockwise rotating vortex formed at 
the top left corner of the hallway. From the size of the vortex, it was inferred to be a 
minor path of transport for this geometry. This mechanism of transport disappeared in 
the absence of an inlet soffit. The impingement of the buoyant plume was observed to 
force some of the gases to take on a trajectory which resulted in their transport down the 
center of the hallway. This constitutes the second path of transport. The majority of the 
external burning was observed to be present both across the hallway and along the upper 
right corner of the hallway, see Fig. 4.1. Based on this observation, the bulk of the 
combustion gases were driven across the hallway by the momentum of the jet exiting the 
compartment. The bulk flow of combustion gases then impinged upon the opposite wall 
causing a large clockwise rotating vortex to develop in the upper right corner of the 
hallway, see Fig. 4.1. The size of the vortex was estimated to be approximately 0.20-0.30 
m in diameter. The size of the vortex is determined by the velocity of the gases 
impinging on the wall opposite the compartment. The velocity of the gases impinging on 
the wall is expected to primarily depend on the distance the gases need to travel before 


reaching the wall. The vortex is, therefore, believed to be scaled to the hallway width. 
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Figure 4.1 An image of the external burning in the hallway when looking down the open 
end of the hallway underneath the fume hood. 
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4.3 EVOLUTION IN HEXANE COMPARTMENT FIRES 


Experiments using liquid n-hexane fires examined the effects of transient time of the 
species concentrations in the hallway and external burning. The evolution of CO, UHC, 
CO,, O, (all reported on a wet-basis) and temperature within the hallway are presented 
for the different situations studied. All experiments were performed with a 0.20 m 
diameter fuel pan, a 0.20 m inlet soffit and a global equivalence ratio of 3.0 + 0.50. Note 
that these conditions are similar to those used in experiments to produce the data showing 
the effect of the exit soffit (upper-layer in the hallway) on the post-hallway yields seen in 
Fig. 3.13, Fig. 3.14 and Fig. 3.15. The 25 sampling locations within the facility were 
shown in Fig. 2.10 and are denoted by the lower case letters on the contour plots to 
follow. Note the y-axis in all contour plots corresponds to the center of the compartment- 


hallway opening. 


Since gas sampling could only be performed at one location during an experiment, 
25 different experiments were performed to generate the data to map the facility. The 
growth rate in each of the 25 fires was slightly different. The experiments had to be 
aligned with one another for the averaged data from each experiment to represent the 
same period in the fire growth. Carbon monoxide concentrations were measured to be 
highest just after the onset of flashover ride (termed post-flashover) of the compartment 
fire. The fire growth parameter, G, defined in Eqn. (2.15) was utilized to determine the 
beginning of the post-flashover period inside the compartment fire. Post-flashover was 


assumed inside the compartment when G was approximately zero. 


The amount of time after flashover for the hallway species concentrations to reach a 
steady-state is termed the transient time. Species concentrations were considered to reach 
a steady-state when a less than 5% change occurred over a 30 second period. The 


transient time was determined by the conditions at the exit of the hallway. Experiments 
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which contained an exit soffit resulted in a short transient time, while longer transient 


times were achieved by blocking the bottom portion of the hallway exit. 
4.3.1 Short Transient Time with No External Burning 


Experiments were performed with a 0.04 m* opening and an upper-layer depth of 
0.60 m (exit soffit of 0.60 m) to investigate the evolution of species concentrations and 
temperatures when no external burning occurred in the hallway. The nondimensional 
upper-layer depth for these conditions was y=1.7. The experimental results previously 
seen in Fig. 3.15 indicate no external burning for this value of y. With the 0.60m exit 
soffit present, species concentrations and temperatures were anticipated to have a short 


transient time. 


The 25 experiments necessary to map the hallway concentrations were aligned with 
on another using the fire growth parameter, G. Compartment temperature data from 10 
experiments was used to verify the validity of the alignment method. The temperatures 
during the post-flashover period of the compartment fires are shown in Fig. 4.2 to be 
aligned with one another. Noting the different times of ignition in Fig 4.2, the alignment 
compensates for differences in the initial growth of the compartment fire. With the data 
aligned, a description of the temporal and spatial evolution of the species concentration 


and temperature in the hallway upper-layer is given. 


Approximately 12-16 seconds before flashover inside the compartment, 1.5 % CO 
and 1.5% UHC, see Fig. 4.3 and Fig. 4.4, respectively, were measured in the 
compartment upper-layer. The CO and UHC levels were measured to rapidly be diluted 
as the gases were transported across the hallway. A CO, concentration of 7.0% was 
measured to be uniform throughout the hallway upper-layer, see Fig. 4.5. An O, 
concentration was also measured to be uniform within the hallway at a level of 8.5%, see 


Fig. 4.6. A plot of the temperature in the hallway is seen in Fig. 4.7. The gases were 
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Adjusted Time Using Fire Growth Parameter, t, [sec] 


The temperature 0.05 m below the ceiling inside the compartment in 10 
different n-hexane fire experiments with a 0.04 m’ opening, a 0.60 m exit 
soffit, and no external burning in the hallway. The experiments were aligned 
with one another in time using the fire growth parameter. 
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measured to enter the hallway at temperatures of 780 K (see Fig. 4.7), but cooled rapidly 
to approximately 450 K by the end of the hallway. 


Post-flashover was reached in the compartment 16-20 seconds after the onset of 
flashover. High concentrations of CO and UHC, see Fig. 4.8 and Fig. 4.9, respectively, 
were measured in the gases flowing across the hallway and along the wall opposite the 
compartment (the path of the bulk flow). The compartment fire gases entered the hallway 
with 2.4% CO, see Fig. 4.8. Carbon monoxide concentrations of 1.6% or greater were 
measured in gases flowing across the hallway and traveling down the side of the hallway 
opposite the compartment. Only 0.9% CO was measured in the gases on the 
compartment side of the hallway which is nearly half the levels measured on the opposite 
side of the hallway. The measured UHC concentrations are shown in Fig. 4.9 to be non- 
uniform in the hallway. Over 2.0 % UHC was measured on the side of the hallway 
opposite the compartment and 1.3% on the compartment side of the hallway. Plots of the 
measured concentrations of 7.0% CO, and 7.5% O, are seen in Fig. 4.10 and Fig. 4.11, 
respectively, to be uniform within the hallway upper-layer. A plot of the measured 
hallway gas temperatures is seen Fig. 4.12. The temperatures are slightly higher than the 


levels 12-16 before flashover, but spatially vary in approximately the same manner. 


At approximately 40-44 seconds after the onset of flashover in the compartment, all 
measured species concentrations were uniform within the hallway upper-layer. The 
hallway species concentrations reached a steady state 56 seconds after CO levels of 
greater than 1% were measured in the gases entering the hallway. The upper-layer gases 
0.05 m below the ceiling were measured to contain concentrations of 1.9% CO and 2.9% 
UHC along the length of the hallway, see Figs. 4.13 and 4.14. The gases were also 
measured to contain concentrations of 6.5% CO, and 7.2% O, along the hallway length, 
see Fig. 4.15 and 4.16, respectively. From the plot of the gas temperature (see Fig. 4.17), 
temperatures are seen to cool from 850 K at the window to less than 500 K by the end of 
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Figure 4.13 The wet CO concentration distribution 0.05 m below the ceiling 40-44 seconds after flashover in n-hexane 
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the hallway. The species concentrations and temperatures within the hallway upper-layer 
at 40-44 seconds after flashover were measured to be relatively constant (+5%) until the 


end of the post-flashover period at 112-116 seconds. 


The hallway upper-layer is considered to be in steady-state between 40 and 112 
seconds after flashover since minimal changes in the species concentrations and 


temperatures were measured during this time period. 


The exposure time, 


1.587 
|" ent a 


where; 


ley - exposure time, [minutes] 


COHb - carboxyhemaglobin level in the blood, [%] 

CO - level of carbon monoxide exposure, [ppm], 
necessary for a particular concentration of CO (in ppm) to result in a fatally high 
concentration of carboxyhemaglobin (COHb) in the bloodstream can be estimated by 
assuming a 50% COHb level causes death (Peterson and Stewart, 1970). The CO 
concentration (1.9% or 19,000 ppm) measured in the hallway upper-layer gases during 


steady-state was estimated to result in death after 3.25 minutes of exposure. 
4.3.2 Short Transient Time with External Burning 


The post-hallway CO yields were measured to be as much as 25% greater than 
measured in-compartment yields (see Figs. 3.13 and 3.14) for experiments with a deep 
hallway upper-layer (1.3<y<1.7) and external burning. In-hallway sampling was 
performed for experiments with a 0.12 m’ opening and a 0.60 m exit soffit (0.60 m upper- 
layer depth and y=1.36) because these experimental conditions resulted in the highest 
post-hallway CO yields. The measurement species concentrations and temperatures in 
the hallway was conducted to investigate the effects of external burning on the temporal 


and spatial distribution of CO. Again, the presence of the 0.60 m exit soffit is anticipated 
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to result in a short transient time for the species concentrations and temperatures in the 


hallway. 


The alignment of the 25 experiments was again performed using the fire growth 
parameter. Using 10 experiments as a representative set of the 25 experiments performed 
for this case, the compartment temperatures during the steady-state of the compartment 
fire are observed in Fig. 4.18 to be aligned with one another after shifting them in time 
using the fire growth parameter. The contour plots of measured species concentrations 


and temperatures in the hallway are presented in the remaining paragraphs of the section. 


Gases containing 1.3 % CO and 1.0% UHC were measured in the compartment 
upper-layer, see Fig. 4.19 and Fig. 4.20, respectively, 12-16 seconds before the onset of 
flashover inside the compartment. The CO and UHC concentrations were measured to 
rapidly be diluted upon entering the hallway. The hallway upper-layer was measured to 
contain a uniform concentration of 8.0% CO,, see Fig. 4.21. Approximately 7.0% O, 
(see Fig. 4.22) was measured in the upper-layer gases within the hallway. From the plot © 
of measured gas temperatures in the hallway shown in Fig. 4.23, the gas temperature is 
seen to be less than 750 K at the opening. Temperatures were measured to cool to less 
than 500 K by the end of the hallway. The measured results shown in Figs. 4.19 - 4.23 
are comparable to those measured in experiments with no external burning at this stage of 


the fire (see Figs. 4.3-4.7). 


The beginning of the post-flashover period in the compartment occurred 
approximately 16-20 seconds after the onset of flashover. The species concentrations and 
temperatures of the gases in the hallway begin to change rapidly in a manner similar to 
that seen in experiments with no external burning, see Figs. 4.8-4.12. Non-uniform 
concentrations of CO and UHC were measured to be present in the hallway. The gases 
traveling on the side of the hallway opposite the compartment were measured to contain 
CO levels of 1.6%, see Fig. 4.24, while gases on the compartment side of the hallway 


contained 0.9% CO. The gases on the side of the hallway opposite the compartment also 
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Figure 4.18 The temperature 0.05 m below the ceiling inside the compartment in 10 
different n-hexane fire experiments with a 0.12 m* opening, a 0.60 m exit 
soffit, and external burning in the hallway. The experiments were aligned 
with one another in time using the fire growth parameter. 
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contained higher UHC concentrations (1.3%) compared with the UHC levels (0.7%) in 
the gases on the compartment side of the hallway, see Fig. 4.25. Plots of the measured 
CO, and O, concentrations are seen in Figs. 4.26 and 4.27. Approximately 8.7% CO, 
and 4.8% O, were measured to be relatively uniform within the hallway upper-layer. The 
gas temperatures were measured to be slightly higher in magnitude (see Fig. 4.28) 
compared with levels at the 16-20 before flashover (see Fig. 4.23). 


Just prior to the onset of external burning, all species concentrations were measured 
to be uniform within the hallway 44-48 seconds after flashover, see Figs. 4.29-4.32. The 
steady-state is observed 60 seconds after gases containing greater than 1% CO were 
measured entering the hallway. The plot of the measured CO and UHC concentrations is 
seen in Figs. 4.29 and 4.30, respectively. Upper-layer gases downstream of the 
compartment were measured to contain 2.3% CO and 2.5% UHC (see Fig. 4.29 and 4.30, 
respectively). Concentrations of approximately 8.5% CO, and 3.6% O, were also 
measured to be uniform within the upper-layer, see Fig. 4.31 and Fig. 4.32, respectively. 
Gas temperatures were measured to be 825 K at the opening (see Fig. 4.33), but were 


only 500 K at the exit of the hallway. 


To this point, the species concentrations and temperatures of the gases within the 
hallway are similar to that observed in experiments without external burning. This is 
expected since the global equivalence ratio and the mass flow of the gases into the 


hallway was similar in the two experiments. 


External burning occurred in the hallway 50-60 seconds after flashover. The burning 
was difficult to observe, since it was almost entirely surrounded by the dark upper-layer 
of combustion gases already present in the hallway. Burning gases were observed 
crossing the hallway and traveling along the wall opposite the compartment for 


approximately 2.0 m. External burning was observed to be attached to the vent. 
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The effect of external burning on species concentrations and temperatures is evident 
72-76 seconds after flashover. The CO concentration was measured to decrease slightly 
to 2.1% (see Fig. 4.34) on the side of the hallway opposite the compartment where the 
burning was visually apparent. The concentration of CO on the compartment side of the 
hallway was, however, measured to be slightly higher at 2.5%. The effect of external 
burning on the UHC concentration in the hallway is seen in Fig. 4.35. The gases along 
the side of the hallway opposite the compartment were measured to contain low UHC 
concentrations (1.5%). The CO, concentration of 9.7% was measured in the gases where 
burning was observed, see Fig. 4.36. The CO, concentration was measured to be nearly 
1.0% higher than those levels on the compartment side of the hallway. The O, 
concentration in upper-layer gases was measured to gradually increase from 1.2% near 
the compartment to 3.6% by the end of the hallway, see Fig. 4.37. Note that the 
concentrations of UHC, CO, and O, measured on the compartment side of the hallway 
are approximately equal to the species levels just prior to external burning 44-48 seconds 
after flashover, see Fig. 4.30-4.32. The oxidation of mainly UHC in the hallway is 
observed in Fig. 4.38 to result in only a 50 K rise in the upper-layer temperature with 
levels near the compartment approximately 750 K. These temperatures were typically 
measured above the flame front and are, therefore, cooler than the flame temperature. 
Slow response times of the aspirated thermocouples, approximately 1 second, prevented 


the measurement of the turbulent flame temperature with flames at the thermocouple. 


The high CO concentrations was attributed to the oxidation of the UHC. According 
to Westbrook and Dryer (1984), UHC have a higher rate of reaction with radicals (i.e. H, 
OH, H,O, O) compared with CO. Also, the oxidation of UHC forms additional CO. 
Through the law of mass action, the ratio of the UHC concentration to the CO 
concentration must be low (i.e. low UHC concentration or high CO concentration) for the 
oxidation reactions of CO to become competitive. In these experiments, the oxidation of 


UHC and CO is believed to occur at or very close to the flame front where temperatures 


137 


in the Hallway 


100NS 


Concentrat 


1€S 


Evolution of Speci 


Chapter 4. Results 


‘BUIUING [LUID}X9 PUL JJJOS XO UI 090 ‘“IJOs JOTI WI QZ‘ ‘Sutuedo ut Z1°O B WIM sjuouttiadxa arty 
OURXOY-U Ul JOAOYSE]J JoIJfe SpUOdES 9/-Z/ BUTIIOD By} MOTEG UW CSQ'Q UONNGIsIP uOTeNUDDUOD OD 19M OUT, PE'p NIy 


Ovo 
09°0 
08'0 
00°} 
02h 
Or'l 
09°t 
08"t 
00°¢ 
02°2 
10) 4 
09'2 
08'e 
00°€ 
02°€ 
Ove 
09°€ 
08’e 
00'r 
0c'v 
Ov'y 
%OD |2A97 


¥~rmn-TrTOuUuwadodvMAKe HARON TAN 


a 


[uw] *X *X.LITIOVA SHL NMOC FONVLSIC 


€ 


AemyeH 


Cc 


[w] ‘A ‘X.LITIOVA AHL SSOUOV AONVLSIG 


138 


in the Hallway 


10NS 


Evolution of Species Concentrat 


Chapter 4. Results 


‘SUTUING [VUIO}XO PUL IJJOS }IXO UW QQ°O ‘WJOS JO[UI WI QZ'Q ‘Sutuado 2H 710 & YIM syuouttodxa ally 
SUBXOY-U Ul JOAOYSe]J JayJe spuodas 9/-7/, BUITIOD ay) MOTAqG UW CO'C uonnamnsip UONVNUIDUOD DHM 10M OUL Sep WNsIy 


000 
0€"0 
090 
060 
oz"t 
os"l 
08"t 
ore 
ovz 
02% 
00'€ 
oe'e 
oe 
06° 
Oz" 
Os"p 
08'r 
ols 
or's 
02'S 
00°9 

% OHN 


<r-TOuUuwaonmnctoawnrnownwtnn - 


oul 


j@Ae7 


[ur] 


¢ 


xX 


‘ 


ALITIOVd FHL NMOd FONVLSIC 


€ 


AemyeH 


c 


juounrIeduio7) 


O'l- 


S'0- 


S‘} 


[wy] ‘A “ALITIOVA AHL SSOUOV AONVLSIA 


139 


in the Hallway 


10nS 


f Species Concentrat 


ion o 


Evolut 


Chapter 4. Results 


“SUIULING [BUIO}XO PUL IJJOS XO W (9'0 “WJOS yo[UT W QZ'O ‘3uruado 2 Z1'0 & YM syuourtiodxa oly 
OURXOY-U UI JOAOYSE]J JoIJe Spuodas 9/-7/, BUIIIAD OY} MOTIq WI CQ'G UOTNGLNSIP UONRNUDDUOD Z7_ 9M 


OL 
8 

v8 

88 

2'6 

9°6 

OL 

v0l 

80 

AD" 

O14 
Zh 

vzk 

sZt 

Ze 

9° 

bh 
% 209 


oI Xa —Torvrwoondtanar ow +o wv = 


jeAe7 


oe 
to 


ES 


[wi] “x “ALITIOVA SH.L NMOd FONVLSIC 


€ 


Aem[eH 


c 


juounreduios 


S'0- 


00 


O'l 


UL OC p wns 


[wi] “4 *X.LITIOVA AHL SSOUOV AONVLSIG 


140 


in the Hallway 


10NS 


Concentrat 


1€S 


Evolution of Speci 


Chapter 4. Results 


‘SUIUING [BUII}IX9 PUP IJJOS ITX9 WI (9'( ‘IJOS Jo[UT WI QZ'C ‘SuTuado UW Z1'0 & YIM syuouttodxe 
ONL} OUBXOY-U Ul JOAOYseTJ 19]Je SpUOdAS 9/-Z/ ZUT[I9D OY} MOTAqG UW CQ'Q UONNGNsIp UoMeNUSDUOD 7H 190M DUT, LE'p INI 


0 
9°0 
Zh 
8h 
v2 

€ 
oe 
Zp 
8p 
ys 

9 
99 
ZL 
gL 
v8 

6 
9°6 

Z'0l 

80k 

vb 
Zt 
% ZO |eAe7 


4 ¥Yrn-Touwuwao0oddtoawnownrnonwtonrane 


[wi] *X “ALITIOVA SHL NMOd FONVLSIC 


€ 


Aem[eH 


c 


quounreduro; 


[w] “A ‘XTITIOVA FHL SSOUWOV AONVLSIG 


14] 


in the Hallway 


10nS 


Concentrat 


1€S 


Evolution of Speci 


Chapter 4. Results 


‘BUIUING [BUIO}XO PUL IIJJOS 1X9 WI 9° ‘IIJJOS JO[UI W QZ” ‘Suluado 2 Z1°0 B YIM sjuouttIodxa 
OJ UBXOY-U Ul J9AOYSeIJ Ja\je spuodses 9/-7/ BUI[IZD dy) MO[9Q WI CQ'Q UONNGsIp oimpesoduiay sy gE'p aNsIy 


00°00S 
00°S¢S 
00°0SS 
00°SZS 
00°009 
00°S¢9 
00°0S9 
00°SZ9 
00°002 
00°S2Z 
00°0SZ 
00°SZZ 
00°008 
00°S28 
00°0S8 
00°SZ8 
00°006 
00°S26 
00°0S6 
00°SZ6 
00°0001 


A dW3l 18A87 


oY*¥rn-Ttrouwaoddoaownreonwtrarae 


[w) “x “A LITIOVA AHL NMOC FONVLSIC 


€ 


AemeH 


re 


yuounredwio, 


O'l- 


30° 


[ut] ‘A ‘ALITIOVA AHL SSOWOV AONVLSIG 


142 


Chapter 4. Results: Evolution of Species Concentrations in the Hallway 


are greater than 1100 K. Temperatures were typically measured above the flame front, 
and were estimated to be less than 800 K within 3 seconds after entering the hallway. No 
oxidation is believed to be occurring above the flame front where the temperatures were 
measured since a residence of at least 6 seconds is necessary for UHC or CO oxidation to 
occur at 800 K (Pitts, 1992). The deep upper-layer in the hallway reduced the 
concentration of O, in contact with the gases entering the hallway. From Westbrook and 
Dryer (1984), limited O, results in less radicals available for the oxidation of CO and 
UHC. From the data, it appears that the majority of the radicals were utilized to reduce 
the UHC concentration. With the measured increase in CO, concentrations, CO 
oxidation was occurring in the hallway. The CO concentration was, however, only 


reduced to 2.0-2.4% before the heat losses prevented further reaction from occurring. 


The species concentrations and temperatures in the upper-layer gases were measured 
to be uniform within the hallway 96-100 seconds after flashover. The CO concentrations 
were measured to be approximately 2.0%, see Fig. 4.39. This CO level was estimated to 
cause death with less than 3 minutes of exposure. The UHC concentration in the upper- 
layer gases was measured to be reduced from 4.0% to 1.2% after the gases had traveled 
across the hallway, see Fig. 4.40. Downstream of the compartment, the UHC 
concentration was measured to be 1.0%. Upper-layer gases were measured to contain 
9.9% CO, and 1.8% O, along the length of the hallway, see Fig. 4.41 and 4.42. The 
external burning resulted in gas temperatures of approximately 900 K at the opening, see 


Fig. 4.43. These temperatures were measured to cool to 700 K by the hallway exit. 
4.3.3 Long Transient Time with External Burning 


Hallway species concentrations and temperatures were mapped (see Fig. 2.10) 0.05 
m below the ceiling in experiments with a 0.04 m? opening and the bottom portion of the 
hallway exit blocked. The blockage resulted in a 0.20 m high, 1.22 m wide opening at 
the top of the hallway exit. The upper-layer gases containing high combustion product 


concentrations exited the hallway through the opening at the exit and were collected in 
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the fume hood. As a result, the species concentrations were not measured to reach a 
steady-state concentration until the end of the post-flashover period. Effects of external 
burning on species concentrations and temperature were also measured. The location of 
the sampling, previously noted in Fig. 2.10, is denoted by the lower case letters on the 


contour plots. 


The fire growth parameter was used to align the compartment fire post-flashover 
period of 23 experiments. The alignment was verified using the temperature 0.05 m 
below the ceiling inside the compartment from ten experiments. The compartment 
temperatures are shown in Fig. 4.44 to coincide with one another during the steady-state 
period. With the experiments aligned, the temporal and spatial variation in the species 
concentrations and temperatures measured in the hallway are described in the paragraphs 


to follow. 


The gases entering the hallway were measured to contain 2.4% CO and 3.9% UHC 
approximately 12-16 seconds after the onset flashover inside the compartment, see Fig. 
4.45 and Fig. 4.46. The flow of gases across the hallway diluted CO and UHC to less 
than 1.0 %. With no external burning in the hallway, CO, concentrations were measured 
remain low at 6.5-7.5% as shown in the plot in Fig. 4.47. A plot of the O, concentration 
is seen in Fig. 4.48. Approximately 9.8% O, was measured downstream of the 
compartment. Figures 4.47 and 4.48 show a region (y=0.9 and x=1.4-4.6 m) of low CO, 
concentration and high O, concentration not present in the two cases discussed 
previously. This is attributed to the presence of a more shallow upper-layer (6~0.45 m) 
in these experiments compared with the previous two cases (6~0.60 m). It is believed 
that lower layer gases containing high (18-21%) O, concentrations were being pulled into 
the upper-layer by the vortical structure traveling along the corner of the wall (opposite 
the compartment) and the ceiling, see Fig. 4.1. The addition of the lower-layer gases in 
this region dilute the upper-layer gases resulting in lower CO, concentrations. The high 


O, concentrations at the bottom right corner of the hallway are attributed to ambient air 
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Figure 4.44 The temperature 0.05 m below the ceiling inside the compartment in 10 

different n-hexane fire experiments with external burning, a 0.04 m? 
opening, and the exit of the hallway blocked leaving a 0.20 m high, 1.22 m 
wide opening. The experiments were aligned with one another in time 


using the fire growth parameter. 
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entering the hallway through the exit opening. The temperatures were measured to be 
approximately 500 K throughout the hallway, see Fig. 4.49. This is expected since no 


external burning was observed to be occurring in the hallway at this point in the fire. 


External burning was observed to occur in the hallway for these experiments 20-24 
seconds after flashover. The external burning in the hallway was not attached to the 
compartment window, but instead burning was observed to begin after the gases had 
traveled approximately half way across the hallway. This type of burning, where the 
flame is not attached to the window, was previously observed in the experiments 
performed by Gottuk (1992a). The downstream yields were measured to be higher in the 
experiments with external burning unattached to the opening compared with experiments 
with burning attached to the opening (Gottuk, 1992a). The burning gases impinged upon 
the wall opposite the compartment and traveled in a clockwise swirling motion along the 
corner of the ceiling and the side wall, see Fig. 4.1. The flame was visually observed to 
extend down the hallway approximately 2.5-3.0 m, slightly further than observed in the 


previous cases. 


Approximately 48-52 seconds after flashover with external burning in the hallway, 
gases traveling along the side of the hallway opposite the compartment were measured to 
contain 1.7% CO, see Fig. 4.50. Only 0.8% CO was measured in the gases on the 
compartment side of the hallway. The UHC were measured to be oxidized from 3.9% to 
1.5% after the gases flowed from the opening to the opposite side of the hallway as 
shown in the plot in Fig. 4.51. External burning resulted in an increase in CO, 
concentration to approximately 10.5%, see Fig. 4.52. The O, concentrations were 
measured to fall to 1-2% range, see Fig. 4.53, another indication that external burning 
was occurring in the hallway. Gas temperatures were measured to rise significantly to 


approximately 850 K at the compartment, see Fig. 4.54. The temperature was measured 


155 


Chapter 4. Results: Evolution of Species Concentrations in the Hallway 


‘BUIUING [eUIN}x9 pur }Ixo 
Aemyjey Jo do} ye Sutusdo opim Wi ZZ] YSIY W OZ'O ‘IJOs Jo[UI W QZ‘ ‘Sutuado 2 P00 & YIM sjusuIodxo 
oJ WUBXOY-U UT JOAOYsely Joye spuosss Q[-Z] BUITIID oy) MOTAq UW CO'Q UONNGIsIP ainjesoduia) oY] 6p'p WNBA 


[w) *x “ALITIOVA SHL NMOd FONVLSIG 


v € c L 
00°00 


00°Sz2S 
00°0SS 
00°SZS 
00°009 
00°Sz9 
00°0S9 
00°S2Z9 
00°002 
00°Se2z 
00°0SZ 
00°SZZ 
00°008 
00°Sz8 
00°0S8 
00°Sz8 
00°006 
00°Sz6 
00°0S6 
00°SZ6 
00°0001 


A dWal |8Ae7 


quounredui0os; 
Ww 
% 


Aem[eH 


4IO¥rn- TrTouwacd dt aowmwnrnonwtooanrae 


[w] ‘A ‘X.LITIOVA AHL SSOUOV AONVLSIA 


156 


“SUIUING [eUIO}X9 pure yIxo AeMeY 
Jo do} ye Sutusdo opim wi ZZ] YysIY WI QZ'O ‘IYJos yo[UI WE QZ'O ‘Sulusdo Wi p0'0 & YIM syuoUtodxe ory 
OUBXOY-U Ul JOAOYSeIJ JOYJB SpUOdES Z7C-gy BUITI90 oY) MO9q WI CO'CO uonnqmsip uONe.NUIDUOD OD 10M YL OG p WNsI4 


[w) *x “ALITIOVA HHL NMOd AONVLSIA 


v € re L 
Ov'0 
09°0 
08°0 
00°} 
Och 
Ov'l 


Evolution of Species Concentrations in the Hallway 


. 
. 


Chapter 4. Results 


09°+ 
Os'l 
00°¢ 
0c'e 
Ov'e 
09'e 
08'e 
00'e 
02°€ 
Ove 
0g'€ 
o8'e 
00'F 
Oc'v 
Orr 
% OD 


4o ¥r-TtTouwaddmdtoaoandnonwtoraran ee 


joAe7 


AemyeH 


juounreduiog 


S0- 


Te) 
ny 
[w] ‘A *X.LITIOVA AHL SSOUOV AONVLSIA 


= 
an 


< 
- 


157 


Chapter 4. Results: Evolution of Species Concentrations in the Hallway 


00°0 
0€'0 
09°0 
06°0 
Oc'l 
ost 
08"} 
Ol’? 
0) 
OL'c 
00°€ 
Oc’e 
09°€ 
06°E 
02'r 
OS'p 
O8'r 
OLS 
Ov's 
O'S 
00°9 
% DHN 


jOAe7 


oyYrn-TtTouwaoadgdtiononrnowonwt nae 


[ut] *x *A.LIMIOVA HL NMOC FONVLSIG 


€ 


AemyeH 


| “SUIUING [eUIO}Xx9 pue yx AVMI[eY 
Jo do} ye Buruedo opim WI ZZ] YSIY WE QZ'C ‘IJos yo[UI W CQZ'Q ‘Suruedo 2 70° & WIM syuouttiadxe arty 
sueXdY-U Ul JOAOYseTJ JajJe SpUOSIS Z7¢-gp SUTIISD SY} MOTOG WI CQ'Q UONNGINSIP UONeNUDDUOD DH 19M SUL [Sp ons 


C 


juounreduios) 


S0- 


S'k 


158 


[wi] ‘A ‘X.LITIOVA AHL SSOMOV AONVLSIG 


“SUIUING [eUIO)X9 pur yIxo ABMTTeY 
Jo do} ye Suruodo opim wi ZZ] YsIY W QZ'O ‘WWJos yo[UI W CZ‘ ‘SuTUado {WI p00 & YM sjuauttedxe arly 
ouBXdY-U Ul JOAOYSe]J J9]Je SpUOdIS Z7C-gp BUI[IZD dy) MOTAq WI CQ'O uorNquIstp UOTJVINUIDUOD “OD 19M YT, ZC'p WNSBI 


ions in the Hallway 


Evolution of Species Concentrat 


Chapter 4. Results 


9 
v9 
8'9 
ol 
92 

8 
v8 
88 
c'6 
9°6 

OL 
vOL 
8'OL 
rm 
OLE 

cl 
v'ch 
Sch 
cel 
9'EL 

vl 


% ZOD |8Ae7 


¥~rxr- TOuUw.wao emit ownwnrnwonw MAA 


| 


[wr] 


¢ 


xX 


¢ 


ALITOVdA FHL NMOC FONVLSIG 


€ 


AemyyeH 


c 


juounreduros 


[w] “A ‘ALITIOVA AHL SSOUOV AONVLSIA 


159 


Chapter 4. Results: Evolution of Species Concentrations in the Hallway 


‘SUIUING [BUIO}XO pur 3X9 
Aemyfey Jo doy ye Sutusdo opim Wi ZZ] YsIYy W QZ'O “IJOs [UT WI QZ'9 ‘Sutuado {WI 70'0 B YIM sjuourtIedxe 
ON} OUBXOY-U Ul JOAOYSeIJ Jolye SpuodasS ZC-gp SUI[I9D SY} MO[3Q UW CQ'Q UONNQINSIP UOTeNUDDUOD 7~Q OM DUT, ESP BNI 


[ut] “x “XLITIOVA AHL NMOd AONVLSIC 
v € C L 


0 
90 
a 
8h 
v2 

e 
oe 
ab 
Sy 
ys 

9 
9'9 
ray 
BL 
v8 

6 
9°6 

a 

8‘! 

a 

ZI 

% ZO 19Ae7 


= 
= 
1] 


iD 
Te 


quounseduio-) 
[w] ‘A ‘ALITIOVA AHL SSOMOV AONVLSIG 


Aem [eH 


a “Von. 2 OO WO OOO 10 i} fOlm SIO oN LS 
°o 
(o) 


160 


‘SUIUING [BUI9}Xo pur JIX9 
Aemypey Jo doy ye Burusdo opim wi ZZ] YsIy WI OZ'C ‘WJJOS JOTUI WI QZ'O ‘suluado 2W 0'0 & YIM sjuauttedxe 
O11J OUBXOY-U UT JOAOYsSeLJ Joye spuoosss Z¢-gpy BuUI[Is0 oy) MOTAq UI CQ'C uonnqinsip ainjeloduia} ayy, Sp wnsiy 


in the Hallway 


10nS 


. 


Concentrat 


1€S 


Evolution of Speci 


Chapter 4. Results 


00°00S 
00°S2eS 
00°0SS 
00°SZS 
00°009 
00°S29 
00°0S9 
00°SZ9 
00°002 
00°S¢Z 
00°0S2 
00°SZZ 
00°008 
00°S28 
00°0S8 
00°SZ8 
00°006 
00°S26 
00°0S6 
00°S26 
00°0001 


A dW3AL 1927 


¥~¥~rmr-TtTOuUuwwao«gadtoawndwnrnoonwtonrnne 


a 


[ur] 


6 


xX 


‘ 


ALYPTOVA AHL NMOC FONVLSIC 
re 


€ 


AemyeH 


juounsJedui07) 


O'L- 


S'0- 


00 


s'0 


[w] ‘A ‘ALITIOVA AHL SSOUOV AONVLSIG 


161 


Chapter 4. Results: Evolution of Species Concentrations in the Hallway 


to rapidly fall with distance away from the compartment to 650 K by the end of the 
hallway. 


Concentrations of 1.9% CO were measured in the gases traveling along the side of 
the hallway opposite the compartment 64-68 seconds after flashover despite the presence 
of external burning in the hallway, see Fig. 4.55. The CO concentration on the 
compartment side of the hallway was measured to be approximately 1.0%, nearly half of 
the concentration present on the opposite side of the hallway. The UHC, however, 
concentration was measured to be reduced from 3.0% at the opening to 1.2% at the wall 
opposite the compartment, see Fig. 4.56. As described in the previous section, the rate of 
reaction of UHC and radicals is higher than that of CO and radicals (Westbrook and 
Dryer, 1984). The UHC are therefore oxidized more readily than CO. The CO 
concentration in the present experiments was measured to be slightly less than in 
experiments described in section 4.3.3. The difference is attributed to the less deep 
upper-layer in these experiments. The more shallow upper-layer allows additional O, to 
mix with the hot gases promoting the production radicals for oxidation. Concentrations 
of 10.5% CO, and 1.2% O, were measured to be uniform in the hallway upper-layer, see 
Fig. 4.57 and Fig. 4.58. The high temperatures in the hallway (greater than 750 K) were 


measured at locations where external burning was observed to be present, see Fig. 4.59. 


An area of low CO concentration was measured to be present 2.0 to 4.0 m down and 
0.7-1.0 m across the hallway 80-84 seconds after flashover, see Fig. 4.60. The area of 
low concentration is at the same location where low CO, and high O, concentrations were 
measured 12-16 seconds after flashover, see Figs. 4.47 and 4.48. The gases on the side of 
the hallway opposite the compartment were measured to still contain 1.9% CO. The 
concentration of CO in gases on the compartment side of the hallway increased to 1.2%. 
The increase was attributed to gases being forced to circulate back toward the 
compartment. The distribution in the UHC concentration within the hallway is similar to 
that measured 64-68 seconds after flashover. The gases entered the hallway containing 
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Chapter 4. Results: Evolution of Species Concentrations in the Hallway 


2.3 % UHC, but contained only 1.2% UHC after flowing across the hallway, see Fig 4.61. 
The gases flowing along the side of the hallway opposite the compartment were measured 
to contain UHC concentrations of nearly 10 times lower than the UHC concentration in 
the gases entering the hallway. The external burning continued to deplete the hallway 
gases of UHC, but the CO concentration in the hallway continued to rise. The upper- 
layer of the hallway was measured to contain a uniform distribution of approximately 
10.5% CO, and 1.5% Oy, see Fig. 4.62 and Fig. 4.63 respectively. The temperatures were 
measured to be greater than approximately 775 K, see Fig. 4.64, and corresponded to 
where the external burning was visually observed in the hallway. The species 
concentrations were measured to reach a steady-state shortly after this period at 


approximately 96 seconds after flashover. 


Approximately 108-112 seconds after flashover, external burning was no longer 
occurring in the hallway. The low concentration cell of CO in the upper-layer was 
measured to still be present, see Fig. 4.65. Slightly lower levels of CO (1.5%) were 
measured in the gases along the side of the hallway opposite the compartment. This was 
attributed to the fire entering the decay period. The CO concentration on the 
compartment side of the hallway was measured to be 1.4% CO, see Fig. 4.65. The UHC 
concentration was measured to increase slightly to 0.6% downstream of the compartment, 
see Fig. 4.66. The absence of external burning was also apparent by the lower levels of 
CO, (8-9%) and higher levels of O, (4-5%) measured in the upper layer gases, see Figs. 
4.67 and 4.68, respectively. The temperatures in the hallway were measured to be 
approximately 50 K lower in magnitude after external burning has ceased in the hallway, 
see Fig. 4.69. 
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4.4 EVOLUTION IN POLYURETHANE FOAM COMPARTMENT FIRES 
4.4.1 Short Transient Time with No External Burning 


Compartment fire experiments with polyurethane foam as fuel were performed to 
investigate the effects of burning a fuel more commonly found in health care facilities 
and office buildings (i.e. upholstered furniture and bedding). Experiments were 
conducted with a 0.04 m’ opening and a 0.60 deep hallway upper-layer. This 
corresponds to a nondimensional upper-layer depth of y=1.7. As stated in section 3.4, no 
external burning was visibly apparent in these experiments. The compartment global 


equivalence ratios were, on average, 2.10 + 0.07. 


The fire growth parameter was also used in these experiments to align the post- 
flashover period of the 11 experiments performed in this study, see Fig. 2.11. The 
compartment temperatures of the experiments are shown in Fig. 4.70 to be aligned with 
one another. Note that the time scale on Fig. 4.70 is different than that of Fig. 4.2, Fig. 
4.18 and Fig. 4.44. The post-flashover period in these experiments was nearly twice as 
long as the post-flashover in the n-hexane experiments. The location of the sampling, 
previously noted in Fig. 2.11, is denoted by the lower case letters over top of the contour 


plots to follow. 


With no external burning in the hallway during these experiments, the variation in 
the combustion gases, both spatially and temporally, is expected to be similar to the 
n-hexane experiments with no external burning discussed in section 4.3.1. The presence 
of the 0.60 m exit soffit, as previously seen in section 4.3.1 and 4.3.2, resulted in a short 


transient time for the species concentrations in the hallway upper-layer. 


At 32-36 seconds after the onset of flashover inside the compartment, the 
concentration of CO was measured be 1.9% in the compartment, see Fig. 4.71. As shown 
in the plot of CO concentrations in Fig. 4.71, the CO concentration was diluted rapidly to 
less than 0.5% by the time the gases had traveled across the hallway. High 
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Figure 4.70 The temperature 0.05 m below the ceiling inside the compartment in 10 
different polyurethane foam fire experiments with no external burning, a 
0.04 m* opening, and a 0.60 m exit soffit. The experiments were aligned 
with one another in time using the fire growth parameter. 
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concentrations of CO, (11%) were measured in the gases entering the hallway, see Fig. 
4.72. These high concentrations were also measured in the gases flowing across the 
hallway and down the side of the hallway opposite the compartment. Similarly, lower 
levels of O, were measured in gases traveling across the hallway and down the side of the 
hallway opposite the compartment, see Fig. 4.73. The temperatures in the hallway just 


downstream of the compartment remained less than 525K as shown in Fig. 4.74. 


At 64-68 seconds after flashover in the compartment, 1.7% CO was measured in the 
gases across the hallway from the compartment, see Fig. 4.75. The CO concentration in 
the gases along the compartment side of the hallway, however, was measured to be 1.0%, 
approximately half the level on the opposite side of the compartment. The CO, 
concentrations were measured to te uniform at a level of 6.7 %, see Fig. 4.76. The O, 
concentrations in the gases near the compartment were measured to be approximately 
7.2% near the compartment and gradually increased in concentration as the gases flowed 
down the hallway, see Fig. 4.77. Temperatures were measured in Fig. 4.78 to remain low 


at levels less than or equal to 500 K downstream of the compartment. 


The species concentrations in the hallway reached a steady-state at approximately 
92-96 seconds after post-flashover inside the compartment. This is approximately 60 
seconds after 1% CO was measured entering the hallway. The concentration of CO in the 
upper-layer was measured to be uniform in the upper-layer at approximately 1.8%, see 
Fig. 4.79. The CO, concentration was measured to be 7.0% within the upper-layer as 
shown in Fig. 4.80. Oxygen concentration in the hallway was measured to remain 
uniform at approximately 7.0%, see Fig. 4.81. The gas nating also remained low at 
temperatures of approximately 500 K or less throughout the hallway upper-layer, see Fig. 
4.82. 
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4.5 SUMMARY OF SPECIES CONCENTRATION EVOLUTION RESULTS 


The temporal and spatial evolution of CO, UHC, CO, and O, concentrations and 
temperatures in the hallway were measured for experimental situations where high CO 
post-hallway yields were measured. The experiments were conducted to determine the 


effect of external burning, species concentration transient time and fuel type. 


In all experiments with the compartment on the side of the hallway, the movement of 
fire exhaust gases away from the compartment was non-uniform in the hallway. Gases 
first flowed across the hallway until they impinged on the wall opposite the compartment. 
The impingement of the hot gases on the wall resulted in the development of a large 
vortex on the side of the hallway opposite the compartment. The majority of the gases 
were visually (through observation of external burning) observed and measured to travel 


down the side of the hallway opposite the compartment. 


External burning was measured to have a significant effect on the UHC, CO, and O, 
concentrations, but a minimal effect on CO concentrations in the hallway. With external 
burning in the hallway, 1.9-2.1 % CO were measured in the upper-layer of the hallway. 
The UHC were, however, measured to decrease to 25% of the UHC concentration exiting 
the compartment. The decrease in the UHC concentration more rapidly than CO 


concentration agrees with the results reported by Westbrook and Dryer (1984). 


The time before the species concentrations reach a steady-state (transient time) in the 
hallway varied depending on the blockage present at the exit of the hallway. With a 
0.60 m soffit, the transient time was approximately 60 seconds. Longer transient times 
were achieved by blocking the bottom portion of the hallway. This resulted in a transient 
time of 96 seconds. Once a steady-state was reached in the hallway, CO levels were 
approximately 2.0%, a level which could result in death with about 3 minutes of 


exposure. 
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Both n-hexane and polyurethane foam fires were burned inside the compartment. 
The amount of CO transported to remote locations was measured to be independent of the 


type of fuel burned in the compartment. 
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CHAPTER 5 


DISCUSSION: ANALYSIS OF RESULTS AND 
MODEL DEVELOPMENT 


5.1 INTRODUCTION 


The goal of this research is to provide fire safety engineers with a procedure to 
estimate CO levels in building fires. Four major variables have been identified through 
this research as being critical to the prediction of CO at a location downstream of the 
burning compartment. The variables are: the stoichiometry of the compartment fire, the 
oxygen entrainment into the gases entering the hallway, the occurrence of external 
burning, and the transient time of species concentration in the hallway. In this chapter, 
models which predict the occurrence of external burning and the transient time of species 
concentrations in the hallway are developed and verified using experimental data. The 
stoichiometry and the oxygen entrainment are used in developing models for the 
prediction of post-hallway yields. A procedure is then given to determine the hazards 
(i.e. flame spread, toxic gas levels) which are present during a building fire upon the 


onset of flashover inside the burning compartment. 


5.2 PREDICTION OF EXTERNAL BURNING 


The occurrence of external burning in the hallway has been shown in Chapters 3 and 
4 to have a significant effect on the amount of CO transported away from the burning 
room. The absence of external burning resulted in post-hallway yields equal to the levels 
measured inside the compartment. The prediction of external burning within the hallway 
was attempted using the ignition index, given as Eqn. (2.16), defined by Beyler (1984). 
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When the ignition index rises above 1.0, the amount of energy available for ignition 
exceeds the amount of energy necessary for ignition. Ignition is, therefore, predicted to 


occur when I.I.>1.0 assuming a pilot ignition source exists. 


In chapter 3, experimental results were presented where an ignitable mixture was 
produced by a compartment fire, but ignition did not occur in the hallway. The absence 
of external burning was attributed to the presence of a deep oxygen deficient upper-layer 
in the hallway. The ignition of hallway combustion gases was determined to be a 
function of the nondimensional upper-layer depth, y=6/z. With a nondimensional upper- 
layer depth less than 1.7, it was experimentally determined that a pilot ignition source 
was present to ignite the hallway gases. Above a value of 1.7, the oxygen deficient 
upper-layer in the hallway was deep enough to prevent the entrainment oxygen into the 
compartment fire exhaust gas jet. This was hypothesized to inhibit the promotion of 


oxidative radicals and to prevent the ignition of the hallway gases. 


The proposed method for the predicting the ignition of gases entering the hallway has 
two parts. First, determine whether a pilot ignition source is present by calculating 
whether the nondimensional upper-layer depth is less than 1.7. If a pilot source is 
present, calculate the ignition index (I.I.). If I.I. is greater than 1.0, then the ignition is 


predicted to occur. 


The method is verified using numerous experimental results which contained various 
nondimensional upper-layer depths. The gases for the calculation of the II. were 
sampled in the middle of the compartment opening, at location “b” in Fig. 2.10. 


The results of the application of the model to predict external burning in an 
experiment is tabulated in Table 5.1. The temporal variation in the ignition index for an 
experiment y<1.7 and one with y>1.7 are shown in Fig. 5.1 and 5.2, respectively. These 
plots are a representative sample of the data presented in Table 5.1. The experiment 


where external burning’ occurred contained a _é 0.12 m? opening, 
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Table 5.1 The use of the ignition index and the nondimensional upper-layer depth to 
predict the occurrence of external burning in the hallway. 
Bottom of External of External 


Opening Upper-Layer 
Size, m Depth, 5, m 
Opening, z, m Burning, t,sec | Burning 


0.08 0.16 1.56 319 
0.04 0.16 1.56 317 
0.12 0.24 1.04 307 
0.04 0.36 1.67 N/A 
0.04 0.30 0.16 1.88 N/A N/A* 
0.04 0.35 0.16 2.19 N/A N/A** 


[Rees ava ws fran eo cr 


*[.J..>1.0 but no ignition occurred since y>1.7. 
**] 1.<1.0 and no ignition occurred. 


LI. at Start 


Start of 


Distance to y=5/z 


N/A* 


o=2.70 and y=1.04. The experiment where no external burning occurred in the hallway 


contained a 0.04 m? opening, ¢=3.0 and y=1.88. 


In the experiment with external burning, the ignition index predicted the ignition of 
the hallway gases at 307 seconds, see Fig. 5.1. The value of y was also determined to be 
less than 1.7, signifying the existence of a pilot source to ignite the hallway gases. As 


predicted, external burning occurred in the hallway at approximately 307 seconds. 


The experiment without external burning had a y=1.88. A pilot source for ignition 


was, therefore, predicted not to exist. The magnitude of the ignition index is seen in 
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Start of External Burning 


Ignition Index, I.I. 


0 100 200 300 400 500 


Time, t, [sec] 


Figure 5.1 The ignition index for an experiment with external burning in the hallway. 
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Ignition Index, I.I. 


0 100 200 300 400 500 
Time, t, [sec] 


Figure 5.2 The ignition index for an experiment with no external burning. 
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Fig. 5.2 to rise above 1.0 at approximately 300 seconds. As predicted by y, no ignition 
source was present and external burning did not occur in the hallway. 

This analysis suggests that the ignition index should, therefore, be used in 
combination with the nondimensional upper-layer depth, y, to predict the occurrence of 
external burning. External burning is expected occur when 

e the nondimensional upper-layer depth, y, is less than 1.7 and 

e the ignition index, I.I., rises above 1.0. 

No external burning is expected occur in the hallway with y21.7. This is attributed to the 
production of oxidative radicals being inhibited by the extremely deep oxygen deficient 
upper-layer present in the hallway. The ignition index in this case gives information on 


whether (I.I.>1.0) or not (I.I.<1.0) the gases entering the hallway are ignitable. 


5.3 PREDICTION OF SPECIES CONCENTRATION TRANSIENT TIME 


The yield of species formed inside the compartment was found to correlate well with 
the global equivalence ratio of the compartment fire (Beyler, 1984, Tewarson, 1984 and 
Gottuk, 1992). Correlations for the estimation of species levels in the space adjacent to 
the compartment are formulated in section 5.4 and are based on the post-hallway species 
yields. The limitation of using a correlation based on the post-hallway yields is that the 
species levels need to be uniform within the hallway upper-layer. The results in Chapter 
4 suggest that yield correlations are valid after the gases have reached a steady-state 
concentration in the hallway. A well-stirred model is developed in this section to provide 


an estimate the transient time for the gas concentrations in the hallway upper-layer. 


The hallway is modeled as a well-stirred volume as seen in Fig. 5.3. Gases 
containing high CO concentrations are assumed to enter the volume and instantaneously 


mix within the volume. From the assumption of instantaneous mixing, the concentration 
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Ci) 
D 


Figure 5.3. The well-stirred volume model for the estimation of species concentration 
transient time in the hallway. 


of the gases exiting the volume is equal to the gas concentration within the volume. The 


variation in the concentration of CO, C(t), inside the volume with time is, 


aD)  PajaP2c 
= +5) C,, (5.1) 


with the initial condition being, 
C(0)=0. (5.2) 


The volumetric flow rate of the gases into the volume, D, the volume, V, and the initial 


concentration of CO in the volume, C;, are all assumed to be known. 


The concentration of CO as a function of time is therefore: 
C(t)= CG 1- exp( 21) : | (3230 


Solving Eqn. (5.3) for t, 


—Ayaif? ] C,-C(t) (5.4) 


bran = (7 ) fears ea [ 
. D\ Thats] Teomp CG; 


and substituting the volume of the hallway upper-layer for V, an expression to estimate 
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Table 5.2 The conditions in experiments with both short and long steady-state times in 
the hallway. 


Transient Transient Depth of 
Timesines Time (from Layer in | D,m’/s Roe RO TK 
sec Ee ; Hallway, 
6, m 


ee A 
[| ie «toase** [0.45 | 0.136 | 950 | 350 
*Exit soffit at end of hallway. 

**Bottom blockage at end of hallway. 

***Determined using effective hallway length. 


the transient time of the gas concentrations in the hallway is given. The flow rate of the 
high temperature compartment fire gases entering the hallway is known. A density 
correction factor for the flow rate in the hallway is added in the denominator of Eqn. (5.4) 


so that the flow rate is representative of the average flow rate in the hallway. 


The transient time of experiments reported in Chapter 4 is estimated using Eqn. (5.4). 
The calculations assume that the gases entering the hallway contain 3.5% CO (C=3.5%), 
and that the upper-layer gases at steady-state contain 2.0% CO (C(t) =2.0%). All times 


are referenced from flashover (i.e. 60 seconds is 60 seconds after flashover). 


A transient time of 48 seconds was calculated using the conditions present in 
experiments with an exit soffit connected to the hallway shown in Table 5.2. This is 


within 20% of the transient time determined from the experimental data. 


The experiments with a blockage at the bottom the hallway exit simulated toxic gas 
movement in a hallway longer than 5.18 m. The effective length of the hallway was 
determined using the time required for gases containing 1.0% CO to circulate back to the 


compartment opening after traveling down the hallway. Approximately 24 seconds was 
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required for gases to circulate in experiments with an exit soffit (5.18 m long hallway). 
With the blockage at the bottom of the hallway exit, approximately 64 seconds was 
required for the gases to circulate back to the opening. The effective length of the 
hallway was determined by multiplying the actual length of the hallway (5.18 m) by the 
ratio of the bottom blockage circulation time and the exit soffit circulation time. The 
assumption here is that the velocity of the gases in the two cases is approximately the 
same. The effective hallway length was determined to be 13.81 m. Using the effective 
hallway length, the transient time was calculated to be 82 seconds. The calculated time is 
within 15% of the experimental time of 96 seconds. 


The transient time of the species concentrations in the upper-layer was calculated to 
within 20% of the experimentally determined value. The transient time calculation 
provides a method for determining when the species concentrations in the hallway have 
reached a steady-state. After such time, post-hallway species yield correlations can 
accurately be used to determine the species levels transported away from the 


compartment. 


5.4 CORRELATIONS FOR THE PREDICTION OF CO 


The post-hallway combustion product yields were shown in Chapter 3 to be strongly 
dependent on the occurrence of external burning in the hallway, the oxygen entrainment 
into the jet exiting the compartment and the stoichiometry of the compartment fire. As 
mentioned in Section 5.2, the absence of external burning results in post-hallway yields 
equal to compartment yields. Therefore, the post-hallway yields in these cases can be 
predicted using correlations based on the compartment global equivalence ratio. With 
external burning occurring in the hallway, post-hallway yields were shown in Chapter 3 
to be highly dependent on oxygen entrainment into the gases entering the hallway and the 
stoichiometry of these gases. A relationship between the post-hallway yields and the 
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oxygen entrainment into the gases entering the hallway is developed. To account for the 
stoichiometry in the experiments, a global equivalence ratio is developed for the 
compartment-hallway system. The global equivalence ratio is then used to predict post- 


hallway yields. 
5.4.1. Oxygen Entrainment into Gases Entering the Hallway 


The data in Chapter 3 showed the inlet soffit height, the compartment opening size 
and the upper-layer depth all had a significant effect on the post-hallway CO yields with 
external burning in the hallway. Variation in the soffit height and opening size was 
shown to alter the type of jet entering the hallway (1.e. ceiling or buoyant) and the jet flow 
characteristics (i.e. shape and velocity). This resulted in a variation in the mass flow of 
gases entrained into the jet. The magnitude of the entrainment was quantified using the 
entrainment function, E. The upper-layer depth in the hallway determined the 
concentration of oxygen present in the gases being entrained into the jet. The degree to 
which the jet was surrounded by hallway upper-layer gases ae quantified by the 
nondimensional upper-layer depth, y (defined as the upper-layer depth divided by the 
distance between the ceiling and the bottom of the opening). A relationship accounting 
for the effect of these parameters on the oxygen entrainment into the compartment fire 
gases is developed. The relationship is then correlated to the post-hallway CO and UHC 
yield data. 


In experiments with external burning, the CO and UHC yields were shown in Fig. 
3.13-3.15 to rise nearly exponentially with nondimensional upper-layer depth for a 
particular opening size. Assuming the oxygen entrainment was the cause of the rise, the 


oxygen entrained into the hallway gases was assumed to decrease exponentially with the 
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nondimensional upper-layer depth. The oxygen entrainment parameter, 


e = E| Xo, amolexp(-8/z))|100, (5.5) 
was developed to account for both the entrainment into the jet and the level of oxygen 


entrained into the jet. 


The data from Fig. 3.13, Fig. 3.14 and Fig. 3.15 where external burning occurred and 
the values of E shown in Table 3.4 for all three opening sizes were used to verify the 
model. The post-hallway CO yields are shown in Fig. 5.4 plotted against the oxygen 


entrainment parameter. The curve fit to the data is given by the following equation: 
Yoo = 0.42 exp(—1. 5) (5.6) 


From Eqn. (5.6), a decrease in the entrainment into the jet or an increase in the 
nondimensional upper-layer depth results in higher CO yields. This trend can also be 


seen in Fig. 5.4. 


The UHC yields are plotted against the oxygen entrainment parameter in Fig. 5.5. A 


curve fit to the data was determined to be 
Youc = 0.38 exp(—-2.0e) (5.7) 


As expected from the results in Chapter 3, the trend in the data is similar to that observed 
for CO yields. 


The oxygen entrainment parameter does not account for the amount of fuel entering 
the hallway. This is governed by the compartment equivalence ratio. The yield of CO 
and UHC produced in the compartment was shown in the study of Gottuk (1992a) to be 
constant at a compartment global equivalence ratio greater than 1.8. All experiments 
considered in this section having a compartment global equivalence ratio of 3.0 +0.5. 


Therefore, the yield of CO and UHC entering the hallway was assumed to be equal in all 
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Figure 5.4 The use of the oxygen entrainment parameter, Eqn. (5.5), to predict the 
hallway exit CO yield with external burning occurring in the hallway. Data 
from opening sizes of @-0.12 m?, A-0.08 m’, and M-0.04 m’. The solid line 
is a curve fit to the data. 
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Figure 5.5 The use of the oxygen entrainment parameter, Eqn. (5.5), to predict the 
hallway exit UHC yield with external burning occurring in the hallway. Data 
from opening sizes of @-0.12 m’, A-0.08 m*, and Ml-0.04 m’. The solid line 
is a curve fit to the data. 
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experiments. A more general approach is taken in the next section to account for the 


stoichiometry. 
5.4.2 Global Equivalence Ratio of the Compartment-Hallway System 


The formation compartment fire products has been shown to be related to the 
compartment global equivalence ratio in the studies of Beyler (1983), Gottuk (1992a) and 
others. A global equivalence ratio for the compartment and hallway is developed in this 


section. The system global equivalence ratio is then used to predict post-hallway yields. 


The system global equivalence ratio, 


m ‘fuel 
Mair. comp H Mair hall 
Olga (5.8) 


m fut =i, : 
Mair 


st 


is of the same form as the global equivalence ratio calculated for the compartment. The 
difference in the two is the air entrained into the flaming portion of the hallway 
gases, M,,, (as The air entrained into the flaming portion of the hallway gases is believed 
to control the oxidation of combustion products in the hallway. The rest of the air 
entering the hallway is merely diluting the upper-layer gases, and has no effect on post- 


hallway yields. 


The measurement of the exact amount of air being entrained into the flaming portion 
of the gases is very difficult. As a first order approximation, the measurements of the 
mass flow rate of air entrained into overventilated gases entering the hallway, see Table 
3.2 and 3.4, are used in the analysis. The air entrainment mass flow rate given in the 
tables is the air entrainment which had occurred after the compartment fire gases had 


flowed 0.60 m across the hallway. The entrainment over the entire flaming region was 
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determined by first estimating the flame extension from the compartment. The flame 


length under the ceiling was estimated using the correlation of Gross (1989) 
H 
f= 1 _ 0.3(1 -- 2)) aren: (5.9) 


where F is the flame extension, L is the free flame height, H is the height of the ceiling 
above the pool fire (height of the compartment), and F,,,,, is the portion of the flame 
extension which occurs inside the compartment (0.61 m). The free flame height was 


determined using the correlation of Heskestad (1983) 


L=-1.02d +0.230"”, (5.10) 


where d is the diameter of the pool fire and Q is the fire size. Using the calculated flame 


lengths, the air entrainment into the flaming portion of the hallway, 


ee Laer 
m.. =m,| — | — |, 5.11 
air,hall {Al Obes ) ( ) 


was determined by multiplying the measured air entrainment into the overventilated 
compartment fire gases by the ratio of the flame length and 0.6 m and a density correction 


factor. 


The model was verified using the data of Ewens (1994b) where the effects of the 
upper-layer depth on the oxygen entrainment were considered to be minimal. The post- 
hallway CO yield is plotted against the compartment equivalence ratio in Fig. 5.6. The 
compartment CO yield data from the study of Gottuk is also shown in the figure. The 
post-hallway data is completely uncorrelated to the compartment equivalence ratio. A 
correlation is, however, seen to exist between the post-hallway CO yield and the system 
global equivalence ratio, see Fig. 5.7. In fact, the majority of the data collapses to the 
data of Gottuk (1992a). This is expected since the system global equivalence ratio in the 
study of Gottuk was the compartment equivalence ratio. A small portion of the data 
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Figure 5.6. The post-hallway CO yields (@) plotted against the compartment 
equivalence ratio. The compartment yields (LJ) from the study of 
Gottuk (1992a) are also shown. 
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Figure 5.7 The post-hallway CO yields (@) plotted against the system global 
equivalence ratio. The compartment yields (LJ) from the study of 
Gottuk (1992a) are also shown. The solid line is a curve fit to the data of 


Gottuk. 
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did not collapse to the data of Gottuk. This may be attributed to an error in the estimation 
of air entrainment into the flaming portion of the gases. These points may also indicate 
that post-hallway yields are not well described by compartment data at higher global 
equivalence ratios. Assuming the error is due to a poor estimation of entrainment, the CO 
yield can be predicted from using the global equivalence ratio for the system and the 


curve fit to Gottuk’s data, 


2 el (5.12) 


The system global equivalence ratio is accurate for y<1.4. With y>1.4, CO yields were 


measured to rise above in-compartment yields making the correlation invalid. 


5.5 PROCEDURE FOR PREDICTING CO LEVELS IN BUILDING FIRES 


Using the models developed in this chapter, a procedure for estimating the CO levels 


in building fires is now possible. The procedure, see Fig. 5.8, is as follows: 
1. Determine the transient time of the gas concentrations in the adjacent space. 
2. Determine whether external burning occurs in the adjacent space. 


3. If no external burning occurs and the time after flashover is greater than the transient 


time, use the compartment yield correlations to estimate CO yields. 


4. If external burning occurs and the time after flashover is greater than the transient 


time, predict the CO yields in the hallway using: 


a) the oxygen entrainment parameter model; valid for upper-layers with y<1.7 and 


compartment global equivalence ratio greater than 1.8, or 
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Figure 5.8 The proposed procedure for estimating CO levels in building fires. 
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b) the system global equivalence ratio model; valid for any compartment global 


equivalence ratio, but for y<1.4. 


5.6 APPLICATION OF RESULTS TO FULL-SCALE SITUATIONS 


The post-hallway CO yield has been shown here to be strongly related to the 
entrainment into the gases entering the hallway. These fluid mechanics are sensitive to 
the scale of the experiment. Suggestions of how the results of the 1/4 scale experiments 


presented in this thesis can be used in full scale situations are given below. 


The scaling laws for transport of gases away from the hallway are mainly taken from 
the study performed by Steckler et al. (1986), with similar relationships developed by 
Emori and Saito (1979) and Daikou and Saito (1995). Steckler et al. (1986) investigated 
smoke movement from a burning room on the side of a hallway using salt-water 
experiments. Dimensions were scaled using the height of the room, while the 


compartment fire was scaled using the Froude number, 


Ay 


Ss 


2g. (p, -p,) 


(5.13) 


which is simply the inverse of the Richardson number, Ri. The term p, is the density of 


the gas at the source. 


It should be noted that the Ri number used in the entrainment correlation, Eqn. (3.3), 
uses a different form of the densimetric term, p,/(p;-P,), compared to p,/(p.-P,) Shown in 
Egn. (5.13). The reasoning for the difference is to account for the physics driving the 
process. The physics driving the fire plume rising to the ceiling in the compartment is the 
inertia of the plume and the buoyancy induced by the density difference in the hot plume 


and the surrounding gases. The physics driving the entrainment into the gases entering 
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the hallway is the gravitational potential of the gases, the kinetic potential of the gases, 
and the density of the hot gases relative to the entraining gases (density defect). 


To scale the present results, the dimensions of the compartment, the hallway, the — 
layer-depth, the distance between the bottom of the opening and any other physical 
distance should be scaled to the distance between floor and ceiling inside the 
compartment, h. Born out of these relationships is the nondimensional upper-layer depth 
(y=8/z) which was previously used to correlate downstream results with the O, 
entrainment into the gases entering the hallway. The compartment fire should be scaled 
relative to the Froude number. The opening size should be scaled to Acmetig to ensure 


the flow into the hallway is similar (Kawogoe, 1959). 


These scaling laws are suggested, but no experiments were performed to validate 
them. In addition, no experimental data was available to determine whether these laws 
would give similar CO yields in full scale situations. The laws should be experimentally 


tested to determine their range of validity. 
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CHAPTER 6 


ENGINEERING TOOLS FOR THE ESTIMATION OF 
CO LEVELS IN BUILDING FIRES 


6.1 INTRODUCTION 


The development of engineering tools to predict the formation, oxidation and 
transport of CO in building fires is essential for designing a more fire safe building. A 
chemical kinetic model was developed to estimate species concentrations formed inside 
the compartment. The compartment is modeled as a homogeneous reactor using the 
FORTRAN version of SENKIN. The initial species mole fractions are determined 
utilizing the conditions present inside the compartment before the onset of flashover. 
Mechanisms for both n-hexane (”-C,H,,) and methane were utilized in the study. Results 
of the model were compared to experimental measurements from n-hexane and methane 


compartment fires. 


6.2 PROGRAMS AND CHEMICAL KINETIC MECHANISMS 
6.2.1 CHEMKIN, STANJAN and SENKIN 


The FORTRAN programs CHEMKIN, developed by Kee et al. (1994) at Sandia 
National Laboratory, and SENKIN, developed by Lutz et al. (1994) also at Sandia, were 
used to perform the chemical kinetic modeling. The program STANJAN was also used in 
the modeling process. The programs were run on a UNIX based Hewlett-Packard 


workstation. 


CHEMKIN organizes the thermodynamic and chemical kinetic data to be input into 
STANJAN and SENKIN. STANJAN is used to determine the chemical equilibrium of a 
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mixture of gases. The program SENKIN is used to determine species mole fractions in a 
homogeneous reactor at a given condition. With a given temperature and pressure, 
SENKIN solves the equations for the species conservation, 

dmf,; _® jMW, 


: 6.1 
a - (6.1) 


to determine the temporal variation in the species mass fractions in the reactor. The 
values for @ ; are determined from the rate of destruction and formation for a species j. 
The rate equations of destruction and formation of species are determined using the law 


of mass action and the modified Arrenhius rate coefficients, 


k= AT* ep = (6.2) 


Constants for equation (6.2) are given for each reaction in the mechanism. 
6.2.2 Chemical Kinetic Mechanisms for n-Hexane and Methane 


A detailed chemical mechanism for n-hexane, in addition to the thermodynamic data 
for all species involved, was obtained from Lawrence Livermore National Laboratories 
(LLNL), Curran and Westbrook (1996), and was then adapted to CHEMKIN format. The 
mechanism contains 398 chemical species and 1910 chemical reactions, and is valid in 


the low to middle temperature range, 300 K to 1600 K. 


Since the mechanism only considered gaseous species, the energy necessary for the 
conversion of the liquid hexane to the gaseous state was accounted for by determining the 
difference in the adiabatic flame temperature for liquid hexane and gaseous hexane. The 
temperature difference, approximately 15 K, was accounted for in the calculations by 
setting the initial conditions in the STANJAN calculations, which model the plume, to 285 
K instead of 300 K. 
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A propane mechanism was used to model the methane compartment fire chemical 
kinetics. This mechanism was used so higher order hydrocarbons reactions could be 


considered. 


6.3 ESTIMATION OF SPECIES PRODUCTION IN COMPARTMENT FIRES 


6.3.1 Background 


Most of the knowledge of species production and behavior in the upper-layer of a 
compartment fire has come from experimental investigations. Hood experiments 
performed by Beyler (1986a., 1986b.), Toner et al. (1987) and Morehart et al. (1990) 
experimentally modeled the upper-layer environment of the compartment, and correlated 


the species produced to the global equivalence ratio (GER). 


The temperatures measured in the hood experiments were lower than those measured 
in actual compartment fires. Compartment fire experimental investigations were, 
therefore, conducted by Gottuk (1992a) and Bryner, et al. (1994) to verify the species 
levels measured in hood experiments. These investigations were also used to test the 
validity of the GER concept for predicting species levels at higher temperatures (Pitts, 
1994a). 


The success of estimating the species levels measured in the compartment upper- 
layer with chemical kinetic modeling has been limited. In the study by Morehart (1990, 
1991), the qualitative effect of temperature on the production of species in the upper-layer 
of the hood experiments was determined, and compared with his data. The chemical 
kinetic modeling of Pitts (1992) used the species levels from Morehart’s experiments to 
investigate the effects of temperature and GER on the species levels in the upper-layer 
environment. The study by Pitts (1992) showed a drastic increase in the level of CO at 
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elevated temperatures (1300 K), and related this rise to a shift in the water-gas shift 


reaction: 
CO, + H,=CO+ H,O0 (6.3) 


toward chemical kinetic equilibrium at high temperatures. The prediction of species 
levels in a compartment upper-layer with a GER less than 1.5 was performed by 
Gottuk(1992a) using CHEMKIN and SENKIN. Gottuk modeled the compartment upper- 
layer as a constant pressure, constant temperature homogeneous reactor. The reactor 
temperature was equal to his experimentally measured compartment fire upper-layer 
temperatures, and the initial species levels were taken from the experiments of Beyler 
(1983). The model accurately predicted the species levels for overventilated 
compartment fires (0<1.1), but over-predicted CO levels by as much as 400% for 
underventilated compartment fires (6>1.1). In a later publication, Gottuk et al. (1995) 
attributed the discrepancy in the CO levels in underventilated compartment fires to the 
oxidation which occurs in the plume. | Better agreement with experimental results was 
achieved when the initial species compositions were adjusted to account for more 
complete oxidation in the plume. This was done by oxidizing all the unburned 
hydrocarbons (in this case C,H,) to CO, and H,O. The initial O, concentration was also 
reduced to 1%. The results of both modeling cases and of the experimental results of 
Gottuk (1992a) and Beyler (1986b) are shown in Fig. 6.1. 


The production of CO in compartment fires is now believed to occur within the fire 
plume which rises to the ceiling (Pitts, 1995, Gottuk, et al., 1995). The vortical structures 
of the fire are hypothesized to entrain pockets of air. This entrained air is subsequently 
carried into the upper-layer of the compartment fire to burn in an extremely 
underventilated environment. Modeling the fluid mechanics and chemical kinetics 


occurring in this process is extremely complex and, hence, is not attempted here. A 
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The CO experimental results of O-Beyler and A-Gottuk compared with the 
modeling results of Gottuk. Using O-Beyler’s data as the initial 
concentrations in the model, the dotted line is the result using the upper-layer 
temperature of Beyler and the A-dashed line is the result using the upper- 
layer temperature of Gottuk. The W-solid line corresponds to CO levels 
where all of the UHC was assumed to be consumed in the plume. 
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simple method which provides an estimate of the upper-layer species concentrations 


expected to be present in a compartment fire during post-flashover is developed. 


The basis of the method is to use the conditions in the compartment just prior to 
flashover as initial conditions in a homogeneous reactor to estimate the species 
concentrations expected after flashover. The GER of the compartment fire and upper- 
layer temperature are the two input parameters for the model. Hexane compartment fire 
modeling results were compared with the experimental measurements of Gottuk (1992a). 
Modeling results of methane fires were compared with the species concentration 
measurements in natural gas compartment fires conducted by Bryner et al. (1994, 1995) 


The effects of temperature on species levels is also investigated. 
6.3.2 Description of the Compartment Fire Model 


Two approaches were considered when modeling the combustion products in the 
upper-layer of the compartment. The compartment upper-layer temperature and GER are 


the necessary inputs in both approaches. 


Equilibrium calculations performed using STANJAN were determined to be adequate 
to model overventilated (GER<0.9) compartment fires. STANJAN was run with the 
constant pressure and constant enthalpy conditions invoked (in other words a constant 


pressure, adiabatic system). 


A homogeneous reactor was used to estimate species concentrations in an 
underventilated (GER>0.9) compartment fire upper-layer. The assignment of the initial 
species mass fractions in the homogeneous reactor was critical. The initial species levels 
were determined using a model of the fire structure inside the compartment present just 


prior to the onset of flashover, see Fig. 6.2. 


A compartment fire, as seen in Fig. 6.2, was considered to consist of two distinct 
parts: a fire plume and an accumulated upper-layer. The plume was assumed to be 
divided into two zones, an external flame sheet and a pure fuel core. Since the plume was 
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Figure 6.2 A phenomenological representation of the conditions present inside an 
idealized compartment fire just prior to the onset of flashover. 
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considered to rise to the upper-layer in an overventilated environment, the flame sheet 
was considered to have an equilibrium composition with an initial equivalence ratio of 
,=0.9. The remainder of the fuel in the compartment was assumed to be convected into 


the upper-layer through the core. 


Experimental data was utilized to develop a basis for the general estimation of 
species concentrations accumulated in the compartment upper-layer before flashover. 
The O, and CO, concentrations present in the compartment upper-layer just ptior to the 
onset of flashover were measured in n-hexane fires to be approximately 4 % and 10%, 
respectively. These levels are consistent with O, and CO, concentrations measured in the 
upper-layer by Gottuk (1992a) at a GER of approximately 0.8. Hence, the upper-layer 
prior to flashover was modeled as an equilibrium mixture with an initial equivalence ratio 


of by, =0.8. 


To predict the species concentrations in the upper-layer during the steady-state 
period (i.e. post-flashover) of the compartment fire, the upper-layer was modeled as a 
homogeneous isothermal reactor. The initial composition in the reactor was obtained by 
adding the moles of products from the flame sheet, the moles of fuel in the plume core 


and the moles of species accumulated in the upper-layer prior to flashover. 


The detailed steps for modeling underventilated compartment fire upper-layers are 
listed below. 
1. Simulate upper-layer accumulated during pre-flashover by running STANJAN at 


,,=0.8 with the constant pressure and enthalpy conditions invoked. 


2. Simulate the burning portion of plume (the flame sheet) rising into upper-layer by 
running STANJAN at o,=0.9 with the constant pressure and the constant enthalpy 


conditions invoked. 


224 


Chapter 6. Engineering Tools for the Estimation of CO Levels in Building Fires 


3. Subtract the moles of fuel burned in step 2 from the global equivalence ratio of the 
upper-layer during the steady-state, and add this amount of fuel (in moles) to the 


result of step 2. This simulates the species in the plume entering the upper-layer. 


4. Add the result of step 1 and step 3 and use this mixture as the initial condition in 
SENKIN with the constant pressure and constant temperature (steady-state gas 


temperature) conditions invoked. 


The code was run until the species concentrations came to a steady state. This was 
assumed to occur when a less than 0.05% change in the CO level was observed, 


approximately 15 seconds of computer run time. 


The simulations considered both highly overventilated (GER=0.1) and 
underventilated (GER=4.0) compartment fires. The initial molar concentrations in the 
homogeneous reactor was determined using the procedure previously described. The 
SENKIN calculations were performed at various upper-layer temperatures to include the 
conditions of all available experimental data. In particular, temperatures of 900 , 1000 
and 1100 K were used in the computations to model the reduced-scale n-hexane 
compartment fires. Temperatures of 1000 K and 1300 K were used to model the reduced- 


scale methane fires. 
6.3.3. n-Hexane Compartment Fire Modeling Results 


In the experiments performed by Gottuk (1992a), the average upper-layer 
temperature inside the compartment during steady-state was approximately 1000 K with a 
variation of +100 K. Measurements of CO, CO, and O, concentrations in the 


compartment upper-layer were made for a GER ranging from 0.3 to 3.0. 


The CO concentration data of Gottuk (1992a) is shown in Fig. 6.3 plotted against the 
GER. Also shown in Fig. 6.3 are the predictions of the current model and the modeling 


results of Gottuk (1992a). Using a 1000 K_  as_ the  upper-layer 
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Comparison of the O-experimental and —modeling results of Gottuk(1992a) 
with the proposed model estimations of CO concentrations for n-hexane fires 
during the steady-state time with temperatures of ----900K, «+*+ 1000 K, and 
—°—1100K. 
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temperature, the proposed model estimated the data within +10% for GER greater than 
1.8 and less than 1.4. With the GER between 1.4 and 1.8, the model over-predicted the 
experimental results by as much as 50 %. Though this error seems high, the estimations 
of CO concentration using the model of Gottuk is seen in Fig. 6.3 to be in error by as 


much as 700% in this range. 


The model estimations of the CO, and O, concentrations are plotted in Fig. 6.4 and 
Fig. 6.5, respectively. The model estimates the experimental data to within +10% over 


the entire range of GER. 


The model estimations of other species concentrations formed in a compartment fire 
at 1000 K are plotted in Fig. 6.6 versus the global equivalence ratio. The ignition index 
calculations assumed the H, concentration was half the concentration of CO. The 
modeling results verify this assumption to be accurate to within +20% over all GER. The 
concentration of different UHC in the upper-layer is plotted in Fig. 6.7. The UHC with 
the highest concentration was ethylene. This verifies the assumption made when using 


ethylene as a calibration gas for the FID analyzer. 
6.3.4 Methane Compartment Fire Modeling Results 


Species concentration measurements in natural gas compartment fire experiments 
were performed by Bryner et al. (1994). Measurements were performed in both the rear of 
the compartment (near the wall on the opposite end of the room from the doorway) and in 
the front of the compartment (by the wall containing the doorway). Temperatures of 
approximately 1000 K were measured in the rear while temperatures of 1300 K were 
measured in the front. The CO concentrations estimated using the model at the 
temperatures in the rear (1000 K) are seen in Fig. 6.8 to be within +25% of the 
compartment rear experimental data. The CO concentrations estimated for the front of 
the compartment (1300 K) are seen in Fig. 6.8 within +30%. The model consistently 
under-predicted the experimental data from a GER of 0.6 to 1.2. 


227 


Chapter 6. Engineering Tools for the Estimation of CO Levels in Building Fires 


el ee 
Go fF WN 


CO, Concentration, [%-wet] 
N U2 - N Oo ~“ io) Ne) 


pe 


f=) 


0.0% 0.5, 91.0 leo 22.0 2S 3 03 Samed. Oa 
Global Equivalence Ratio, o 


Figure 6.4 Comparison of the O-experimental and — modeling results of Gottuk (1992a) 
with the proposed model estimations of CO, concentrations for n-hexane fires 
during the steady-state time with temperatures of ----900K, +++ 1000 K, 
and —-— 1100 K. 
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Figure 6.5 Comparison of the O-experimental and — modeling results of Gottuk (1992a) 
with the proposed model estimations of O, concentrations for n-hexane fires 
during the steady-state time with temperatures of ----900K, *-+* 1000 K, 
and —*-— 1100 K. 
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Figure 6.6 Model estimation of the major species concentrations present in the 1000 K 
compartment upper-layer during the steady-state time of a n-hexane fire. 
Symbols: O-N,, O-H,O, O)-CO,, A-O,, V-CO, 0-H, and @-UHC. 
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Figure 6.7 Model estimation of the UHC concentrations present in the 1000 K 
compartment upper-layer during the steady-state time of a n-hexane fire. 
Open symbols, O-C,H,, D1-CH,, A-C3H,, and V-C,H>, correspond to the left 
y-axis while closed symbols, @-C,H,, M-C,H, and A-C;H; correspond to the 
right y-axis. 
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Comparison of the experimental results of Bryner, et al. (1994), in both the 
O-front (1300 K) and the DJ-back (1000 K) of their compartment, with the 
proposed model estimations of CO concentrations for methane compartment 
fires during the steady-state time with temperatures of ----1000K and 
°° ¢1300K. 
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The estimations of CO, and O, concentrations in the upper-layer are plotted in 
Fig. 6.9 and 6.10, respectively. The model estimated the CO, and O, concentrations to 
within +10% of the experimental data over the entire range of GER study. 


The model estimations of other species concentrations in a compartment upper-layer 
at 1000 K are shown in Fig. 6.11. The H, concentrations were measured to be nearly 
equal to the concentration of CO. The difference in CO-H, ratio between the methane 
experiments and the n-hexane experiments was attributed to the carbon-hydrogen ratio in 
the fuel. Methane has a 1:4 carbon-hydrogen ratio while n-hexane has a 1:2.3 carbon- 
hydrogen ratio. The estimations of the concentrations of different unburned 
hydrocarbons in the upper-layer is plotted in Fig. 6.12. Methane is the dominant 
hydrocarbon in these compartment fires. 
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Figure 6.9 Comparison of the experimental results of Bryner, et al. (1994), in both the 
O-front (1300 K) and the O-back (1000 K) of their compartment, with the 
proposed model estimations of CO, concentrations for methane compartment 
fires during the steady-state time with temperatures of ----1000K and 
¢¢°1300K. 
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Figure 6.10 Comparison of the experimental results of Bryner, et al. (1994), in both the 
O-front (1300 K) and the D-back (1000 K) of their compartment, with the 
proposed model estimations of O, concentrations for methane compartment 


fires during the steady-state time with temperatures of ----1000K and 
e¢¢ 1300 K. 
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Figure 6.11 Model estimation of the major species concentrations present in the 1000 K 
compartment upper-layer during the steady-state time of a methane fire. 
Symbols: O-N,, O-H,O, LI-CO,, A-O,, V-CO, 0-H, and @-UHC. 
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Figure 6.12 Model estimation of the UHC concentrations present in the 1000 K 
compartment upper-layer during the steady-state time of a methane fire. The 
open symbol, O-CH,, corresponds to the left y-axis while closed symbols, 
W-C,H,, A-C,H,, and W-C,H, correspond to the right y-axis. 


Chapter 7. Conclusions and Recommendations 


CHAPTER 7 


CONCLUSIONS AND RECOMMENDATIONS 


7.1 SUMMARY 


Smoke inhalation accounted for 76% of the deaths in building fires during 1990. 
Two-thirds of these victims were found at locations remote from the burning room. 
Carbon monoxide (CO) has been determined by numerous studies to be the major 
toxicant in the fire gases, and has been deemed responsible for the majority of the smoke 


inhalation deaths in building fires. 


The level of CO transported to remote locations from a burning compartment is a 


function of the 
e global equivalence ratio in the compartment, 
e oxygen entrainment into gases burning in the hallway, and 
e heat transfer. 


The entrainment into the gases entering the hallway strongly affects the last two 
variables. Therefore, it was hypothesized that this entrainment strongly controls the CO 
level transported to remote areas. This experimental study focused on the effects of 
entrainment into the compartment fire gases entering the hallway on the CO levels at 


locations remote from the burning compartment. 


The results of the experimental investigation were discussed in Chapters 3-5. The 
level of CO exiting the hallway (termed post-hallway yields) was measured for different 
experimental conditions. The transport of CO away from a burning compartment on the 


side of a hallway was characterized by sampling in the hallway. Experimental data was 
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used to identify the key parameters controlling the transport of high CO concentrations. 
Correlations were developed using the experimental data to estimate CO levels 
transported to remote locations. Models which identify when these correlations can 
accurately be applied to building fires were also developed. A general procedure was 
proposed for the estimation of CO levels transported in building fires. A chemical kinetic 
model was also developed to estimate the species concentrations present in a 


compartment fire upper-layer after the onset of flashover. 


Through the research, the following parameters were identified to most significantly 


affect the level of CO transported to remote locations. 

1. entrainment rate of surrounding gases into compartment fire gases, 
2. oxygen level in the surrounding gases being entrained, 

3. presence or absence of external burning in the hallway, 

4. stoichiometry of the gases entering the hallway, and 

5. the transient time of the species concentrations in the hallway. 


Overventilated n-hexane compartment fire experiments were performed to measure 
the entrainment rate under different experimental conditions. The entrainment rate was 
most drastically affected by the compartment-hallway opening size and the soffit height 
above the opening. A reduction in the opening size resulted in an increase in the 
entrainment. The entrainment was measured to exponentially increase with a decrease in 
the Richardson number of the flow entering the hallway. An increase in the entrainment 
rate was also measured when the soffit height was increased. This was attributed to the 
presence of the soffit allowing the gases to enter the hallway as a buoyant plume instead 


of a ceiling jet. 


During a building fire, a deep upper-layer of oxygen deficient combustion gases 


develops in the space adjacent to the burning room. The depth of this upper-layer was 
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hypothesized to significantly affect the O, level entrained into compartment fire gases 
entering the space. Underventilated n-hexane experiments were conducted to study the 
effect of gradually increasing the depth of the upper-layer in the hallway on the post- 
hallway species yields. The degree to which the hallway upper-layer surrounds the jet of 
compartment gases entering the hallway was defined by the nondimensional upper-layer 
depth (y=8/z, where 4 is the upper-layer depth and z is the distance between the bottom of 
the opening and the ceiling). The post-hallway CO and UHC yields were measured to 
increase with an increase in y (by increasing the upper-layer depth). This was attributed 
to less oxygen being entrained into the compartment fire gases entering the hallway. The 
study of Westbrook and Dryer (1984) noted that a decrease in O, levels at the flame front 
results in limited radicals being formed to oxidize CO and UHC. The CO yield was 
actually measured to increase to levels 6-25% greater than measured compartment yields. 
This was measured in experiments with a 0.12 and 0.08 m? opening sizes, y between 1.4 
and 1.7, and external burning occurring in the hallway. The jets associated with these 
opening sizes were measured to have the lowest entrainment rates of experimentally 
tested jets. From the results of Westbrook and Dryer (1984), the low entrainment rate 
coupled with the low oxygen level being entrained is believed to have resulted in a 
limited pool of radicals. It is hypothesized that the majority of these radicals were used to 
oxidize UHC thus forming additional CO in the process (Westbrook and Dryer, 1984). 
The increase in CO yields is attributed to limited radicals being able to react with the CO 
before the oxidation reactions were frozen (Westbrook and Dryer, 1984). With y>1.7, no 
external burning was observed to occur in the hallway. Post-hallway CO yields were 


measured to be equal to those measured inside the compartment. 


No external burning was observed to occur when polyurethane foam slabs were 
burned in the compartment. The yields downstream were, therefore, always nearly equal 
to the species yields measured inside the compartment (0.19 for CO) regardless of the 


hallway upper-layer depth. 
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With the compartment on the side of the hallway, the fluid mechanics in the hallway 
are highly 3-dimensional. Flow visualization of the hallway fluid mechanics was 
performed by observing the external burning occurring in the hallway. The burning was 
observed at the open end of the hallway. The gases were observed to flow initially across 
the hallway and were then seen to impinge on the wall opposite the compartment. The 
impingement resulted in large clockwise rotating vortical structure. The vortical structure 
swirled down the hallway along the corner of the ceiling and the wall opposite the 
compartment. From the size of the structure, the majority of the gases flowing away from 
the room appeared to be transported down the hallway within the structure. 


The temporal and spatial variation in species concentrations and temperature in the 
hallway was investigated for experiments where high post-hallway CO yields were 
measured. Concentrations of CO exceeding 1.8% were measured in gases flowing across 
the hallway and down the side of the hallway opposite the compartment shortly after 
compartment flashover. This is precisely the flow pattern observed from the flow 
visualization. The concentration of CO measured in the gases on the compartment side of 
the hallway at this period of the fire was less than 0.8%. The non-uniformity of CO 
concentration in the hallway persisted until the hallway species concentrations were 
measured to reach a steady-state in the hallway. This was approximately 60 seconds after 
flashover in experiments with a 0.60 m exit soffit and 96 seconds in experiments with the 
bottom of the hallway exit blocked. The CO concentrations during the steady-state 
period were measured to be approximately 2.0% within the upper-layer despite the 
occurrence of external burning. The burning was measured to reduce UHC levels to 
approximately 1.0%. Compartment fire exhaust gas temperatures cooled to less than 
800 K approximately 4 seconds after entering the hallway with temperatures never 
measured in excess of 1000 K. The CO oxidation was, therefore, believed to occur at or 


near the flame front where temperatures greater than 1200 K exist. 


241 


Chapter 7. Conclusions and Recommendations 


The effects of fuel type in the compartment was studied by conducting a series of 11 
polyurethane foam fire experiments. The variation in the species concentrations and 
temperatures were measured in the hallway 0.05 m below the ceiling. No external 
burning occurred in these experiments. The flow of gases containing high concentrations 
of CO was measured to be nearly equivalent to that measured in the n-hexane fire 
experiments. The species concentrations were measured to equilibrated approximately 60 
seconds after compartment flashover. Approximately 1.8% CO was measured to be 


uniform within the upper-layer gases 0.05 m below the ceiling. 


External burning in the hallway was measured to have significant effect on the post- 
hallway yields. A model to determine the occurrence of external burning was developed. 
The model uses the ignition index and the nondimensional upper-layer depth to predict 
external burning. The ignition index was used to determine whether the mixture entering 
the hallway was ignitable. The nondimensional upper-layer depth determined whether a 
pilot ignition source (in the form of oxidative radicals) was available. Ignition was 
determined to occur in the hallway when the ignition index exceeded approximately 1.0 


and the nondimensional upper-layer depth was less than 1.7. 


Correlations to predict CO levels transported to remote locations are most 
conveniently expressed in terms of yields. For the post-hallway yields to accurately 
predict the level of species in the hallway, the concentration of the species must be 
uniform within the hallway. A perfectly stirred reactor model was developed to estimate 
the time necessary for the species concentrations to reach a steady-state (termed the 
transient time). The transient time from the model was determined to be within +15% of 


the experimentally determined transient times. 


Post-hallway yields were measured to be equal to compartment yields when no 
external burning occurred in the hallway. In this case, correlations based on the 


compartment global equivalence ratio can be used to predict post-hallway CO yields. 
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The occurrence of external burning in the hallway led to the oxidation of species in the 
hallway. Post-hallway yields were determined to no longer be correlated with the 


compartment global equivalence ratio when external burning occurred . 


Two models are proposed to predict post-hallway CO yields with external burning 
occurring in the hallway. The first model relates the oxygen entrained into the 
compartment fire gases entering the hallway to the post-hallway yields of CO and UHC. 
The model accounts for the entrainment rate of the jet and the concentration of oxygen 
entrained by the oxygen entrainment parameter. The model is accurate for a 
compartment global equivalence ratio greater than 1.8 and for hallway upper-layer 
depths less than y=1.7. A second model relates the post-hallway yields to a global 
equivalence ratio for the hallway and compartment. The oxidation of CO in the hallway 
is most significantly affected by the air entrained into the burning gases in the hallway. 
This entrainment was approximated using the air entrainment measurements and a 
estimation of the flame extension in the hallway. The data was determined to collapse to 
the compartment data of Gottuk (1992). This was expected since the system global 
equivalence ratio in the study performed by Gottuk was the compartment global 
equivalence ratio. The model is accurate for all compartment global equivalence ratios, 


but is only valid for upper-layers with y<1.4. 


A phenomenological model which considers only chemical kinetics was developed 
to provide an estimation of species concentrations inside the compartment during post- 
flashover. Previous models have relied on available experimental data to set the initial 
conditions in a reactor. The proposed model requires knowledge of the compartment 
global equivalence ratio and upper-layer temperature to provide estimations of species 
concentrations in the compartment upper-layer. The model was able to estimate CO, 
CO,, and O, concentrations to within +25% for both n-hexane and methane compartment 


fire data. The model was also able to estimate higher levels of CO measured in the full- 
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scale methane fire experiments conducted by Bryner, et al.(1995). With validation of the 
models agreement with available experimental data, the model verified several 
approximations and assumptions which have been used in the study. The model confirms 
H, levels being half of the CO level in the 900-1100K temperature range, in addition to 
confirming that ethylene, C,H,, is the major unburned hydrocarbon produced in the n- 


hexane fires. 


7.2 CONCLUSIONS 


The goal of the experimental study was to provide fire protection engineers with 
correlations to predict CO levels transported to remote locations. The experiments were 
conducted with the compartment on the side of the hallway. The research identified the 
parameters most significantly affecting the post-hallway CO yields. The temporal and 
spatial variation of the CO concentration in the hallway upper-layer was measured with 
the burning compartment on the side of the hallway. A procedure was developed to 
provide engineers with a sound method for predicting CO yields in building fires. An 
engineering tool was also developed to predict species formation in a compartment fire 


after flashover. 
The experimental investigation has determined that 


e the existence of a deep oxygen deficient upper-layer in the adjacent space 
contributes significantly to the transport of fatally high levels of CO in building 
fires, and 


e with the burning room on the side of the hallway, fire gases are transported non- 


uniformly down the hallway. 


Both of these conditions were present in the unfortunate fire which took place in the 
Hillhaven Nursing Home in Norfolk, Virginia during October of 1989. In this fire, 13 
people died with 10 people perishing on the wing of the floor containing the room with 
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the fire. Nine of the fatalities were located in rooms across the hallway from the burning 
room while only one occurred on the side of the hallway containing the burning room. 
The gas distributions in the hallway adjacent to the burning room are believed to be 
similar to the experiments with a long transient time. Disturbingly, incidents such as 


these are not isolated and occur several times annually. 


Present fire protection engineering computer models would not have predicted such a 
hazard to the Hillhaven patients on the side of the hallway opposite the burning room 
since these codes use zonal models. Zonal models assume a spatially uniform CO yield 
in the upper-layer along the entire length of the hallway. With the burning room located 
on the side of the hallway, results show the transport of gases containing high levels of 
CO are transported to remote locations is highly non-uniform. The duration of time 
which the non-uniformity in the hallway exists depends on the characteristics of the 
adjacent space (i.e. hallway length, number of rooms connected to the hallway). Field 
models allow spatial and temporal variation in the species within the space applied. 
These models can be used to accurately predict the toxic gas flow away from a room on 
the side of a hallway at all times during the fire. However, field models are expensive 


and time consuming. 


Some situations may not require field models to provide good estimates of the toxic 
gas flow. The species concentrations in the adjacent space were measured in this study to 
be uniform after the species concentrations had reached a steady-state. The transient time 
can be estimated using the perfectly stirred volume model. With species concentrations 
uniform in the adjacent space hallway, zonal models may be used to predict the toxic gas 


hazard in the space. 


An engineering tool was also developed to provide a more detailed estimate of the 


concentrations present in a compartment fire upper-layer after post-flashover. The model 


245 


Chapter 7. Conclusions and Recommendations 


requires the upper-layer temperature and the global equivalence ratio of the compartment 


as inputs. 


7.3 RECOMMENDATIONS 


The research included in this thesis has addressed the effects of oxygen entrainment 
on the transport of CO to remote locations when the burning compartment is on the side 
of the hallway. Additional experiments need to be performed to measure the effects of 
heat loss, a door connecting the hallway and compartment, and the compartment-hallway 
orientation on the CO yields and in-hallway concentrations. The effects of these 
parameters on the CO yields need to be compensated for in the correlations developed in 
this work. An engineering tool needs to be developed to estimate species levels in the 


hallway. 


Heat losses in the hallway are responsible for the gases cooling to a point where 
oxidative reactions no longer occur. To gain a full understanding of the phenomena 
controlling the transport of high levels of CO, measurements of heat losses in the hallway 
need to be performed. Specifically, the heat transfer losses to the surroundings (i.e. the 
walls, the area below the upper-layer) and the heat transfer effects of the air entrainment 


need to be determined. 


Experiments measuring the entrainment into the burning hallway gases need to be 
performed. Besides the measurements made in this study, very little data is available on 
the entrainment of air into ceiling jet flows. Measurements made in this study provided a 


good approximation for these levels, but improved measurements are suggested. 


Experiments need to be conducted with the compartment shifted from the end of the 
hallway to the middle of the hallway forming a T shape. In this scenario, the gases would 


split off into two directions after impinging on the wall opposite the compartment. In one 
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direction, the gases would flow toward the open end of the hallway where the fume hood 
is located, possibly leaving the hallway. The gases moving in the other direction are 
forced to circulate within the dead end of the hallway. The dead end of the hallway has 
been speculated by many to be an area where high levels of CO are produced because of 
the long residence times of the gases in this area (Nelson, 1996). The temporal and 
spatial evolution of hallway CO concentrations in this compartment-hallway scenario 


would be particularly interesting. 


The burning room is typically connected to the adjacent space by a door . From the 
results in Appendix C, a doorway connecting the compartment and the hallway has an 
effect on both the level of species produced inside the compartment and the oxidation of 
gases within the hallway. Further experiments should be conducted with a door (where a 
bi-directional flow of cool air entering the room and hot combustion gases exiting the 
room is present) to quantify the effects of deep hallway upper-layers on CO levels. The 
application of the nondimensional upper-layer depth for doorway scenarios should be 
tested. The distance between the neutral plane of the doorway and the ceiling (instead of 
the distance between the bottom of the opening and the ceiling) should be investigated. 
The neutral plane is typically 1.44 m below a 2.44 m high ceiling. External burning will 
most likely occur in all experiments since a maximum value of y is slightly less than 1.7. 
External burning was observed not to occur when y>1.7 in experiments with a uni- 


directional flow opening. 


A simple plug flow reactor model is desired for the estimation of species oxidation in 
the hallway. With no plug flow reactor code available, assumptions were applied to the 
governing equations in Appendix D to allow the momentum equation to be disregarded. 
This leaves just the species conservation equations to be solved; a task which can be 
accomplished using a code such as SENKIN. Since SENKIN solves such equations with 


time and all the data is in spatial coordinates, the transformations necessary to convert the 
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spatial results to a temporal coordinate system, and vice versa, are also developed in 
Appendix D. The plug flow reactor model developed in Appendix D. would use the 
results of the compartment fire model as the species mole fractions initially entering the 
reactor. Experimentally determined temperatures and air entrainment into the upper-layer 


gases would be used to apply the model. 
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APPENDIX A. 


EFFECT OF ORIFICE DIAMETER ON AIR INLET DUCT 
VELOCITY PROFILES 


Different size orifices were attached to the end of the air inlet duct to vary the air 
entrainment rate into the compartment. The velocity profiles in the duct were measured 
using a hot film probe for the four different orifice diameters used in the study: no orifice 
(0.30 m), 0.25 m, 0.20 m and 0.15 m. The velocity profiles were used to determine the 


radial location of the probe to measure the average velocity in the duct. 


The velocity profile was measured by moving the hot film probe along the height of 
the air inlet duct in 0.025 m increments. Data was recorded for 150 seconds. Each data 


point is the result of averaging over a 90 second time window. 


Shown in Fig. A.1 are the velocity profiles generated with each air inlet condition 
listed above. As can be seen in Fig. A.1, the shape of the velocity profiles are very 
similar for all of the orifice diameters used in the study. The position of the average 
velocity of the profiles, therefore, does not change significantly, and was set at 0.044 m 


(1.75 inches) below the top of the duct. 
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Distance from Top of Duct, x/L 


0.0 0.2 0.4 0.6 0.8 1.0 
Velocity, v/V___.. 


Figure A.1 The velocity profiles in the air inlet duct for with @-no orifice (0.30 m), i- 
0.25 m, A-0.20 m and W-0.15 m diameter orifices attached to the duct. The 
vertical lines represent average velocities of the profiles. 
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APPENDIX B. 


DATA REDUCTION PROGRAMS 


FIRERED3.FOR 


C MODIFIED BY: BRIAN LATTIMER LAST UPDATE: 8/25/94 

C PREVIOUS VERSION: FIRERED2.FOR LAST VERSION BY: DAVID EWENS, 11/6/93 

C PREVIOUS VERSION: FIREDATA.FOR LAST VERSION BY: DAN GOTTUK, 7/31/92 

C 

C PURPOSE: 

C THIS PROGRAM REDUCES THE EXPERIMENTAL DATA FROM THE RAW DATA FILE 

C CREATED BY PROGRAM FIRES.BAS. 

C THIS PROGRAM COMPENSATES FOR ANALYZER AND FLOW LAG TIME (FID,CO,CO2,02). 
C DATA FILE CAN (SHOULD) INCLUDE DATA PAST END OF FIRE, AT LEAST 15 SECS. 

C THE REDUCED DATA MAY BE AVERAGED USING FIREAVG3.FOR 


c 


C VARIABLES: 


Scone hevet one ee ie "ele ie vel eketekekelekel eye a° lane ie 


AIRRAT 


AIRVEL(1) 


AREA 


CMOLIN 
CMOLOT 
COAMB 
COCST 
CORNG 
COYLD 
CO2AMB 


DENSTY 


DRYCO() 
DRYCO2(1) 
DRYO2(I) 


DCTFLW 


- INLET DUCT AIR MASS FLOW RATE [Kg AIR/S] 
- INLET DUCT AIR VELOCITY [M/S] 
- AREA OF AIR INLET DUCT [M’2] 
- CARBON COUNT, # C ATOMS IN FUEL MOLECULE [#] 
- RELATIVE CARBON ERROR [%] (Cout-Cin)/Cin 
- # OF CARBON MOLES IN FROM FUEL [KMOLES] 
- # OF CARBON MOLES OUT FROM CO, CO2 AND THC [KMOLES] 
- AVG AMBIENT ATMOSPHERE CO MOLE FRACTION [KMOL CO/KMOL AIR] 
- TO ADJUST CO [MOLE FRACTION/CONC.(ppm OR %)] 
- RANGE OF CO ANALYZER [ppm OR %] 
- CO YIELD [KG CO PRODUCED/KG FUEL BURNED] 
- AVG AMBIENT ATMOSPHERE CO2 MOLE FRACTION [KMOL CO2/KMOL AIR] 
- TO ADJUST CO2 [MOLE FRACTION/CONC.(ppm OR %)] 
- RANGE OF CO2 ANALYZER [ppm OR %] 
- CO2 YIELD [KG CO2 PRODUCED/KG FUEL BURNED] 
- DATE OF EXPERIMENT [TEXT] 
- HEAT OF COMBUSTION OF FUEL [KJ/Kg] 
- INTERMEDIATE VALUE TO CALCULATE FUEL RATE [S] 
- INTERMEDIATE VALUE TO CALCULATE FUEL RATE [S] 
- 2nd DERIVATIVE OF FUEL W.R.T. TIME [KG/S”2] 
- DENOMINATOR OF DRY TO WET CONC CONVERSION [UNITLESS] 
- INLET DUCT AIR DENSITY [KG/M’3] 
- DRY CO CONCENTRATION [ppm OR %] 
- DRY CO2 CONCENTRATION [%] 
- DRY 02 CONCENTRATION [%] 
- EXHAUST DUCT VOLUME FLOW RATE [M’3/S] 
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C EQUIV - PLUME EQUIVALENCE RATIO [UNITLESS] 
C EXMOLS —_- TOTAL NUMBER OF MOLES IN COMPARTMENT [KMOL EXHAUST/S] 
ASSUMES EXHAUST HAS AVG MOLE WT OF AIR!!! 
EXCOEF  - EXTINCTION COEFFICIENT FROM BEER'S LAW [1/m] 
EXTDAT _- REDUCED DATA FILE EXTENSION [TEXT] (ie. ' DTA’) 
EXTRAW  - RAW DATA FILE EXTENSION [TEXT] (ie. ' RAW’) 
FOWT — - BEGINNING WEIGHT OF FUEL [Kg] 
FID(I) - SIGNAL FROM FID [mV*FIDRANGE] 
FIDCAL _- FID SPAN GAS READING [mV * FIDRANGE] 
FIDCST -TO ADJUST THC [MOLE FRACTION/ppm C2H4] 
FIDRNG  - FID RANGE [MULT FACTOR] 
FIDSPN _- FID SPAN GAS CONCENTRATION [ppm C2H4] 
FIDZER _- FID ZERO (ATMOSPHEREIC SAMPLE) READING [mV * FIDRANGE] 
FNAME _ - DATA FILE NAME [7 CHAR'S, ie. HALL001] 
FNAMEX - TOTAL NAME FOR OPENING ALL FILES [TEXT] (ie. 'HALL001 RAW’) 
FTOA — - MEASURED MASS FUEL TO MASS AIR RATIO [KG FUELKG AIR] 
FUEL(I) - FUEL MASS [KG] 
FUELRT  - FUEL MASS BURN RATE [KG FUEL/S] 
GF - FUEL VOLATILIZATION RATE [g FUEL/S] 
HTOC  -HTOCRATIO OF FUEL [# H## C, UNITLESS] 
I - TOTAL NUMBER OF RECORDEDD DATA POINTS [#] 
ICORG  -CORANGE (1,2, OR 3) [#] 
ICO2RG  -CO2 RANGE (1, 2, OR 3) [#] 
IFIDRG _- FID RANGE [POWER OF 100] 
IFUELT  - FUEL TYPE [#] 
= 1 HEXANE 
=2 PMMA 
=3 WOOD 
= 4 POLYURETHANE 
IPROBE - PROBE LOCATION [#] 
= 1 IF COMPARTMENT SAMPLED 
= 2 IF EXHAUST DUCT SAMPLED 
= 3 IF HALLWAY SAMPLED 
IRUNNO  - RUN NUMBER [#] 
ISTEP = - DATA POINT USED STEP SIZE [#] (1 = USE ALL) 
IQUIT - DUMMY ABORT REDUCTION FLAG VARIABLE [NONE] 
JJ. | - DUMMY COUNTER INTEGER [#] 
K - DATA POINT COUNTER [#] 
Kl - SPECIES ANALYZER LAG COMPENSATED TIME DATA POINT COUNTER [#] 
K2 - FID LAG COMPENSATED TIME DATA POINT COUNTER [#] 
LASPTH -EXHAUST DUCT SMOKE LASER PATH LENGTH [m] 
MASODL _ - MASS OPTICAL DENSITY OF SMOKE PER UNIT PATH LENGTH [m“2/g] 
C MOLWT  - MOLECULAR WEIGHT OF FUEL [KG FUEL/KMOLE FUEL] 
CN - DUMMY INTEGER COUNTER [#] 
C NOTE] | - ANY NOTE ON EXPERIMENT RUN [TEXT] (60 CHAR'S MAX) 
C NOTE2 -ANY NOTE ON EXPERIMENT RUN [TEXT] (60 CHAR'S MAX) 
C NOTE3  -ANY NOTE ON EXPERIMENT RUN [TEXT] (60 CHAR'S MAX) 
C NOTE4 - ANY NOTE ON EXPERIMENT RUN [TEXT] (60 CHAR'S MAX) 
C ODL - OPTICAL DENSITY OF SMOKE PER UNIT PATH LENGTH (1/m) 


‘Mele eNOlelete @tekohe ene he veleroteleve te kelren¢ else). erelenshelekecelomel eile). 
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C O2AMB_—_- AVG AMBIENT ATMOSPHERE 02 MASS FRACTION [Kg 02/Kg AIR] 
C O2CST  -TO ADJUST 02 [MOLE FRACTION/CONC.(ppm OR %)] 

C O2RNG —_ - MAXIMUM MEASURABLE 02 CONCENTRATION [%] 
C O2YLD  - O02 YIELD [KG 02 USED/KG FUEL BURNED] 

C PA - AMBIENT PRESSURE [mm Hg] (TYP. 710 mmHg) 

C PBELOC  - DESCRIPTION OF PROBE LOCATION IN HALLWAY [TEXT] 
C PRESS _- EXHAUST DUCT ORIFACE PRESSURE DROP [mmHg=Torr] (TYP. 15.8) 
C PROBE - PROBE LOCATED DESCRIPTION [TEXT] 

C PNSIZE  - FUEL PAN SIZE [inches] 

C PNWT _ - WEIGHT OF FUEL PAN [Kg] 

C PTSPST -#OF DATA POINTS PAST END OF FIRE [#] 

C PTSTRT  -LAG TIME AS# OF DATA POINTS [#] (FUNC. OF PROBE LOCATION) 
CQ - HEAT RELEASE RATE [KW] 

C RESTIM _ - RESIDENCE TIME OF GASSES IN COMPARTMENT [S] 
C RUNTIM- - TOTAL RUNTIME [S] 

C SEC - SPECIFIC EXTINCTION COEFFICIENT FOR SMOKE [m’2/g] 
C SMK10 - SMOKE0 CALCULATED FROM IST 5 SMOKE] DATA POINTS 
e [RELATIVE LIGHT INTENSITY] 

C SMKMAX _- MAX SMOKE VOLUME (MOLE) FRACTION DURING FIRE [ppb] 
C SMKYLD- - SMOKE YIELD [KG SMOKE PRODUCED/KG FUEL BURNED] 
C SMOKEO(I) - MEASURED I0 FOR BEER'S LAW (NOT WORKING/NOT USED) 
G [RELATIVE LIGHT INTENSITY] 

C SMOKE1(I) -IFOR BEER'S LAW [RELATIVE LIGHT INTENSITY] 
C SMKVF -_- SMOKE VOLUME FRACTION IN DUCT [ppm] 

C SOFF1  - HEIGHT OF SOFFIT FROM COMPARTMENT TO HALLWAY [inches] 

C ALSO CEILING HEIGHT ABOVE TOP OF COMPARTMENT VENT 

C SOFF2. - HEIGHT OF SOFFIT AT EXIT OF HALLWAY [inches] 

C SSTIME- - RATIO OF RESIDENCE TIME TO D2FUEL [UNITLESS] 

C STOIC  - STOICIOMETRIC AIR TO FUEL RATIO BY MASS [KG AIR/KG FUEL] 

C TAMBA — - AMBIENT ATMOSPHERIC TEMPERATURE INPUT BY USER [K] 

C TAMBC _- AMBIENT COMPARTMENT TEMPERATURE INPUT BY USER [K] 

G AVERAGE THE TOP 3 COMPARTMENT TC RAKE THERMOCOUPLES 

C TCC(1) | - COMPARTMENT THERMOCOUPLE TREE TEMPERATURES [K] (1 = TOP) 
C 4TCC(s) 

C TCC3. -3RD TC OF COMPARTMENT TC TREE TEMPERATURE [K] 

C TCH(1) -HALLWAY THERMOCOUPLE TREE TEMPERATURES [K] (1 = TOP) 

C - TCH(9) 

C TCEV  - EXHAUST VENT THERMOCOUPLE [K] 

C TCED  -EXHAUST DUCT THERMOCOUPLE [K] 

C TCID - AIRINLET DUCT THERMOCOUPLE [K] 

C TDMOLS'~— - TOTAL NUMBER OF MOLES IN EXHAUST DUCT [KMOL EXHAUST/S] 
C TFUEL  _- FUEL TYPE NAME [TEXT] 

C THC - UNBURNED HYDROCARBONS [KMOLES UHC/KMOLES TOTAL] 

C TIME - TIME FROM START OF FIRE [S] 

C TMEPST ~ - TIME PAST END OF FIRE THAT DATA WAS COLLECTED [S] 

C VENT _ - VENT SIZE [TEXT: inches] 

C VOLC  -AVG VOLUME OF UPPER LAYER [M’3] 

C WETCO = - WET CO CONCENTRATION [ppm OR %] 

C WETCO2 = - WET CO2 CONCENTRATION [%] 
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C WETO2 - WET 02 CONCENTRATION [%] 

Cc 

C 

C FILES ACCESSED: 

C ( HALL####.*, WHERE #### IS THE RUN NUMBER, ie. HALL0001 ) 
C 

C READ: 

C 

C HALL####.RAW - RAW EXPERIMENTAL SERIAL DATA FILE CREATED BY FIRE3.BAS 
C 

C WRITTEN: 


HALL####.DTI - EXPERIMENTAL INFORMATION DATA FILE (CONTAINS: FILE NAME 
TEST DATE, AMBIENT CONDITIONS, ANALYZER SETTINGS, ETC) 
HALL###.DTA - EXPERIMENTAL CHRONOLOGICAL DATA FILE (CONTAINS: TIME, 
WET CONC'S, YIELDS) 

HALL###.DTB - EXPERIMENTAL CHRONOLOGICAL DATA FILE (CONTAINS: TIME, 
FUEL & INLET AIR DATA, CARBON ERROR, EQUIVALENCE RATIO) 

HALL####.DTC - EXPERIMENTAL CHRONOLOGICAL DATA FILE (CONTAINS: TIME, 
TIME SCALES, DRY CONC'S, THC, HEAT RELEASE, DILUTION) 

HALL####.DTD - EXPERIMENTAL CHRONOLOGICAL DATA FILE (CONTAINS: TIME, 
SMOKE DATA, WET H20 & H2 CONC'S) 

HALL####.DTE - EXPERIMENTAL CHRONOLOGICAL DATA FILE (CONTAINS: TIME, 
HALLWAY THERMOCOUPLE DATA) 

HALL###.DTF - EXPERIMENTAL CHRONOLOGICAL DATA FILE (CONTAINS: TIME, 
COMPARTMENT THERMOCOUPLE DATA, INLET AIR & EXHAUST 
THERMOCOUPLE DATA) 


Co CC) COG) OCC) 0 Cit) CsG) OO) &) 


C2345-7-9012345678901234567890123456789012345678901234567890123456789 72 
C * INITIALIZATION 
IMPLICIT DOUBLE PRECISION (A-H,O-Z) 
DIMENSION TIME(1300), DRYCO(1300),DRYCO2(1300),DRY02(1300) 
DIMENSION FUEL(1300), AIRVEL(1300), TCC(8), TCH(9),TCC3(1300) 
DIMENSION FUELRT(1300),AIRRAT(1300),RESTIM(1300),Q(1300),Q2(1300) 
DIMENSION D2FUEL(1300),SMOKE1(1300),FULRT2(1300), FUEL2(1300) 
DIMENSION FID(1300),FTOA(1300), TCEV(1300), TCED(1300), TCID(1300) 
DIMENSION TOTRT(1300) 
INTEGER PTSTRT,PTSPST,PNSIZE,SOFF1,SOFF2 
DOUBLE PRECISION MASODL,MOLWT 
REAL LASPTH 
CHARACTER*4 EXTDAT,EXTRAW 
CHARACTER*8 FNAME 
CHARACTER*9 DATE 
CHARACTER* 12 FNAMEX 
CHARACTER*20 TFUEL 
CHARACTER*30 VENT,PROBE 
CHARACTER*60 PBELOC 


ee ee eee Fp 
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C * INITIALIZE VARIABLES 


IQUIT =1 
FIDRNG =1.0E-11 
ISTEP =1 
IWOOD =0 
IAFUEL =0 
SMKVF =0.0 
EXCOEF =0.0 
SMK10 =0.0 
DO 10 N=1,1300 
TCC3(N) =0.0 
SMOKE1(N) = 0.0 
FID(N) =0.0 


FUELRT(N) = 0.0 
FULRT2(N) = 0.0 
RESTIM(N) = 0.0 
Q2(N)  =0.0 
10 CONTINUE 
WETH20 = 0.000 
WETH2 = 0.000 
C * INITIALIZE CONSTANTS 
AREA = 0.0699573 
FIDSPN= 615.0 
PA=713.0 
PRESS=14.3 
VOLC=1.137 
COAMB=0.000000 
CO2AMB=0.000400 
O2AMB=0.232 
eng ke Be oe lL LB 2 Se 72 
C * INPUT RUN PARAMETERS 
WRITE(*,*)'ENTER FILE NAME W/O EXT, 8 CHARACTORS (i.e. HALL0001)' 
READ(*,21)FNAME 
21 FORMAT(A8) 
WRITE(*,*) ' ENTER RUN #' 
READ(*,*) IRUNNO 
WRITE(*,*)' ENTER DATE (ie. "12-23-93") 
READ(*,21) DATE 
WRITE(*,*) 'HEXANE (1), PMMA (2), WOOD (3) OR POLYURETHANE (4) FIR 
SE? 
READ(*,*)IFUELT 
40 GOTO(41,42,43,44)IFUELT 
4] TFUEL="HEXANE' 
ISTEP=1 
DELHC=44735.0 
HTOC=2.333 
STOIC=15.2010 
MOLWT=86.18 
CCT=6.0 
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GOTO 45 
42 TFUEL="PMMA' 
WRITE(*,*)' DO YOU WISH TO KEEP EVERY Ist, 2nd OR 3rd POINT?’ 
WRITE(*,*)' ENTER 1, 2, OR 3 RESPECTIVELY’ 
READ(*,*)ISTEP 
DELHC=25200.0 
HTOC=1.6 
STOIC=8.2824 
MOLWT=100.0 
CCT=5.0 
GOTO 45 
43 TFUEL='WOOD' 
ISTEP=1 
DELHC=17900.0 
HTOC=3.584 
C STOICIOMETRIC RATIO ASSUMES 15% CHAR 
STOIC=3.823 
MOLWT=40.428 
CCT=1.0 
GOTO 45 
44 TFUEL='POLYURETHANE' 
ISTEP=1 
DELHC=26570.0 
HTOC =1.74 
C STOICIOMETRIC RATIO ASSUMES NO CHAR 
STOIC=8.83 
MOLWT=19.91 
CCT=1.0 
C * END OF FUEL SELECTION 
45 WRITE(*,*)TFUEL 
WRITE(*,*) 'ENTER PAN SIZE (6, 8, 9, 11) [inches]' 
READ(*,22) PNSIZE 
22 FORMAT(I2) 
23 FORMAT(A]13) 
WRITE(*,*) 'ENTER FUEL PAN WEIGHT [Kg]' 
READ(*,*) PNWT 
WRITE(*,*) ' ENTER BEGINNING WEIGHT OF FUEL ONLY [Kg]' 
READ(*,*) FOWT 
C** ADDITIONAL FUELS 
C*** WOODEN CEILING INSIDE COMPARTMENT 
WRITE(*,*)' WAS A WOODEN CEILING PRESENT INSIDE THE COMPARTMENT?’ 
WRITE(*,*) ' YES(ENTER 1) AND NO(ENTER 2)' 
READ(*,*) IANS] 
IF (IANS1 .EQ. 1) THEN 
WRITE(*,*) ' ENTER THE INITIAL WEIGHT OF THE WOODEN CEILING', 
+ ' [kg]' 
READ(*,*) FO2WT 
IWOOD = 1 
HCWOOD = 26656.9 
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ENDIF 
C*** AUTOMATED FUEL SUPPLY ACTIVATED 

WRITE(*,*) ' WAS THE AUTOMATED FUEL SUPPLY USED? 

WRITE(*,*) ' YES(ENTER 1) AND NO(ENTER 2)' 

READ(*,*) LANS2 

IF (IANS2 .EQ. 1) THEN 
WRITE(*,*)' ENTER THE INITIAL WEIGHT OF THE RESERVOIR FUEL', 

+ ' [kg]’ 
READ(*,*) FO2WT 
IAFUEL = 1 

ENDIF 


WRITE(*,*) ' ENTER EXHAUST VENT SIZE (ie. 20 X 6.25) [inches]' 
READ(*,23) VENT 

WRITE(*,*) ' ENTER COMPARTMENT/HALLWAY SOFFIT HEIGHT or' 

WRITE(*,*) ' HALL CEILING HEIGHT FROM TOP OF COMPARTMENT EXH VENT" 
WRITE(*,*) ' (0, 4, 8) [inches] 

READ(*,22) SOFF1 

WRITE(*,*) ' ENTER HALLWAY EXIT SOFFIT HEIGHT (0, 4, 8) [inches] 
READ(*,22) SOFF2 


WRITE(*,*) ' ENTER CO RANGE 1 (1000 ppm), 2 (1%) or 3 (10%)' 
READ(*,*) ICORG 
WRITE(*,*) ' ENTER CO2 RANGE 1 (5000 ppm), 2 (2%), 3 (15%) or’, 
+' 4 (20%) 
READ(*,*) ICO2RG 

30 WRITE(*,*) ' (1) COMPARTMENT, (2) EXHAUST DUCT, OR (3) HALLWAY’ 
WRITE(*,*)' ENTER 1, 2 or 3' 
READ(*,*) IPROBE 
GOTO(31,32,33)IPROBE 
GOTO 30 

C COMPARTMENT SAMPLED 

31 PROBE='COMPARTMENT' 
PTSTRT=4 
GOTO 34 

C EXHAUST DUCT SAMPLED 

32 PROBE='EXHAUST’ 

C EXHAUST DUCT 
PTSTRT=14 
GOTO 34 

C HALLWAY SAMPLED 

33 PROBE='HALLWAY' 
PTSTRT=4 

34 WRITE(*,*)PROBE 
WRITE(*,*) 'ENTER PROBE LOCATION (60 CHARACTORS)' 
WRITE(*,*)' (_ ie. EXH DUCT: L=132, H=2, W=0 IN., -OR-' 
WRITE(*,*)' HALL: L=72, H=2(6), W=-1(1) IN. RAKE (SAMPLE) )' 
READ(*,24) PBELOC 

24 FORMAT(A60) 
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WRITE(*,*)' ENTER AMBIENT ATMOSPHERIC TEMPERATURE [sC]' 


READ(*,*)TAMBA 
TAMBA = TAMBA + 273.15 


WRITE(*,*)' ENTER AMBIENT COMPARTMENT TEMPERATURE [sC]' 
WRITE(*,*)' AVERAGE THE TOP 3 COMPARTMENT TC RAKE THERMOCOUPLES' 


READ(*,*) TAMBC 
TAMBC = TAMBC + 273.15 


WRITE(*,*)' ENTER AMBIENT PRESSURE [mm Hg] 


WRITE(*,*)' DEFAULT IS 713.0 mm Hg' 
READ(*,*)PA 


WRITE(*,*)' ENTER THE FID SPAN GAS CONCENTRATION [ppm C2H4]' 


READ(*,*)FIDSPN 
FIDCST=1000000.0 


WRITE(*,*)' ENTER THE FID OPERATING RANGE (for 1E-10, enter 10)' 


READ(*,*)IFIDRG 
FIDRNG=10.0**(-IFIDRG) 


WRITE(*,*)' ENTER THE FID ZERO GAS READING AS [mV]' 


READ(*,*)FIDZER 
FIDZER=FIDZER*FIDRNG 


WRITE(*,*)' ENTER THE FID SPAN GAS READING AS [mV]' 


READ(*,*)FIDCAL 
FIDCAL=FIDCAL*FIDRNG 


WRITE(*,*)' USES A RUN TIME PAST END OF FIRE = 30 sec' 


TMEPST=30.0D0 


WRITE(*,*) ' CONTINUE WITH REDUCTION ? ENTER "0" TO QUIT’ 


READ(*,*) IQUIT 
IF (IQUIT.EQ.0)GOTO 9999 
WRITE(*,*) ' OK, HERE WE GO...' 


er ly aa ta ante aa ciate sce tn 


C * CALCULATIONS BEGIN 


C * ASSIGN THE APPROPRIATE MAXIMUM VALUE FOR CHOSEN ANALYZER RANGE 


C ** O02 ANALYZER 
O2RNG=22.0 
O2CST=100.0 

C ** CO ANALYZER 
GOTO(50,51,52)ICORG 

50 CORNG=1000.0 
COCST=1.000E+6 
GOTO 53 

51 CORNG=1. 
COCST=100. 

GOTO 53 

52 CORNG=10. 
COCST=100. 

C ** CO2 ANALYZER 

53 GOTO(54,55,56,57)ICO2ZRG 

54 CO2RNG = 5000.0 
CO2CST = 1.000E+6 
GOTO 58 
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55 CO2RNG=2. 
CO2CST=100. 
GOTO 58 
56 CO2RNG=15. 
CO2CST=100. 
GOTO 58 
57 CO2RNG=20. 
CO2CST=100. 
Cg eee ee 72 
C * READ DATA AND CREATE AN ARRAY FOR EACH CHANNEL 
C * WRITE TEMPERATURE DATA FILES E AND F 
C # NOTE: I IS THE NUMBER OF DATA POINTS PER CHANNEL 
C OPEN *.RAW DATA FILE 
58 EXTRAW='.RAW' 
FNAMEX=FNAME 
FNAMEX(9:12)=EXTRA W(1:4) 
OPEN(15,FILE=FNAMEX) 
WRITE(*,*)' ' 3 
WRITE(*,*)' READING DATA FILE '.FNAMEX 
C OPEN NEW *.DTE AND *.DTF DATA FILES 
EXTDAT=".DTE’ 
FNAMEX(9:12)=EXTDAT(1:4) 
DO 59 JJ=1,2 
FNAMEX(12:12) = CHAR(68+JJ) 
OPEN(UNIT=244+JJ, FILE=FNAMEX, STATUS="NEW)) 
WRITE(*,*)' AND WRITING DATA FILE '",;FNAMEX 
59 CONTINUE 
I=1 
C BEGINNING OF READING DO LOOP (CURRENTLY 29 DATA BITS / TIME STEP) 
60 READ(15,*,END=65) TIME(I) 
RUNTIM=TIME(I) 
READ(15,*) DRYCO(I) 
READ(15,*) DRYCO2(I) 
READ(15,*) DRYO2(1) 
READ(15,*) FUEL() 
READ(15,*) AIRVEL(I) 
READ(15,*) FID() 
READ(15,*) FUEL2(1) 
READ(15,*) SMOKE1(1) 
READ(15,*) TCH(1) 
READ(15,*) TCH(2) 
READ(15,*) TCH(3) 
READ(15,*) TCH(4) 
READ(15,*) TCH(5) 
READ(15,*) TCH(6) 
READ(15,*) TCH(7) 
READ(15,*) TCH(8) 
READ(15,*) TCH(9) 
READ(15,*) TCC(1) 
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READ(15,*) TCC(2) 
READ(15,*) TCC(3) 
READ(15,*) TCC(4) 
READ(15,*) TCC(5) 
READ(15,*) TCC(6) 
READ(15,*) TCC(7) 
READ(15,*) TCC(8) 
READ(15,*) TCEV(D) 
READ(15,*) TCED(1) 
READ(15,*) TCID(I) 
C TEMPERATURE REDUCTION - CONVERT TEMPERATURES TO KELVIN 
DO 63 JJ=1,8 
TCH(JN=TCH(JN)+273.15 
TCC(IJ=TCC(IN)4273.15 
63 CONTINUE 
TCH(9)=TCH(9)+273.15 
TCEV(I=TCEV(1)+273.15 
TCED(I)=TCED(1)+273.15 
TCID(I)=TCID(1)+273.15 
TCC3(D=TCC(3) 
C WRITE TEMPERATURE DATA TO FILES 
WRITE(25,2005)TIME(I), TCH(1), TCH(2), TCH(3), TCH(4), TCH(5), TCH(6), 
$TCH(7),TCH(8), TCH(9) 
WRITE(26,2006)TIME(I), TCC(1), TCC(2), TCC(3), TCC(4), TCC(5), TCC(6), 
$TCC(7), TCC(8), TCEV(1), TCED(1), TCID(1) 
2005 FORMAT(F7.1,9F6.0) 
2006 FORMAT(F7.1,11F6.0) 
I=1+1 
GOTO 60 
65 I=I-1 
CLOSE (15) 
CLOSE (25) 
CLOSE (26) 
WRITE(*,3001)FNAME,'.RAW READ! ' 
3001 FORMAT(1X,A8,A14) 
WRITE(*,*)' ' 
WRITE(*,3001)FNAME,'DTE WRITTEN!' 
WRITE(*,3001)FNAME,'.DTF WRITTEN" 
WRITE(*,*)' ' 
WRITE(*,*)' TEMPERATURE DATA REDUCED!" 
(OI fanmade a nn a ay 9) 
C * CREATE EXPERIMENT INFORMATION FILE 
EXTDAT=".DTI' 
FNAMEX(9:12)=EXTDAT(1:4) 
OPEN(UNIT=20, FILE=FNAMEX, STATUS="NEW’) 
C * WRITE TO INFO FILE 
WRITE(20,1001)FNAME 
WRITE(20,1002)DATE,IRUNNO,ICORG,ICO2RG,IPROBE, TMEPST 
WRITE(20,1003)PROBE 
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WRITE(20,1004)PA,TAMBA, TAMBC 
WRITE(20,1005)TFUEL 
WRITE(20,1006)IFUELT,ISTEP,DELHC,HTOC,STOIC,MOLWT,CCT 
WRITE(20,1007)FIDSPN,FIDRNG,FIDZER,FIDCAL 
WRITE(20,1008)RUNTIM,FOWT,PNSIZE, VENT,SOFF 1,SOFF2 
WRITE(20,1009)PBELOC 
IF (WOOD .EQ. 1) THEN 

WRITE(20,1010) FO2WT 
ENDIF 
IF (AFUEL .EQ. 1) THEN 

WRITE(20,1011) FO2WT 
ENDIF 

1001 FORMAT(1X,A8) 

1002 FORMAT(1X,A8,I5,312,F9.3) 

1003 FORMAT(1X,A30) 

1004 FORMAT(1X,3F9.3) 

1005 FORMAT(1X,A20) 

1006 FORMAT(1X,213,F10.1,F8.3,F10.4,F10.3,F6.1) 

1007 FORMAT(1X,F12.1,E11.1,2E15.4) 

1008 FORMAT(1X,F7.1,F8.3,13,A14,213) 

1009 FORMAT(1X,A60) 

1010 FORMAT(1X,'A 1/4 INCH THICK, DOUGLAS FIR WOODEN CEILING WAS'/, 
+ 'PRESENT IN THE COMPARTMENT WITH AN INITIAL WEIGHT OF',F6.3, 
+'kg.') 

1011 FORMAT(1X,'THE AUTOMATED FUEL SUPPLY WAS USED WITH AN INITIAL’, 
+/,,FUEL WEIGHT IN THE RESERVOIR OF',F6.3,' kg.') 


CLOSE(20) 
WRITE(*,3001)FNAME,'.DTI WRITTEN!" 
nL 
C * DATA MANIPULATION FOR EACH CHANNEL LOOP J=1 to I 
WRITE(*,*)' ' 


WRITE(*,*)' PLEASE WAIT -- I AM DIGESTING THE DATA' 
PTSPST=NINT(TMEPST/1.97) 
C BEGINING OF ANALYZER DATA CONVERSION DO LOOP 
DO 77 J=1,I-PTSPST 
K1=J+PTSTRT-1 
C ** CONVERT INSTRUMENT MEASUREMENTS FROM VOLTS TO CONCENTRATIONS 
C *** CO ANALYZER [mole fraction] 
IF(DRYCO(K1).LT.0.0)THEN 
DRYCO(K1)=0.0 
ENDIF 
DRYCO(J)=(DRYCO(K1)/5.0)*CORNG/COCST 
C *** CO2 ANALYZER [mole fraction] 
IF(DRYCO2(K1).LT.0.0)THEN 
DRYCO2(K1)=0.0 
ENDIF 
DRYCO2(J)=(DRYCO2(K1)/5.0)*CO2RNG/CO2CST 
C *** O02 ANALZER [mole fraction] 
IF(DRYO2(K1).LT.0.0) THEN 
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DRYO2(K1)=0.0 
ENDIF 
DRYO2(J)=(DRYO2(K1)/5.08)* O2RNG/O2CST 
C *** THC [mole fraction] 
K2 =J+PTSTRT 
74 FID(J)=FID(K2) 
IF (FID(J).LT.0.0) THEN 
FID(J)=0.000 
ENDIF 
FID(J)=FID(J)* 1000.0* FIDRNG 
77 CONTINUE 
WRITE(*,*)) ANALYZER DATA CONVERTED AND COMPENSATED FOR LAG!" 
DO 80 J=1,I 
C *** FUEL WEIGHT [kg] 
FUEL(J)=(FUEL(J)-1.0)* 10.0/4.0 
IF (WOOD .EQ. 1) THEN 
FUEL2(J)=(FUEL2(J)-1.0)* 10.0/4.0 
ENDIF 
C *** INLET AIR VELOCITY [m/s] 
C SHOWN TO BE IDEAL LINEAR TO WITHIN 5% (IE. V=VEL* AREA) 
C EQUATION FROM AIRVEL = (VOLTS*2M/S)/SVOLTS MAX*(760/PA)*(TCID/298) 
AIRVEL(J)=(AIRVEL(J)/5.0D0)*5.098D0*TCID(J)/PA 
ENDIF 
80 CONTINUE 
WRITE(*,*)' FUEL AND AIR DATA CONVERTED AND COMPENSATED FOR LAG"' 
C *** ELIMINATE INCORRECT FUEL DATA POINTS 
DO 85 J=2,I-2 
IF(FUEL(J+1).GT.FUEL(J)) THEN 
FUEL(J+1)=FUEL(J) 
ENDIF 
85 CONTINUE 
IF (IWOOD .EQ. 1) THEN 
DO 86 J=2,I-2 
IF(FUEL2(J+1).GT.FUEL2(J)) THEN 
FUEL2(J+1)=FUEL2(J) 
ENDIF 
86 CONTINUE 
ENDIF 
C *** CALCULATE FUEL BURN RATE USING AVERAGING SMOOTHING [kg/sec] 
DO 90 J=6,(I-5) 
FUELRT(J)=(FUEL(J-5)-FUEL(J+5))/(TIME(J+5)-TIME(J-5)) 
90 CONTINUE 
IF (T(WOOD .EQ. 1) THEN 
DO 92 J=6,(I-5) 
FULRT2(J)=(FUEL2(J-5)-FUEL2(J+5))/(TIME(J+5)-TIME(J-5)) 
92 CONTINUE 
ENDIF ; 
C *** ESTIMATE INITIAL FUEL BURN RATES [kg/sec] 
DO 95 J=1,5 
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FUELRT(J)=FUELRT(8) 
95 CONTINUE 
IF (IWOOD .EQ. 1) THEN 
DO 96 J=1,5 
FULRT2(J)=FULRT2(8) 
96 CONTINUE 
ENDIF 
C *** ESTIMATE LAST FUEL BURN RATES [kg/sec] 
DO 97 J=(I-4),(1-2) 
FUELRT(J)=(FUEL(J-2)-FUEL(J+2))/(TIME(J+2)-TIME(J-2)) 
97 CONTINUE 
IF (IWOOD .EQ. 1) THEN 
DO 98 J=(I-4),(I-2) 
FULRT2(J)=(FUEL2(J-2)-FUEL2(J+2))/(TIME(J+2)-TIME(J-2)) 
98 CONTINUE 
ENDIF 
DTIME=TIME(I)-TIME(I-2) 
FUELRT(I-1)=(FUEL(I-2)-FUEL(I)/DTIME 
IF WOOD .EQ. 1) THEN 
FULRT2(I-1)=(FUEL2(I-2)-FUEL2(1))/DTIME 
ENDIF 
DTIME2=TIME(1)-TIME(I-1) 
FUELRT(I)=(FUEL(I-1)-FUEL(1))/DTIME2 
IF (WOOD .EQ. 1) THEN 
FULRT2(1)=(FUEL2(I-1)-FUEL2(1))/DTIME2 
ENDIF 
DO 100 J=(I-PTSPST+1),] 
FUELRT(J)=0.0 
IF (IWOOD .EQ. 1) THEN 
FULRT2(J)=0.0 
ENDIF 
100 CONTINUE 
WRITE(*,*)' BURN RATE CALCULATIONS COMPLETE’ 
C CALCULATE SMK10 FROM 1ST 5 SMOKEI'S SINCE SMOKE0 NOT WORKING 
SMK10=0.000 
DO 110 J=1,5 
SMK10 = SMK10 + SMOKE1(J)/S. 
110 CONTINUE 
(ee Pe Ee 72 
C * CREATE OUTPUT FILES 
EXTDAT='.DTA' 
FNAMEX(9:12)=EXTDAT(1:4) 
DO 140 J=1,4 
FNAMEX(12:12) = CHAR(64+J) 
OPEN(UNIT=20+J, FILE=FNAMEX, STATUS="NEW’) 
140 CONTINUE 
(Ge Ny PR SRE a a A as 72 
C 1 TIME INITIALIZATIONS BEFORE LOOP (150) 
SMKMAX = 0.0 
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DO 150 K=1,I-PTSPST,ISTEP 
C * CALCULATE HEAT RELEASE RATE (kW) 
Q(K) = FUELRT(K)*DELHC 
IF (IWOOD .EQ. 1) THEN 
Q2(K) = FULRT2(K)*HCWOOD 
ENDIF 
Cc ae cnet ee oe ee ip 
C * CALCULATE AIR FLOWRATE (kg/sec) 
C AND RESIDENCE TIME (sec) 
DENSTY=(PA/TCID(K/0.287)*0.1333 
IF (AIRVEL(K).LE. 0.0) THEN 
WRITE(*,*)' AIR VELOCITY ERROR ! (AIRVEL < 0) 
RESTIM(K)=999. 
ELSE 
RESTIM(K)=VOLC/(AIRVEL(K)*AREA*TCC3(K)/TCID(K)) 
ENDIF 
AIRRAT(K)=AIRVEL(K)*DENSTY* AREA 
(gotta Resweiceobin Atinao tracy Rar ste es RARE AB cy 72 
C * CALCULATE FUEL/AIR RATIO 
IF (AIRRAT(K).LE.0.0) THEN 
FTOA(K) = 0.0 
ELSE 
IF (IWOOD .EQ. 1) THEN 
TOTRT(K)=FUELRT(K)+FULRT2(K) 
ELSE 
TOTRT(K)=FUELRT(K) 
ENDIF 
FTOA(K)=TOTRT(K)/AIRRAT(K) 
ENDIF 


ty ncn lat cect aso ee BB 


C *** CALCULATE WET CONCENTRATIONS, YIELDS AND STORE DATA IN FILES 
Co Tee 
C * VOLUMETRIC FLOWRATE (m‘3/s) IN EXHAUST DUCT 
C EXMOLES ASSUMES AVG MOLECULAR WEIGHT EXHAUSTED IS SAME AS AIR 
EXMOLS = (AIRRAT(K)+FUELRT(K))/28.97 
DCTFLW = 1.2099*DSQRT(PRESS*TCED(K)/PA) 
TDMOLS = DCTFLW*PA*0.1333/(8.3144*TCED(K)) 
Cr eer renee TS 
C * FOR EXHAUST DUCT SAMPLED FIRES 
IF (IPROBE.EQ.2) THEN 
WETCO=DRYCO(K) 
WETCO2=DRYCO2(K) 
WETO2=DRYO2(K) 
ELSE 
C * FOR (UNDILUTED) HALLWAY AND COMPARTMENT SAMPLED FIRES 
DENOM = 1.0+((HTOC/2.0)*DRYCO2(K)) 
WETCO=DRYCO(K)/DENOM 
WETCO2=DRYCO2(K)/DENOM 
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WETO2=DRYO2(K)/DENOM 
ENDIF 
Cnn] SO en 6 ee 79 
C SET YIELDS TO 0 IF FUELRATE = 0 
C OR IF HALLWAY SAMPLED 
IF (FUELRT(K).EQ.0.0.0R.IPROBE.EQ.3) THEN 
COYLD=0.0 
CO2YLD=0.0 
O2YLD=0.0 
ELSE IF (IPROBE.EQ.2) THEN 
C EXHAUST DUCT SAMPLED 
COYLD=(WETCO-COAMB)*TDMOLS*28.01/FUELRT(K) 
CO2YLD=(WETCO2-CO2AMB)*TDMOLS*44.01/FUELRT(K) 
O2YLD=0.0 
ELSE 
C COMPARTMENT SAMPLED 
COYLD=(WETCO-COAMB)*EXMOLS*28.01/FUELRT(K) 
CO2YLD=(WETCO2-CO2A MB)*EXMOLS*44.01/FUELRT(K) 
O2YLD=((O2AMB*AIRRAT(K))-(WETO2*EXMOLS*32.00))/FUELRT(K) 
ENDIF 
Co eee ere ee ee 
C CHECK FOR YIELD ERRORS 
IF (CO2YLD.GT.100.0) THEN 
WRITE(*,*)' ERROR IN CO2 YIELD! > 100° 
ENDIF 
IF (COYLD.GT.100.0) THEN 
WRITE(*,*)' ERROR IN CO YIELD! > 100! 
ENDIF 
IF (02YLD.GT.100.0) THEN 
WRITE(*,*)' ERROR IN 02 YIELD! > 100' 
END IF 
C WRITE(*,*)' WET CONCENTRATIONS AND YIELDS COMPLETE’ 
ee 72 
C * CALCULATE SMOKE VOLUME FRACTION IN [ppb] 
IF (SMOKEI(K).EQ.0.)THEN 
WRITE(*,*)' ERROR IN SMOKE VOLUME FRACTION CALCULATION !' 
WRITE(*,*)' SMOKE1 SIGNAL = 0!" 
EXCOEF= 999, 
SMKVF = 999. 
SMKYLD= 999. 
ELSE 
LASPTH = 0.4572 
IF (SMOKEI(K).GE.SMK10) THEN 
EXCOEF = 0.000 
ELSE 
EXCOEF = DLOG(SMK10/SMOKE1(K))/LASPTH 
ENDIF 
SMKVF = 1.3697E-07*EXCOEF* 1.0E+9 
IF (SMKVF.LT.0.0) THEN 
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SMKVF = 0.0 
ENDIF 
C * CALCULATE SMOKE YIELD ACCORDING TO TEWARSON 
IF (FUELRT(K).LE.0.0)THEN 


WRITE(*,*)' ERROR IN SMOKE YIELD CALCULATION !' 
WRITE(*,*)'FUELRATE = 0 ='",;FUELRT(K),'AT TIME = '", TIME(K) 
SMKYLD = 999. 
ELSE 
GF = FUELRT(K)*1000. 
ODL = EXCOEF/2.303 
MASODL = ODL * DCTFLW / GF 
SEC = 3.213/(0.67*1.1) 
SMKYLD = MASODL / SEC 


ENDIF 
ENDIF 
C  WRITE(*,*)' SMOKE CALCULATIONS COMPLETE’ 


C * EQUIVALENCE RATIO 
EQUIV=FTOA(K)*STOIC 
IF (EQUIV.LE.0.0)THEN 
EQUIV=0.0 


C * CALCULATE DERIVATIVE OF FUEL BURN RATE 
IF (K.LE.5.OR.K.GE.I-PTSPST-4) THEN 


D2FUEL(K)=0.0 
ELSE 
D2FUEL(K)=(FUELRT(K+5)-FUELRT(K-5))/(TIME(K+5)-TIME(K-5)) 
ENDIF 
fa ciate a a ita i ich th hia Lip ns oa =72 


C * STEADY STATE TIME RATIO 
IF (FUELRT(K) .EQ. 0.0) THEN 
FUELRT(K) = 0.00001 
ENDIF 
SSTIME=RESTIM(K)*D2FUEL(K)/FUELRT(K) 
IF (ABS(SSTIME).GT.10.0) THEN 
SSTIME = 10.0 


C * CALCULATE UNBURNED HYDROCARBON CONCENTRATIONS AND YIELDS 
IF (FID(K).LE.0.0)THEN 
THC=0.0 
THCYLD=0.0 
ELSE 
THC=(FIDSPN/FIDCST)*(FID(K)-FIDZER)/(FIDCAL-FIDZER) 
IF (IPROBE.EQ.1)THEN 
THCYLD=(THC*(EXMOLS)*28.05)/FUELRT(K) 
ELSE IF (IPROBE.EQ.2) THEN 
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THCYLD=(THC*(TDMOLS)*28.05)/FUELRT(K) 
ELSE 
THCYLD=0.0 


C * CARBON BALANCE ERROR CHECK 
IF (IPROBE.EQ.3) THEN 
C  _WRITE(*,*)' OR SAMPLED IN HALLWAY ' 
CERROR=999.0 
GOTO 145 
ENDIF 
IF (FUELRT(K).LE.0.0) THEN 
WRITE(*,*)' ERROR IN CARBON BALANCE CHECK"' 
WRITE(*,*)' FUELRATE = 0 AT TIME = '", TIME(K) 
CERROR=999.0 
GOTO 145 
ENDIF 
fae pe ONES Oe ly 9) 
C * MOLES OF CARBON IN FROM FUEL 
CMOLIN=CCT*FUELRT(K)/MOLWT 
IF (CMOLIN.LE.0.0)THEN 
WRITE(*,*)' ERROR IN CARBON BALANCE CHECK W/ CMOLIN !" 
CERROR=999.0 
GOTO 145 


(C(----- Joann on nn nnn nn nnn nnn mene nn nnn nnn nnnnn nn nnn nnn nn nnnnncnnnncns 72 
C * MOLES ACCOUNTED FOR BY MEASUREMENTS 
C IF COMPARTMENT SAMPLED 
IF (IPROBE.EQ.1) THEN 
CMOLOT=(WETCO+WETCO2+2.0*THC)*EXMOLS 
ELSE 
C IF EXHAUST DUCT SAMPLED 
CMOLOT=(WETCO-COAMB+WETCO2-CO2AMB+2.*THC)*TDMOLS 
ENDIF 
CERROR=(CMOLOT-CMOLIN)* 100.0/CMOLIN 
145 IF (ABS(CERROR).GT.999.1) THEN 


CERROR=999.0 
ENDIF 
C WRITE(*,*)' CARBON ERROR CALCULATIONS COMPLETE’ 
(qutet, REE EAReeee y 9) 


C * STORE DATA (80 CHARACTORS / LINE MAX) 
WRITE(21,2001)TIME(K), WETCO,WETCO2,WETO2,COYLD,CO2YLD,O2YLD, 
$SMKYLD,THCYLD 
WRITE(22,2002)TIME(K),FUEL(K),FUELRT(K),AIRVEL(K),AIRRAT(K), 
$D2FUEL(K),CERROR,EQUIV 
WRITE(23,2003)TIME(K),RESTIM(K),SSTIME,DRYCO(K),DRYCO2(K), 
$DRYO2(K), THC,Q(K), TDMOLS/EXMOLS 
WRITE(24,2004) TIME(K),SMK VF,SMOKE1(K),EXCOEF, WETH20, WETH2, 
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$FUEL2(K), FULRT2(K),Q2(K) 
2001 FORMAT(F6.1,3E10.3,F7.4,2F10.4,F9.3,F7.4) 
2002 FORMAT(F6.1,F10.3,F10.4,F10.2,F10.3,F9.5,F11.2,F8.3) 
2003 FORMAT(F6.1,F7.1,F7.2,3E10.3,F10.6,F8.2,F8.1) 
2004 FORMAT(F6.1,F8.2,6F9.3,F10.2) 
150 CONTINUE 
CLOSE (21) 
CLOSE (22) 
CLOSE (23) 
CLOSE (24) 
WRITE(*,*)' DATA REDUCTION COMPLETED!’ 
WRITE(*,*)' ' 
WRITE(*,*)' ' 
WRITE(*,*)' FILES WRITTEN:’ 
WRITE(*,*)' 
WRITE(*,3000)FNAME,'.DTI' 
WRITE(*,3000)FNAME,'.DTA' 
WRITE(*,3000)FNAME,'.DTB' 
WRITE(*,3000)FNAME,'.DTC’ 
WRITE(*,3000)FNAME,'.DTD' 
WRITE(*,3000)FNAME,'.DTE' 
WRITE(*,3000)FNAME,'.DTF' 
3000 FORMAT(1X,A8,A4) 
9999 STOP 
END 
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FIREAVG3.FOR 

C PROGRAM: *** FIREAVG3.FOR *** 

C REVISED BY: BRIAN LATTIMER LAST UPDATE: 08/31/94 

C WRITTEN BY: DAVID EWENS LAST UPDATE: 11/17/93 

Cc 

C PURPOSE: 

C THIS PROGRAM FINDS THE AVERAGE VALUES FOR THE VARIABLES OVER A GIVEN 

C TIME RANGE IN THE DATA FILES CREATED BY PROGRAM FIRERED2.FOR. 

C THIS PROGRAM ALSO CALCULATES 95% C.I. FOR EACH VARIABLE (+/- 2 SIGMA) 

C OUTPUTS AVG, % ERROR BAND (FROM C.I.), 95% C.I. BAND 

Cc 

C VARIABLES: 

Cc 

C FILES ACCESSED: 

C (HALL####.*, WHERE #### IS THE RUN NUMBER, ie. HALLO0001 ) 

Cc 

C READ: 

CG 

C HALL###4.DTA - EXPERIMENTAL CHRONOLOGICAL DATA FILE (CONTAINS: TIME, 
WET CONC'S, YIELDS) 

HALL####.DTB - EXPERIMENTAL CHRONOLOGICAL DATA FILE (CONTAINS: TIME, 
FUEL & INLET AIR DATA, CARBON ERROR, EQUIVALENCE RATIO) 

HALL##4#4.DTC - EXPERIMENTAL CHRONOLOGICAL DATA FILE (CONTAINS: TIME, 
TIME SCALES, DRY CONC'S, THC, HEAT RELEASE, DILUTION) 

HALL####.DTD - EXPERIMENTAL CHRONOLOGICAL DATA FILE (CONTAINS: TIME, 
SMOKE DATA, WET H20 & H2 CONC'S) 

HALL####.DTE - EXPERIMENTAL CHRONOLOGICAL DATA FILE (CONTAINS: TIME, 
HALLWAY THERMOCOUPLE DATA) 

HALL####.DTF - EXPERIMENTAL CHRONOLOGICAL DATA FILE (CONTAINS: TIME, 
COMPARTMENT THERMOCOUPLE DATA, INLET AIR & EXHAUST 
THERMOCOUPLE DATA) 

C WRITTEN: 

6 

C TOPRINTER - AVG AND 95% C.I. FOR ALL FILES ABOVE 

C 

C2345-7-9012345678901234567890 1234567890 1234567890 1234567890 123456789 72 

C * INITIALIZATION 

IMPLICIT DOUBLE PRECISION (A-H,O-Z) 
DOUBLE PRECISION AVG(15),A(15,200),A2SIG(15),APER(15) 
CHARACTER*4 EXTDAT 
CHARACTER*8 FNAME 
CHARACTER*10 CNAME(12) 
CHARACTER* 12 FNAMEX 
C CHARACTER*12 FDAT 
CHARACTER*27 PATHNM 


OOO GG Gr Gara: O).cy O 


C * INPUT FILE NAME 
C  WRITE(*,*)' IS THE FOLLOWING PATH TO DATA FILES CORRECT ?' 
C  WRITE(*,*)' YES = ENTER ANY NUM, NO = CHANGE PATHNM IN PROGRAM ' 
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C MAX 15 CHAR'S, USE ALL 15, PUT BLANKS AS FRONT CHAR'S, NOT LAST CHAR'S! 
C >>>>>>>>123456789012345<<<<<< 
PATHNM="' 
C WRITE(*,*)PATHNM 
C READ(*,*)DUMMY 
WRITE(*,*)' ENTER THE FILE NAME TO AVERAGE" 
WRITE(*,*)' WITH OUT EXTENSION (ie. HALL0001)' 
READ(*,10)FNAME 
C  WRITE(*,*)PATHNM,FNAME 
10 FORMAT(A8) 
C  WRITE(*,*)' ENTER THE TEST NUMBER' 
C  READ(*,*)ITEST 
WRITE(*,*)' ENTER START TIME FOR AVERAGE! 
READ(*,*)START 
END=START+60.0D0 
C  WRITE(*,*)' ENTER END TIME FOR AVERAGE!’ 
C READ(*,*)END 
C OPEN PRINTER FILE 
EXTDAT="'.AVG' 
FNAMEX=FNAME 
FNAMEX(9:12)=EXTDAT(1:4) 
C OPEN(UNIT=10, FILE=FNAMEX, STATUS="NEW’) 
OPEN(UNIT=10, FILE='PRN') 
C SET UP PRINTER PAPER 
WRITE(10,30)' AVERAGE TIME SPAN (SEC):',START,' - END 
30 FORMAT(4X,A25,F7.1,A3,F7.1) 
C SET UP FOR READ FILES 
EXTDAT=".DTA' 
FNAMEX(9:12)=EXTDAT(1:4) 
PATHNM(16:27)=FNAMEX(1:12) 
WRITE(*,*)FNAMEX,PATHNM 
tl es et Le tan i et 72 
C * READ DATA 
C * WRITE TEMPERATURE DATA FILES E AND F 
C # NOTE: I IS THE NUMBER OF DATA POINTS PER CHANNEL 
C OPEN *.RAW DATA FILE 
DO 400 I=1,6 
FNAMEX(12:12) = CHAR(64+1) 
PATHNM(27:27) = CHAR(64+1) 
WRITE(*,*)FNAMEX,PATHNM 
WRITE(*,*)' ' 
WRITE(*,*)' READING DATA FILE '"PATHNM 
WRITE(*,*)' ' 
OPEN(15,FILE=PATHNM) 
DO 40 KK=1,15 
AVG(KK)=0.0D0 
DO 39 KKK=1,50 
A(KK,KKK)=0.0D0 
39 CONTINUE 
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40 CONTINUE 
C KOUNT=0 
K=1 
C * READ DATA FROM CORRECT FILE 
50 GOTO (51,52,53,54,55,56)I 


C 51 READ(21,2001)TIME,WETCO,WETCO2, WETO2,COYLD,CO2YLD,O2YLD, 
C $SMKYLD,THCYLD 


gy feed Ne al ek ent on tren eel B) 

51 READ(15,2001,END=70)TIME,A(1,K),A(2,K),A(3,K),A(4,K),A(5,K),A(6,K) 
$,A(7,K),A(8,K) 
JJ=8 
CNAME(1)='WET CO' 


CNAME(2)="WET CO2' 
CNAME(3)='WET 02' 
CNAME(4)='CO YIELD' 
CNAME(5)='CO2 YIELD' 
CNAME(6)='02 YIELD’ 
CNAME(7)="SMK YIELD’ 
CNAME(8)="THC YIELD' 
GOTO 57 


C 52 READ(22,2002)TIME,FUEL,FUELRT,AIRVEL,AIRRAT, 
C $D2FUEL,CERROR,EQUIV 
Ch Lo ee ea aE, 5 
52 READ(15,2002,END=70)TIME,A(1,K),A(2,K),A(3,K),A(4,K),A(5,K),A(6,K) 
$,A(7,K) 
JJ=8 
A(8,K)=A(2,K)+A(4,K) 
CNAME(1)='FUEL WT" 
CNAME(2)='FUEL RATE’ 
CNAME(3)='AIR VEL' 
CNAME(4)='AIR RATE' 
CNAME(5)='D2FUEL' 
CNAME(6)='C ERROR' 
CNAME(7)="EQ. RATIO’ 
CNAME(8)='VENT FLOW' 
GOTO 57 


C 53 READ(23,2003)TIME,RESTIM,SSTIME,DRYCO,DRYCO2, 
C $DRYO2,THC,Q,TDMOLS/EXMOLS 
ea pen het el ori Re ene =72 
53 READ(15,2003,END=70)TIME,A(1,K),A(2,K),A(3,K),A(4,K),A(5,K),A(6,K) 
$,A(7,K),A(8,K) 
JJ=8 
CNAME(1)="RES TIME' 
CNAME(2)='SS TIME’ 
CNAME(3)="DRY CO' 
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CNAME(4)="DRY CO2' 
CNAME(5)="DRY 02' 
CNAME(6)='THC' 
CNAME(7)='Q, HEAT’ 
CNAME(8)='DIL RATIO’ 
GOTO 57 


C 54 READ(24,2004)TIME,SMK VF,SMOKE0,SMOKE1,EXCOEF,WETH20, 
C $WETH2 
(ely faecal nln, Fe ea BF 
54 READ(15,2004,END=70)TIME,A(1,K),A(2,K),A(3,K),A(4,K),A(5,K),A(6,K) 
$,A(7,K),A(8,K) 
JJ=8 
CNAME(1)='SMK VF' 
CNAME(2)='SMOKE 1' 
CNAME(3)="EXCOEF' 
CNAME(4)='WET H20' 
CNAME(5)="WET H2' 
CNAME(6)='FUEL WT2' 
CNAME(7)="FUEL RATE2' 
CNAME(8)='Q2, HEAT2' 
GOTO 57 


C 55 READ(25,2005)TIME(I), TCH(1), TCH(2), TCH(3), TCH(4), TCH(S), TCH(6), 
C $TCH(7),TCH(8),TCH(9) 


Cece cen 
55 READ(15,2005,END=70)TIME,A(1,K),A(2,K),A(3,K),A(4,K),A(5,K),A(6,K) 
$,A(7,K),A(8,K),A(9,K) 
JJ=9 


CNAME(1)='TC HALL 1' 
CNAME(2)='TC HALL 2' 
CNAME(3)="TC HALL 3' 
CNAME(4)='TC HALL 4' 
CNAME(S5)='TC HALL 5" 
CNAME(6)="TC HALL 6' 
CNAME(7)="TC HALL 7' 
CNAME(8)="TC HALL 8' 
CNAME(9)='TC HALL 9 
GOTO 57 


C 56 READ(26,2006)TIME(I), TCC(1), TCC(2), TCC(3), TCC(4), TCC(S), TCC(6), 
C $TCC(7),TCC(8), TCEV(I), TCED(),TCID() 
Cm ee 72 
56 READ(15,2006,END=70)TIME,A(1,K),A(2,K),A(3,K),A(4,K),A(5,K),A(6,K) 
$,A(7,K),A(8,K),A(9,K),A(10,K),A(11,K) 
J=11 
CNAME(1)="TC COMP 1' 
CNAME(2)='TC COMP 2' 
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CNAME(3)='TC COMP 3' 
CNAME(4)='TC COMP 4' 
CNAME(5)='TC COMP 5' 
CNAME(6)='TC COMP 6' 
CNAME(7)='TC COMP 7' 
CNAME(8)='TC COMP 8' 
CNAME(9)='"TC EX VENT 
CNAME(10)="TC EX DUCT 
CNAME(11)="TC IN DUCT 


57 IF (TIME.LT.START) GOTO 50 
IF (TIME.GT.END) GOTO 70 
C KOUNT=KOUNT+1 
K=K+1 
DO 60 J=1,JJ 
AVG(J) = AVG(J)+A(,K-1) 
60 CONTINUE 
GOTO 50 


Fray nt ort pn, an er enhnale 72 
70 WRITE(10,*)FNAMEX;' AVERAGE +/-% ERROR +/- 95% CI. 
$ % ERROR (2 SIGMA)’ 
WRITE(*,*)FNAMEX;' AVERAGE +/-% ERROR +/- 95% C.L. 
$% ERROR (2 SIGMA) 
75 FORMAT(2X,A10,2X,12,2X,F16.10,2X,F9.2,2X,D12.6) 
K=K-1 
DO 80 J=1,JJ 
C CALCULATE AVG AND 2 SIGMA 
A2SIG(J)=0.0D0 
AVG(J}=AVG(J)/DBLE(FLOAT(K)) 
DO 78 KS=1,K 
A2SIG(J)=A2SIG(J+((A(J,KS)-A VG(J))**2) 
78 CONTINUE 
A2SIG(J)=A2SIG(J)/DBLE(FLOAT(K-1)) 
A2SIG(J)}=2.0D0*SQRT(A2SIG(J)) 
APER(J)=0.0D0 
IF (ABS(AVG(J)).GE.1.0D-9) THEN 
APER(J)=100.0D0*A2SIG(J)/AVG(J) 
ENDIF 
C WRITE OUTPUT 
WRITE(10,75)CNAME(J),J,A VG(J), APER(J),A2SIG(J) 
WRITE(*,75)CNAME(J),J,A VG(J),APER(J),A2SIG(J) 
80 CONTINUE 
C WRITE(10,*)' ' 
CLOSE(15) 


(C7111) 7 III III IO IOI KT 
400 CONTINUE 
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2001 FORMAT(F6.1,3E10.3,F7.4,2F10.4,F9.3,F7.4) 
2002 FORMAT(F6.1,F10.3,F10.4,F10.2,F10.3,F9.5,F11.2,F8.3) 
2003 FORMAT(F6.1,F7.1,F7.2,3E10.3,F10.6,F8.2,F8.1) 

2004 FORMAT(F6.1,F8.2,6F9.3,F10.2) 

2005 FORMAT(F7.1,9F6.0) 

2006 FORMAT(F7.1,11F6.0) 


WRITE(10,*)' ' 
WRITE(*,*)' ' 
WRITE(*,*)' DATA REDUCTION COMPLETED!" 
WRITE(*,*)' ' 
WRITE(*,*)' FILE PRINTED!" 
WRITE(*,*)' ' 
9999 STOP 
END 
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Appendix C. Enhanced CO Generation in the Compartment 


APPENDIX C. 


ENHANCED CO GENERATION IN THE COMPARTMENT 


C.1 INTRODUCTION 


Two experimental situations were investigated where CO levels inside the 
compartment were measured to be significantly higher than levels measured in the study 
by Gottuk (1992a). The effects of bi-directional flow opening (a door) connecting the 
compartment and hallway on CO levels in the compartment and hallway is initially 
discussed. The effect of a wood ceiling inside the compartment on the CO levels 
measured in the compartment and hallway is then discussed. Results of these two 
situations are compared with results from experiments containing a 0.12 m’ opening with 


unidirectional flow. 


C.2 DOORWAY EXPERIMENTS 
C.2.1 Experimental Methodology 


A study was performed to investigate the effect of the presence of a bi-directional 
flow opening on the species concentrations. To produce this using the current 
compartment, shown in Fig. 2.2, the plenum was blocked off from both the ambient air 
and the inside of the compartment, see Fig. C.1. Note that the compartment-hallway 
configuration is different than experiments previously discussed. This compartment- 
hallway orientation is the same as that used in the study by Ewens (1994a). The opening 
connecting the compartment and hallway was enlarged to its maximum dimensions of 
0.50 m wide, 0.75 m high (0.38 m’). The experiments were conducted with a 0.20 m 
inlet soffit and a 0.23 m diameter fuel pan in the room. Comparisons are made with 


experiments having a 0.12 m* opening, a 0.20 m inlet soffit, and a 0.20 m diameter fuel 
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pan. A different fuel pan diameter was used so the experiments would have 
approximately the same fire size, 576 kW with a door and 542 kW with the 0.12 m’ 


opening. 


Sampling was conducted at two locations in the compartment and three locations 
along the hallway with all sampling being performed 0.05 m below the ceiling. The in- 
compartment sampling was performed 0.20 m from the wall containing the door on the 
temperature rake side of the compartment (termed the front), and 1.02 m from the wall 
containing the door on the same side of the compartment (termed the back). The in- 
hallway sampling locations were performed 0.45 m, 1.80 m and 3.20 m downstream of 
the compartment opening as shown in Fig. C.1. The door (or opening) is taken as the 


spatial origin in the hallway and is shown in the plots as 0.0 m. 
C.2.2 Flow Visualization 


The differences in the flow of hallway gases in the experiments with a door and a 
0.12 m? opening is expected to be significantly different. The external burning in the 


hallway was observed to determine the differences in the flow patterns. 


With a door connecting the compartment and the hallway, two large vortical 
structures were present in the first half of the hallway as shown in Fig. C.2 (a). The 
vortical structures originated near the bottom of the door approximately 0.9 m below the 
ceiling. The structures were observed to rise nearly vertically in the corridor until the 
flame impinged on the ceiling. These vortical structures were also observed in the 0.12 
m’ opening experiments, but were much smaller, see Fig. C.2 (b). With the 0.12 m* 
opening, the gases were only permitted to rise 0.44 m before impinging upon the ceiling. 
In both cases, the gases were seen to travel down the hallway as a ceiling jet after 
impinging upon the ceiling. Vortical structures were also evident downstream of the 


compartment along the corner of the ceiling and the walls in both cases. The vortical 
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(b) 


Figure C.2 The external burning n-hexane fire experiments with (a) a door and (b) a 
0.12m* opening connecting the compartment to the hallway. 
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structures do not separate as observed by Hinkley et al. (1984), but are instead connected 


by a sheet of fire. 


C.2.3. Results and Discussion 


Axial profiles of the combustion gases from experiments with a door and with a 
0.12 m” opening connecting the compartment to the hallway are shown in Fig. C.3 and 
Fig. C.4, respectively. The normalizing values are shown in parentheses beside the 
appropriate species in the caption. Negative axial locations indicate sampling inside the 


compartment. 


The levels of CO inside the compartment with a door are seen in Fig. C.3 to be 6.4% 
in the back and 4.5% in the front. This is a significant increase over the levels measured 
inside the compartment with a 0.12 m opening of 2.8% in the back and 3.7% in the front, 
see Fig. C.4. The difference in the two experiments is not only the path of the air 
entrainment into the compartment, but also the momentum of the air entering the 
compartment. When air is drawn into the compartment through the plenum from two 
sides of the plume, the plume (with a 0.12 m? opening) appears to be relatively 
undisturbed, and rises vertically in the compartment, see Fig. C.5 (a). With the air 
entering through the door, the momentum of the air is apparent by the flame from the 
pool fire being bent backward. The flame eventually was observed to impinge on the 
back wall of the compartment as shown in Fig. C.5 (b). It is believed that the air which 
does not get entrained into the plume impinges upon the back wall. The air is believed to 
then be drawn upward into the compartment upper-layer. This “injection” of air into the 
hot, fuel rich environment of the compartment upper-layer results in the oxidation of 


UHC leading directly to the formation of additional CO. 


Concentrations of CO and UHC are both seen in Fig. C.3 as being oxidized rapidly 
in the first portion of the hallway, but the cool temperatures in the later 1.80 m of the 


hallway prevented further oxidation. Significant dilution of the gases was, however, 
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Figure C.3 The normalized species concentrations (wet) and temperatures within the 
compartment and hallway for n-hexane fire experiments with a door. 
Symbols and normalization values: @-CO (8.0%), HI-UHC (10%), A-CO, 
(20%), W-O0,(10.5%), and O-Temperature. 
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Figure C.4 The normalized species concentrations (wet) and temperatures within the 
compartment and hallway for n-hexane fire experiments with a 0.12 m” 
opening. Symbols and normalization values: @-CO (8.0%), H-UHC (10%), 
A-CO, (20%), W-0,(10.5%), and O-Temperature. 
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(b) 


Figure C.5 The effect of air being entrained into the compartment through (a) the plenum 
and (b) the doorway on the fire plume inside the compartment. 
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measured to occur in the later portion of the hallway. This is evident by the rise in the O, 
concentrations to approximately 7.5% by the end of the hallway. Oxidation of CO and 
UHC is also significant in the first half of the hallway in experiments with a 0.12 m* 
opening, but CO oxidation was measured to be much slower than that measured in the 
door experiments. With a 0.12 m’ Opening, temperatures in the hallway were measured 
to remain above 900 K the entire length of the hallway. This allowed CO oxidation to 
occur for the majority of the length. The O, concentration of 2% in the upper-layer at the 
end of the hallway was measured to be considerably lower compared to the 7.5% 


measured in the experiments with a door. 


The post-hallway yields of CO and UHC with a door were 0.10 and 0.029, 
respectively. In experiments with a 0.12 m’ opening, the post-hallway yield of CO of 
0.12 and UHC of 0.038 are slightly higher than those measured with a door present. The 
reduction in the concentration of CO and UHC at the end of the hallway was, however, 


higher in the experiments with a door. The concentration reduction 


Thea -(1-( Fata na) )r00 [%] (C.1) 
i/comp 


in CO and UHC was 95% and 99%, respectively for experiments with a door and 85% 
and 96%, respectively, with a 0.12 m’ opening. Noting that approximately the same 
amount of CO and UHC leave the hallway, the higher concentration reduction efficiency 
with a door present is the result of a dilution of the combustion gases via increased air 


entrainment into the upper-layer, which is unaccounted for in the yield data. 


C.2.4 Summary and Conclusions 


Experiments with a door resulted in approximately the same post-hallway CO and 
UHC yields compared with levels measured with a 0.12 m” opening. Experiments with a 
door resulted in 50% lower CO and UHC concentrations and 300% higher O, 


concentration exiting the hallway compared with measurements in experiments with a 
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0.12 m? opening. The higher dilution in the door case is attributed to more efficient 
entrainment into the compartment fire gases exiting the compartment and rising to the 


ceiling. 


C.3 COMPARTMENT WITH A COMBUSTIBLE CEILING 


C.3.1 Background 


In 1987, three people died due to smoke inhalation in a townhouse in Sharon, 
Pennsylvania. Extremely high levels of carboxyhemoglobin, 91%, were present in one of 
the victim’s bloodstream. This prompted an investigation by NIST/BFRL to simulate the 
townhouse fire (Levine and Nelson, 1990). The source of the fire was in the kitchen of 
the townhouse where a large amount of wood (wood paneling and cabinets) was located. 
The CO levels exiting the kitchen were found to be as high as 8.5%-dry, while the levels 
upstairs where the bedrooms were located was 5.0%-dry (Levine and Nelson, 1990). 
Experiments performed by Pitts et al. (1994a) in a reduced-scale compartment with the 
ceiling and the upper-portion of the walls lined with wood showed the CO concentrations 
were 6%-dry or greater in the front and 12%-dry in the rear of the compartment. This 
was a dramatic increase from the CO concentrations, approximately 4%-dry in the front 
and 3%-dry in the rear, measured in a compartment without combustible surfaces (Pitts, 


1994a). 


C.3.2 Experimental Methodology 


The compartment and hallway were oriented in the same manner shown in Fig. C.1. 
A 0.20 m soffit was present at the hallway entrance while no soffit was present at the exit. 
Three sheets of 1.22 m long, 0.30 m wide, 6.35 mm thick Douglas fir plywood were 
suspended approximately 0.05 m below the ceiling of the compartment, see Fig. C.6. The 


mass loss rate of the wood was monitored using a A&D load cell located on the top of the 
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compartment. During the fire, both sides of the wood (total surface area of 2.20 m’) were 
exposed to the hot upper layer gases and were allowed to pyrolyze. This is less surface 
area than was used by in the experiments performed by Pitts, et al. (1994) where 3.19 m* 


of wood was pyrolyzing. 


Experiments were performed with a 0.12 m’ window, with and without a 
combustible ceiling in the compartment. Results from sampling inside the compartment 
and along the length of the hallway were compared for the two cases. Sampling inside 
the compartment was performed in the back (0.10 m away from the wall opposite the 
window) and in the front (0.10 m away from the wall containing the window). Both 


sampling probes were 0.10 m below the ceiling and 0.23 m from the side wall. 


The sampling in the hallway was performed 0.46 m, 1.83 m and 3.20 m down the 
hallway from the window. The sampling was done 0.05 m below the ceiling at an equal 
distance from both side walls. The sampling locations and wood ceiling are shown in 


FiggG7/) 


Sampling Sampling 
Location in Location in 
Back of Front of 
Compartment Compartment 


Douglas Fir Plywood Hallway Sampling Locations 


Figure C.7 A top view of the facility showing the orientation of the wood in the upper 
layer of the compartment and the gas sampling locations. 
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C.3.3 Results and Discussion 


The mass loss rate of the wood ceiling was experimentally measured and compared 
with the estimation of the wood mass loss rate found using the oxygen depletion method 
as determined by Pitts, et al. (1994). The burning of the wood during flashover, shown in 
Fig. C.8, is seen to be divided into three regions: (1) gasification stage, (2) transition 
stage and (3) char stage. The first stage of gasification (volatiles escape) is indicated by 
the first peak of 13 g/s m” seen in Fig. C.8. During this stage, as noted by Pitts, et al. 
(1994), gases are formed directly from the wood. Char began forming on the wood 
during the transition stage causing a drop in the mass loss rate to an average value of 
8 g/s m’. Eventually, enough energy was available in the system (i.e. high enough 
temperatures) to pyrolyze the char, and the mass loss rate was measured to rise to a peak 
value of 17 g/s m’. The mass loss rate of 10-15 g/s m” determined by Pitts, et al. (1994) 
is similar to these measured values. As previously noted by Pitts, et al. (1994), the 
pyrolysis of the wood is an endothermic process causing the upper-layer temperatures to 
be lower when a wood ceiling was present compared with results from a compartment 


with no combustible surfaces. 


The species concentrations and the temperatures within the facility having a non- 
combustible ceiling are shown in Fig. C.9. The variations in species concentration and 
temperature for a compartment containing a wooden ceiling are shown in Fig. C.10. The 
CO and UHC concentrations in the experiments with a non-combustible ceiling were 
lower in the rear of the compartment compared with the levels in the front. This was 
attributed to the gases in the rear of the compartment having a longer residence time in 
the compartment. The species concentrations were fairly uniform in the upper-layer of 
the compartment with the wood ceiling. This uniformity is a result of the burning of the 


wood dominating the production of species in the upper-layer. The in-compartment 
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Wood Mass Loss Rate, [g/s m2] 
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Figure C.8 The mass loss rate of the wooden ceiling as a function of time. The sampling 
locations of the five different experiments: Experimental sampling in the 
compartment: —Rear, ----Front, and down the hallway: ¢°0.46 m, -*- 1.83 m, 


-°e- 3.20 m. —Average. 
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Figure C.9 The normalized species concentrations (dry) and temperatures within the 
compartment and hallway for n-hexane fire experiments with a 0.12 m 
opening. Symbols and normalization values: @-CO (15.0%), H-UHC (10%), 
A-CO, (20%), W-0,(10.5%), and O-Temperature. 
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Figure C.10 The normalized species concentrations (dry) and temperatures within the 
compartment and hallway for n-hexane fire experiments containing a wood 
ceiling with a 0.12 m? opening . Symbols and normalization values: @-CO 
(15.0%), HI-UHC (10%), A-CO, (20%), W-O0,(10.5%), and O-Temperature. 
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levels of CO were found to drastically increase to 13.6% in the rear and 12.6% in the 
front compared to 3.1% and 4.1%, respectively, in a compartment without wood. The 
CO concentration produced in the combustible ceiling experiments are consistent with 


those measured in the experiments reported by Pitts, et al. (1994). 


As the exhaust gases exited the compartment into the hallway, there was a 
competition between CO and UHC for the entrained air. The competition is explained 


through the application of the ‘law of mass action’ 
RR = K juei| Cjuet | Or] (C.2) 


where RR is the rate of reaction, kg,.; is the reaction rate constant for the fuel and [Cg,,.] is 
the concentration of the fuel. In the non-combustible ceiling case, both the UHC 
concentration and the reaction rate constant are higher than those of CO resulting in a 
higher oxidation rate of UHC. This is evident in Fig. C.9 by the faster reduction in UHC 
relative to CO. When a wooden ceiling was present in the compartment, the 
concentration of CO was significantly higher than the UHC. This is believed to result in 
a higher CO oxidation rate of reaction compared with that of UHC. This is evident in the 
first part of the hallway as seen in Fig. C.9, by the CO being reduced more rapidly than 
UHC. As the concentrations of CO and UHC became equivalent, the reduction of UHC 
occurred at a higher rate. This was attributed to the competition between the oxidation of 


the two species becoming more a function of the reaction rate constant. 


By the exit of the hallway, the highest in-compartment levels of CO had been 
reduced by 86% in the non-combustible ceiling case and by 79% in the wood ceiling 
case. The UHC exhibited a similar trend with a measured reduction of 96% in the non- 
combustible ceiling case and 82% in the wood ceiling case. The presence of the wood 
ceiling in the compartment led to higher concentrations of CO and UHC entering the 


hallway. The air entrained into the gases flowing in the hallway was, therefore, not 
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sufficient to reduce the CO and UHC as effectively as that measured in experiments with 


no combustible ceiling. 


C.3.4 Summary 


The presence of a Douglas fir plywood ceiling in the compartment had a dramatic 
effect on the level of CO in both the compartment and the hallway, while the UHC 
showed a slight change in their in-compartment levels but a dramatic change in the 
hallway. The CO concentration in the compartment was seen to increase to levels in 
excess of 12%. This is 3-5 times higher (depending on the sampling location) than that 
measured in a compartment without combustible surfaces. The oxidation of the 
compartment fires gases in the hallway in the two cases is also different, and was 
explained through the application of the “law of mass action”. The high concentrations 
of CO generated in the compartment caused the CO to initially be oxidized at a faster rate 
(compared with the UHC) until the concentrations became nearly equivalent. Afterward, 
the UHC were oxidized at faster rate. For the non-combustible ceiling case, the UHC 
levels were always measured to be reduced at a faster rate than CO. Most important, the 
presence of the wood ceiling allowed the CO to leave the hallway at 2.5% which is lethal 


after approximately 2 minutes of exposure. 
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APPENDIX D. 


CHEMICAL KINETIC MODELING OF THE HALLWAY 
GASES 


D.1 INTRODUCTION 


A plug flow reactor model for the oxidation and transport of compartment fire gases 
entering and traveling down the hallway is developed in this section. The equations 
describing the fluid mechanic, thermodynamic and chemical behavior of the species in a 
plug flow reactor are derived from basic principles. Application of assumptions reduces 
the set of equations to only the equation for species conservation. Finally, 
transformations for converting spatial experimental data to time data, and vice versa, is 


given. 


D.2 DERIVATION OF GOVERNING EQUATIONS 


In the following derivations, consider a plug flow reactor in steady-state and steady- 


flow, see Fig. D.1. One-dimensional (1-D) flow is considered with uniform properties 


et eae eet Ge A 


— — i — 


Figure D.1 The plug flow reactor being used to model the hallway gases. 
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present at each cross section. No mixing is assumed to occur in the axial direction (x- 
direction), and the flow is assumed to be frictionless. The gases within the reactor are 


assumed to behave as ideal gases. 


Figure D.2 shows a control volume of gases in the plug flow reactor of length dx. 
Summing the mass flowing into and out of the control volume, the law of mass 


conservation yields: 


oud dIPHA) a oud =0 : (D.1) 
ax 
and reduces to 
d 
—\(pud)=0 2 
Fat ) (D.2) 


The momentum of the gases entering and exiting the control volume shown in Fig. D.2 is 


equivalent to the forces acting upon the control volume. Assuming no body forces, 


dx aie | (PA) 

cai jal -— PA + dx 

PA ———> 
> d(u (upA 
u (upA) ———_ u (upA) + eta) dx 
MES 

ms | 
BS eal 
Ss 


Figure D.2 The control volume showing the fluxes of momentum. 
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d\u’pA A 
Ay gi ec ey apr 5 
dx dx 
which is equivalent to 
d\u’pA 
(up ) _ pd4_a(PA) D4) 


ax axe ids 


Expanding the first term on the left hand side and the second term on the right hand side 


using the chain rule, 


ior tao i aug Sala Cl (D.5) 
ds ee te 


u 
Applying continuity, Eqn. (D.2), and a little algebra transforms Eqn. (D.5) into 
up— = -— (D.6) 


To derive the energy and species equations, consider the control volume shown in 


Pigtins. 
q’’ Q dx 
a 
ve | 
PHaRS t d(m (u? /2+ h.)) 

wh eee edi (u? /2+ h,) + —————— dx 
m (u2/2 + h,}———> | 
| 
| 
| 
\ | 


Figure D.3 The control volume showing the fluxes of energy and species. 
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By summing the energy into, out of and stored within the control volume 


e d(pu(h, +22/2)) 


- dx — pud(h, +u?/2)=-q"'Qdx (D7) 


pud(h, + u?/2 


After canceling terms one and three on the left side and using the chain rule to expand the 
second term on the same side, the application of continuity reduces Eqn. (D.7) to 
waiter l2) =—g"'Q (D.8) 
dx 
Similarly, the summation of the species into, out of and stored within the control volume 
yields the following relation 


d (pudmf,) 
ax 


puAdmf; + dx — puAmf, = o ; MW, Adx (D.9) 


Through the cancellation of terms, the expansion of the second term on the left side , and 
the application of continuity, Eqn. (D.9) is reduced to 

anf; _ © ,MW, (D.10) 
dk pu 


To model the 1-D system of equations, Eqns. (D.2), (D.6), (D.8) and (D.10), a 
computational fluid dynamics (CFD) code must be written. Instead, further manipulation 
was performed on Eqn. (D.6) to reduce the complexity of the problem. If the plug flow 
reactor is consider to have a constant cross sectional area, then the (pu) term in equation 


(D.6) is constant. This allows Eqn. (D.6) to be integrated, 
1 1 
(pu) | du=—| dP, (D.11) 


yielding the following expression 


(pu)(u,—u,) = P, - P, (D.12) 
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Applying continuity with the constant cross sectional area assumption invoked, 


Egn. (D.11) is transformed to 


Aitat (D.13) 


wi 

[(Z/e,)-U/e.) | 

Equation (D.13) represents the equation of the negatively sloped Rayleigh line. If q’’ is 
considered to be zero, then Eqn. (D.8) can be integrated and yields the following 


relationship 
h, -h, =|(u2/2)-(u7/2)| (D.14) 


If the velocity is put in terms of mass flow again and the mass flow from Eqn. (D.13) is 
substituted into the equation, the Eqn. (D.14) becomes 


P,-P. 
p= |(W/ex)- 7/0.) (D.15) 
the Hugoniot equation. For a deflagration wave, the density variation is responsible for 
the majority of the rise in enthalpy and the pressure difference is considered small. By 


assuming no pressure drop, Eqn. (D.6) is no longer a function of x. 


Though the energy equation could be solved to determine the temperature along the 
reactor, a temperature will be specified along the length of the reactor from experimental 


results. Thus, the energy equation need not be solved. 


From this additional analysis and consideration of the problem, the number of 
fundamental equations is reduced down to the species conservation Eqn. (D.10). The 
number of ordinary differential equations which must be solved is actually equivalent to 


the number of species, j, considered. 
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D.3 METHODOLOGY OF MODELING 


D.3.1 Transformation of Experimental Results to a Temporal Coordinate 
System 


Figure D.4 shows the location of the sampling points in the hallway and the spatial 
transformation of the L-shaped path. 


Compartment 


Figure D.4 The bulk flow path of the compartment fire gases. The letters denote 
locations where temperature and air entrainment measurements were 
performed. 


Assuming that the gases behave as ideal gases, the density can be determined from 


P 
= — D.16 
PH oF (D.16) 
A discretized form equation of the velocity, 
Ax 
u=—, Ad 
i¥ (D.17) 


defines the gas velocity between two points in the hallway. Since velocity measurements 


are not made, the conservation of mass is applied 
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zal) (D.18) 
p 


to determine the velocity in terms of known values. From the continuity equation applied 
in the plug flow reactor derivation, the value of (m/ A) in Eqn. (D.18) is constant. A 
value for m and A can be arbitrarily chosen at any location along the length of the 
hallway. The m is, therefore, the mass flow rate exiting the compartment while A is the 


area of the opening. Substituting the right hand side of Eqn. (D.16) for p into 


Eqn. (D.18), an expression containing previously measured values is generated 


_ (m/ A)RT 
To ee (D.19) 


Solving Eqn. (D.17) for At and substituting Eqn. (D.19) in for the velocity term, a 
discretized expression for the transformation of the spatial coordinates into temporal 


coordinates 


PAx 
At =7= D.20 
(m1 /A )RT ( 
is formed. The value for P is equivalent to the ambient pressure, while Ax is the distance 
between sampling positions. The ideal gas constant is divided by the molecular weight of 
the upper-layer gases (assumed to be equivalent to that of air) to determine R, and T is the 


temperature of the gases in the upper-layer of the hallway. 


To determine the temperature distribution in the plug flow reactor, the reactor is 
divided into 10 sections with the center of each section containing a sampling point, as 
seen in Fig. D.5. Each section takes on the temperature which is measured at each point 
and the velocity is determined from the temperature using Eqn. (D.19). Based on the 
velocity and Ax of the section, a At for each section is calculated using Eqn. (D.20). 
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Lise | 
iaRREBBR 


x=0 


Figure D.5 The location of the data points within the plug flow reactor. 


The time at each point after gases enter the reactor was determined through linear 
interpolation between values of x and t. The air entrainment data, taken at the locations 
where the temperatures were measured, was transformed into the temporal coordinate 
system using the same analysis performed to convert the spatial temperature data to a 


temporal temperature data. 


D.3.2 Transformation of Plug Flow Reactor Results to Spatial Coordinate 
System 


To compare experimental results with SENKIN results, the data from SENKIN is 
transformed from a temporal coordinate system to a spatial coordinate system. To do 


this, consider the differential form of the velocity, 


dx 


=— D.21 
ory (D.21) 


Converting the velocity to mass flow and integrating the left side from 0 to t and the right 
side from 0 to x , Eqn. (D.21) becomes 


[ide =(mm/ A) (1/p)at (D.22) 
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By assuming a steady state process, the cross sectional area of the plug flow reactor and 
the mass flow rate are constant with time. The cross sectional area of the plug flow 
reactor and the mass flow can be taken outside of the integrand. Applying the ideal gas 
law to Eqn. (D.22) to rewrite the density in terms of temperature, the following 
formulation is obtained: 


(eee bf a f, Tae (D.23) 


Note that Eqn. (D.20) is simply the integrated form of Eqn. (D.23), having assumed 
constant temperature across the integrated area. Assuming constant pressure and 
integrating the left hand side of the equation, the following transformation from spatial to 


temporal coordinates is derived 


ae 


NEES 7 (D.24) 
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APPENDIX E. 


EXPERIMENTAL ERROR ANALYSIS 


E.1 INTRODUCTION 


The experimental error of the measured and calculated values in the research is 
determined in the following sections. The error into parameters measured directly is 
addressed first. Error in the parameters calculated using the measured values is then 


calculated. 


E.2 CALCULATION OF ERROR 


The calculation of error for a parameter where all the variables are assumed to be 


linear dependent with one another is 
e? =e, +e, +6, +82... (FE) 


Equation (E.1) is used to calculate the error for measured quantities where the 


dependence of each source of error is proportional to the other sources of error. 


When determining the error for parameters calculated from equations, a more 


fundamental method is applied: 


2 2 2 
w-( 2 (2m, [2 )+- (E22) 
2 


where w is the variation in a particular variable, y is the parameter being calculated, and 
the x’s are the variables in the equation. For proof sake, consider a function y whose 


variables are linearly dependent on one another, y=x,x,x3. Taking the appropriate partial 
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derivatives, and then dividing both sides of Eqn. (E.2) by y, the following formulation is 


2), OS EY Any SY 
Seley « 
yy x} x2 x3 


C2, 2 2 
+ Saleh (E.4) 


produced 


or, 


Equation (E.4) is, as expected, equivalent to (E.1). 


E.3 ERROR IN PARAMETERS 
E.3.1 Measured Species Concentrations 


The measurement error of CO, CO,, O, and UHC concentrations was attributed to 


four sources: 
e analyzer repeatability, E,p, 
e calibration gas, Eg, 
e calibration uncertainty, ecy, and 
e experimental repeatability, exp. 


The CO, CO, and O, concentrations were all measured dry while the UHC were 


measured wet. 


The repeatability of the analyzers was noted by the manufactures as being +1% 
(€,p=0.01). The maximum error in the calibration gas concentration, determined by the 
supplier using a gas chromatograph, was +1.5% (€¢g=0.015). The last two sources of 


error depended upon whether the gases were being sampled in the hallway (or 
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Table E.1. The uncertainty induced in the measurement of the species concentrations. 


TmSpecies | “eczat tf es aud eee 
PUHC | 0.01 | 0.025 | 0.036 _| 


compartment) or in the exhaust duct, and the species which was being measured. 
Tabulated values for these errors are given in Table E.1. Utilizing these values for the 
error, the error when measuring species in the exhaust duct, €,,, and in the 


hallway/compartment, €,/,, was determined by applying Eqn. (E.3). 
E.3.2 Wet Species Concentrations 


Assuming the CO, and H,O are present in a stoichiometric proportions, the wet 


concentrations of CO, CO, and O, can be calculated from 


X id 


X, y= (E.5) 
ee TRTTG oe 
Gottuk estimated the assumption of the CO, and H,O being in stoichiometric proportions 
as being valid to within +6%. Utilizing the tabulated values for the error in the dry 
concentration measurement, the error in the calculation of the wet concentrations is 


determined, see Table E.2. 


E.3.3 Gas Temperature 


Gas temperatures were measured in both the compartment and the hallway. The 
error in the temperature measurements was due to the aspiration process (utilized to 


shield the thermocouples from radiation) and the repeatability of the experiment. 
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Table E.2. The uncertainty in the calculation of the wet concentrations of CO, CO, and 
O,. 


Ss ad pee re 


Inside the compartment, the error in the aspiration was approximated to be 1.0%. 
The repeatability of the experimental values was determined to be approximately 3% 
through comparison of two experiments run under the same conditions. The result was 


an error of +3.2%. 


In the hallway, the temperatures in the hallway are more transient compared to the 
compartment. The error due to aspiration was determined to again be 3%. Through the 
comparison of the same experiments utilized in the compartment error analysis, the 
experimental repeatability of the temperatures was within 13%. The error in the hallway 
temperatures was calculated to be +13.2%. 


E.3.4 Fuel Weight 


The error induced on the weight of the fuel was attributed to the thermal drift in the 
load cell, the resolution of the load cell and the A/D conversion. The thermal drift was 
considered to be the most significant of the three sources, and was estimated by 
comparing the weight of the fuel pan before and after the experiment. The error was 
found to be approximately 1%. With the load cell being accurate to within 1 gram, the 
error in the resolution was considered to be negligible as was the A/D conversion. The 


error in the fuel weight was, therefore, +1%. 
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E.3.5 Air Mass Flow Rate into Compartment 


The measurement of air mass flow rate into the compartment was determined to have 
three sources of error: the temperature measurement (used to determine the density) (1%), 
the location of the velocity probe relative to the location average velocity (4%), and the 
use of the average velocity instead of an integrated velocity profile (1%). Based on these 


sources of error, the error in the measurement was calculated to be +4.2%. 


E.3.6 Global Equivalence Ratio 


The global equivalence ratio is defined as 


b= 7 ee (E.6) 
Bee 
Moir ‘¢f 
Since the stoichiometric fuel to air ratio is a constant, the error in the fuel volatilization 


rate and the air entrainment rate into the compartment need only be considered. The error 


in the fuel volatilization rate was determined from 
Wy SEN ii, (E.7) 


which corresponds to the maximum variation in the numerical derivative of the fuel 
weight performed over a +10, h, second band around the point. The maximum value of 
im, of 9.5 x 10° was found to occur just prior to post-flashover. Using a value of h=10, 
the maximum variation in the fuel volatilization rate is expected to be 1.58 x 10°. Using 
the fuel volatilization rate at that instant where the maximum second derivative of the fuel 


volatilization rate was found, the error was determined to be +6%. 


Using the error in the air entrainment rate into the compartment, +4.2%, in 
conjunction with the error in the fuel volatilization rate, the error in the global 


equivalence ratio was calculated to be +7.3%. 
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Table E.3. The uncertainty in the post-hallway yields for the different species measured 
in the study. 


E.3.7 Species Yields 
The sources of error in the species yield 


X; MW he, 
aa (E.8) 
Th, 


calculation were the error of all variables in Eqn. (E.8) and the time shift in the data to 
align the gas measurement data with the instantaneous measurement of the fuel weight. 
Assuming the shift results in at most 1% error, the molecular weight is within 0.5% of 
that of air and that 7,, is within 10% of the actual molar flow rate within the duct; the 


error in the yields for each species is given in Table E.3. 
E.3.8 Fire Size 


The error in the ideal fire size is dependent upon the fuel volatilization rate, since the 
error in the heat of combustion multiplying the fuel volatilization rate is comparatively 


negligible. The error in the fire size is, therefore, calculated to be +6%. 
E.3.9 Nondimensional Hallway Upper-Layer Depth 


The sources of error in the determination of the nondimensional hallway upper-layer 
depth, y=5/z, are from the measurement of the upper-layer within the hallway, 6, and the 
measurement of the distance between the bottom of the opening and the ceiling, z. The 


error in the measurement in the upper-layer depth was found to be 12.5%, while the value 
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of z only varied less than 0.5%. The analysis determined the error the nondimensional 


upper-layer depth to be +12.5% 
E.3.10 Entrainment Function 


With the entrainment function defined as the ratio between the entraining velocity 
and the velocity of the gases entering the hallway. The inlet velocity was known to 
within 10%. The velocity of the entraining air depended on the mass flow of air 
entrained, the density of the entraining air, and the wetted area of the jet. The largest 
errors were calculated for the mass flow of air entrained (15%) and the wetted area of the 
jet (19%), resulting in a 24% error in the entrained velocity. The entrainment function 


error is, therefore, +26%. 
E.3.11 Richardson Number 


The sources of error in the Ri is attributed to the density of the gases entering the 
hallway and the entraining gases (10%), the velocity of the gases entering the hallway 
(10%) and the hydraulic diameter of the opening (15%). The error in the Ri number is 


expected to be no greater than +25%. 
E.3.12 Ignition Index 


The majority of the error in the ignition index is from the estimation of H, from the 
CO, and the assumption that all the UHC are ethylene. The other sources of error are the 
measurement of CO and temperature. The error for the H, was estimated by Gottuk 
(1992a) as being 10%, while the UHC error was taken from Table E.1 to be 10.2%. The 
error in the CO measurement was 11.9%, but the temperature is accurate to within 3%. 


The application of these errors to Eqn. (E.1) produces a 18.8% error in the ignition index. 
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E.4 SUMMARY 


Table E.4 lists the errors associated with the measured and calculated parameters. 
Table E.4 The errors of all the parameters measured and calculated in the study. 


ERROR, (+) 


PARAMETER 


[CO]ary (Ex. Duct/ Hall-Comp.) 4.5% / 10.3% 
[CO2]é, (Ex. Duct/ Hall-Comp.) 4.4% / 7.4% 
[O>]ary (Ex. Duct/ Hall-Comp.) 2% / 5.3% 


Se 
CO Yield (Post-Hallway) 

CO, Yield (Post-Hallway 

: 

Ri 


ele 
Ignition Index, I.I 18.8% 
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